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Abstract 

We introduce a set of operations that we call crystal operations on matrices with entries either in 
{0, 1} or in N. There are horizontal and vertical forms of these operations, which commute with each 
other, and they give rise to two different structures of a crystal graph of type A on these sets of 
matrices. They provide a new perspective on many aspects of the RSK correspondence and its dual, 
and related constructions. Under a straightforward encoding of semistandard tableaux by matrices, the 
operations in one direction correspond to crystal operations applied to tableaux, while the operations 
in the other direction correspond to individual moves occurring during a jeu de taquin slide. For the 
(dual) RSK correspondence, or its variant the Burge correspondence, a matrix M can be transformed by 
horizontal respectively vertical crystal operations into each of the matrices encoding the tableaux of the 
pair associated to M, and the inverse of this decomposition can be computed using crystal operations 
too. This decomposition can also be interpreted as computing Robinson's correspondence, as well as the 
Robinson-Schensted correspondence for pictures. Crystal operations shed new light on the method of 
growth diagrams for describing the RSK and related correspondences: organising the crystal operations 
in a particular way to determine the decomposition of matrices, one finds growth diagrams as a method 
of computation, and their local rules can be deduced from the definition of crystal operations. The 
Schiitzenberger involution and its relation to the other correspondences arise naturally in this context. 
Finally we define a version of Greene's poset invariant for both of the types of matrices considered, and 
show directly that crystal operations leave it unchanged, so that for such questions in the setting of 
matrices they can take play the role that elementary Knuth transformations play for words. 

§0. Introduction. 

The Robinson-Schensted correspondence between permutations and pairs of standard Young tableaux, 
and its generalisation by Knuth to matrices and semistandard Young tableaux (the RSK correspondence) 
are intriguing not only because of their many surprising combinatorial properties, but also by the great 
variety in ways in which they can be defined. The oldest construction by Robinson was (rather cryptically) 
defined in terms of transformations of words by "raising" operations. The construction by Schensted 
uses the well known insertion of the letters of a word into a Young tableau. While keeping this insertion 
procedure, Knuth generalised the input data to matrices with entries in N or in {0, 1}. He also introduced 
a "graph theoretical viewpoint" (which could also be called poset theoretical, as the graph in question is 
the Hasse diagram of a finite partially ordered set) as an alternative construction to explain the symmetry 
of the correspondence; a different visualisation of this construction is presented in the "geometrical form" 
of the correspondence by Viennot, and in Fulton's "matrix-ball" construction. A very different method of 
describing the correspondence can be given using the game of "jeu de taquin" introduced by Lascoux and 
Schiitzenberger. Finally a construction for the RSK correspondence using "growth diagrams" was given 
by Fomin; it gives a description of the correspondence along the same lines as Knuth's graph theoretical 
viewpoint and its variants, but it has the great advantage of avoiding all iterative modification of data, 
and computes the tableaux directly by a double induction along the rows and columns of the matrix. 

The fact that these very diverse constructions all define essentially the same correspondence (or at 
least correspondences that are can be expressed in terms of each other in precise ways) can be shown using 
several notions that establish bridges between them. For instance, to show that the "rectification" process 
using jeu de taquin gives a well defined map that coincides with the one defined by Schensted insertion, 
requires (in the original approach) the additional consideration of an invariant quantity for jeu de taquin 
that is a special case of Greene's invariant for finite posets, and of a set of elementary transformations 
of words introduced by Knuth. A generalisation of the RSK correspondence from matrices to "pictures" 
was defined by Zelevinsky, which reinforces the link with Littlewood-Richardson tableaux already present 
in the work of Robinson; it allows the correspondences considered by Robinson, Schensted, and Knuth 
to be viewed as derived from a single general correspondence. The current author has shown that that 
correspondence can alternatively be described using (two forms of) jeu de taquin for pictures instead 
of an insertion process, and that in this approach the use of Greene's invariant and elementary Knuth 
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operations can be avoided. A drawback of this point of view is that the complications of the Littlcwood- 
Richardson rule are built into the notion of pictures itself; for understanding that rule we have also 
given a description that is simpler (at the price of loosing some symmetry) , where semistandard tableaux 
replace pictures, and "crystal" raising and lowering operations replace one form of jeu de taquin, so that 
Robinson's correspondence is completely described in terms of jeu de taquin and crystal operations. 

In this paper we shall introduce a new construction, which gives rise to correspondences that may 
be considered as forms of the RSK correspondence (and variants of it). Its interest lies not so much in 
providing yet another computational method for that correspondence, as in giving a very simple set of 
rules that implicitly define it, and which can be applied in sufficiently flexible ways to establish links with 
nearly all known constructions and notions related to it. Simpler even than semistandard tableaux, our 
basic objects are matrices with entries in N or in {0, 1}, and the basic operations considered just move 
units from one entry to an adjacent entry. As the definition of those operations is inspired by crystal 
operations on words or tableaux, we call them crystal operations on matrices. 

Focussing on small transformations, in terms of which the main correspondences arise only implicitly 
and by a non-deterministic procedure, our approach is similar to that of jeu de taquin, and to some extent 
that of elementary Knuth transformations. By comparison our moves are even smaller, they reflect the 
symmetry of the RSK correspondence, and they can be more easily related to the constructions of that 
correspondence by Schensted insertion or by growth diagrams. Since the objects acted upon are just 
matrices, which by themselves hardly impose any constraints at all, the structure of our construction 
comes entirely from the rules that determine when the transfer of a unit between two entries is allowed. 
Those rules, given in definitions 1.3.1 and 1.4.1 below, may seem somewhat strange and arbitrary; how- 
ever, we propose to show in this paper is that in many different settings they do precisely the right thing 
to allow interesting constructions. One important motivating principle is to view matrices as encoding 
semistandard tableaux, by recording the weights of their individual rows or columns; this interpretation 
will reappear all the time. All the same it is important that we are not forced to take this point of view: 
sometimes it is clearest to consider matrices just as matrices. 

While the above might suggest that we introduced crystal operations in an attempt to find a unifying 
view to the numerous approaches to the RSK correspondence, this paper in fact originated as a sequel 
to [vLee5] , at the end of which paper we showed how two combinatorial expressions for the scalar product 
of two skew Schur functions, both equivalent to Zelevinsky's generalisation of the Littlewood-Richardson 
rule, can be derived by applying cancellations to two corresponding expressions for these scalar products 
as huge alternating sums. We were led to define crystal operations in an attempt to organise those 
cancellations in such a way that they would respect the symmetry of those expressions with respect to 
rows and columns. We shall remain faithful to this original motivation, by taking that result as a starting 
point for our paper; it will serve as motivation for the precise form of the definition crystal operations on 
matrices. That result, and the Littlewood-Richardson rule, do not however play any role in our further 
development, so the reader may prefer to take the definitions of crystal operations as a starting point, 
and pick up our discussion from there. 

This paper is a rather long one, even by the author's standards, but the reason is not that our 
constructions are complicated or that we require difficult proofs in order to justify them. Rather, it is the 
large number of known correspondences and constructions for which we wish to illustrate the relation with 
crystal operations that accounts for much of the length of the paper, and the fact that we wish to describe 
those relations precisely rather than in vague terms. For instance, we arrive fairly easily at our central 
theorem 3.1.3, which establishes the existence of bijective correspondences with the characteristics of 
the RSK correspondence and its dual; however, a considerable additional analysis is required to identify 
these bijections precisely in terms of known correspondences, and to prove the relation found. Such 
detailed analyses require some hard work, but there are rewards as well, since quite often the results have 
some surprising aspects; for instance the correspondences of the mentioned theorem turn out to be more 
naturally described using column insertion than using row insertion, and in particular we find for integral 
matrices the Burge correspondence rather than the RSK correspondence. We do eventually find a way 
in which the RSK correspondence arises directly from crystal operations, in proposition 4.4.4, but this is 
only after exploring various different possibilities of constructing growth diagrams. 

Our paper is organised as follows. We begin directly below by recalling from [vLee5] some basic 
notations that will be used throughout. In §1 we introduce, first for matrices with entries in {0,1} 
and then for those with entries in N, crystal operations and the fundamental notions related to them, 
and we prove the commutation of horizontal and vertical operations, which will be our main technical 
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tool. In §2 we mention a number of properties of crystal graphs, which is the structure one obtains by 
considering only vertical or only horizontal operations; in this section we also detail the relation between 
crystal operations and jeu de taquin. In §3 we start considering double crystals, the structure obtained 
by considering both vertical and horizontal crystal operations. Here we construct our central bijective 
correspondence, which amounts to a decomposition of every double crystal as a Cartesian product of two 
simple crystals determined by one same partition, and we determine how this decomposition is related to 
known Knuth correspondences. In §4 we present the most novel aspect of the theory of crystal operations 
on matrices, namely the way in which the rules for such operations lead to a method of computing 
the decomposition of theorem 3.1.3 using growth diagrams. The use of growth diagrams to compute 
Knuth correspondences is well known of course, but here the fact that such a growth diagram exists, 
and the local rule that this involves, both follow just from elementary properties of crystal operations, 
without even requiring enumerative results about partitions. In §5 we study the relation between crystal 
operations and the equivalents in terms of matrices of increasing and decreasing subsequences, and more 
generally of Greene's partition-valued invariant for finite posets. Finally, in §6 we provide the proofs of 
some results, which were omitted in text of the preceding sections to avoid distracting too much from 
the main discussion. (However for all our central results the proofs are quite simple and direct, and we 
have deemed it more instructive to give them right away in the main text.) 

0.1. Notations. 

We briefly review those notations from [vLee5] that will be used in the current paper. We shall use 
the Iverson symbol, the notation [ condition ] designating the value 1 if the Boolean condition holds 
and otherwise. For n G N we denote by [n] the set {i G N | i < n} = {0, . . . , n — 1}. The set C of 
compositions consists of the sequences (aj)j e N with q, 6 N for all i, and en = for all but finitely 
many i; it may be thought of as UneN ^™ w here each N" is considered as a subset of N" +1 by extension 
of its vectors by an entry 0. Any a G C is said to be a composition of the number \a\ = XieN a i- We 
put C [2] = { a £ C | Vi: £ [2] }; its elements are called binary compositions. The set V C C of partitions 
consists of compositions that are weakly decreasing sequences. The operators '+' and ' — ', when applied 
to compositions or partitions, denote componentwise addition respectively subtraction. 

The diagram of A £ ?, which is a finite order ideal of N 2 , is denoted by [A], and the conjugate 
partition of A by A*. For K,AeP the symbol A/re is only used when [k] C [A] and is called a skew shape; 
its diagram [X/k] is the set theoretic difference [A] — [k], and we write |A//i| = |A| — \fj,\. For a, (3 G C 
the relation a -j— (3 is defined to hold whenever one has (3i + \ < Qj < for all i G N; this means that 
a, (3 G V, that [a] C \(3], and that [(3 /a] has at most one square in any column. When /j, A, the skew 
shape A//x is called a horizontal strip. If /i* v- A* holds, we call A//i a vertical strip and write /j, A; 
this condition amounts to n, A G V and A-fi e C l2] (so [X/n] has at most one square in any row). 

A semistandard tableau T of shape X/k (written T G SST(A//t)) is a sequence (A^)i 6 N of partitions 
starting at k and ultimately stabilising at A, of which successive members differ by horizontal strips: 

^_ f or a n { e N. The weight wt(T) of T is the composition (|A( J+1 ) /A^ |) ieN . Although we shall 

work mostly with such tableaux, there will be situations where it is more natural to consider sequences in 
which the relation between successive members is reversed (A^ — A^ +1 ') or transposed (A^ A^ +1 ^), 
or both (AW —r A^ +1 '); such sequences will be called reverse and/or transpose semistandard tableaux. 
The weight of transpose semistandard tableaux is then defined by the same expression as that of ordinary 
ones, while for their reverse counterparts it is the composition (|A ( ^/^ < ' i+1 ' ) l)ieN- 

The set M is the matrix counterpart of C: it consists of matrices M indexed by pairs (i, j) G N 2 , 
with entries in N of which only finitely many are nonzero (note that rows and columns are indexed 
starting from 0). It may be thought of as the union of all sets of finite matrices with entries in N, 
where smaller matrices are identified with larger ones obtained by extending them with entries 0. The 
set of such matrices with entries restricted to [2] = {0,1} will be denoted by M l2] ; these are called 
binary matrices. For matrices M G A4, we shall denote by Mj its row i, which is (Mij)j e w G C, while 
Mj = (Mj ; j)j e N G C denotes its column j. We denote by row(M) = (|Mj|)j e N the composition formed 
by the row sums of M, and by col(M) = (|M- |)j 6 N the composition formed by its column sums, and we 
define M a . fj = {M&M\ row(Af) = a, col(M) = (3 } and M 1 ^ = M a ,p n M l2] . 

In the remainder of our paper we shall treat M [2] and M as analogous but separate universes, in 
other words we shall never consider a binary matrix as an integral matrix whose entries happen to be < 1 
or vice versa; this will allow us to use the same notation for analogous constructions in the binary and 
integral cases, even though their definition for the integral case is not an extension of the binary one. 
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§1. Crystal operations on matrices. 

The motivation and starting point for this paper are formed by a number of expressions for the scalar 
product between two Schur functions in terms over enumerations of matrices, which were described 
in [vLce5]. To present them, we first recall the way tableaux were encoded by matrices in that paper. 

1.1. Encoding of tableaux by matrices. 

A semistandard tableau T = (AW)i £ N of shape X/k can be displayed by drawing the diagram [X/k] in 
which the squares of each strip [A^ +1 ^ /A^] are filled with entries i. Since the columns of such a display 
are strictly increasing and the rows weakly increasing, such a display is uniquely determined by its shape 
plus one of the following two informations: (1) for each column Cj the set of entries of Cj , or (2) for each 
row Ri the multiset of entries of Ri. Each of those informations can be recorded in a matrix: the binary 
matrix M G M [2] in which Mij G {0, 1} indicates the absence or presence of an entry i in column Cj 
of the display of T will be called the binary encoding of T, while the integral matrix N € A4 in which 
Ni j gives the number of entries j in row of the display of T will be called the integral encoding of T. 
In terms of the shapes X^ these matrices can be given directly by M t = (A( 4+1 ))* - (A«)' for all i and 
Nj = (A^' +1 )) - (A^) for all j, cf. [vLee5, definition 1.2.3]. Note that the columns M| of the binary 
encoding correspond to the columns Cj, and the rows Ni of the integral encoding to the rows Ri. While 
this facilitates the interpretation of the matrices, it will often lead to an interchange of rows and columns 
between the binary and integral cases; for instance from the binary encoding M the weight wt(T) can 
be read off as row(M), while in terms of the integral encoding N it is col(A). Here is an example of 
the display of a semistandard tableau T of shape (9, 8, 5, 5, 3)/(4, 1) and weight (2, 3, 3, 2, 4, 4, 7), with its 
binary and integral encodings M and N, which will be used in examples throughout this paper: 





















f° 


i 








1 







































2 


4 


■5 


5 




1 


i 


1 



















( l 





1 





1 


2 









1 


3 


4 


G 


6 


6 




1 





1 








1 













1 


1 





1 


1 





3 


1 


1 


2 


4 


5 








, M: 





1 





1 

















, N: 





2 


1 





1 


1 





2 


3 


4 


G 


G 
















1 


1 


1 





1 
















1 


1 


1 





2 

2/ 


■5 


6 


6 














1 











1 








1 


1 




Vo 














1 




















Vo 


1 


1 


1 


1 


1 


1 


I 


oJ 

















To reconstruct T from its binary or integral encoding, one needs to know the shape X/k of T, which 
is not recorded in the encoding; since A and k are related by A' — k 1 = col(M) in the binary case and by 
A — k = row(N) in the integral case, it suffices to know one of them. Within the sets Ai l2] and M. of all 
binary respectively integral matrices, each shape X/k defines a subset of matrices that occur as encodings 
of tableaux of that shape: we denote by Tabl [2] (A/fv) C M l2] the set of binary encodings of tableaux 
T € SST(A/k), and by Tabl(A/«) C M. the set of integral encodings of such tableaux. The conditions 
that define such subsets, which we shall call "tableau conditions", can be stated explicitly as follows. 

1.1.1. Proposition. Let X/k be a skew shape. For M G M [2] one has M e Tabl [21 (A/ft) if and only if 
col(M) = A' - k\ and k* + J2 i<k M l e V for all fc G N. For M G M one has M G Tabl(A//c) if and only 
if row(M) = A - k, and (k + J2 3<i Mj) ^(k + J2j<i Mj) for all I G N. 

Proof. This is just a verification that an appropriate tableau encoded by the matrix can be reconstructed 
if and only if the given conditions are satisfied. We have seen that if M G M [2] is the binary encoding of 
some (A«) ieN G SST(A/k), then M l = (X^+^Y - (A«)* for all i, which together with A^ = k implies 
(AW)* = K* + E 4 <fc M * for k e N. A sequence of partitions A^ satisfying this condition exists if and only 
if each value K t + J ^2 i<k Mi is a partition. If so, each condition A^ -j— A^ +1 ^ will be automatically satisfied, 
since it is equivalent to (A^)* (A (l+1) )*, while by construction (X^+^Y - (A w )* = M t 6 C [2] ; therefore 
(A^')i e N will be a semistandard tableau. Moreover col(M) = A* — k* means that k 1 + ^2 i<k Mi — A* for 
sufficiently large k, and therefore that the shape of the semistandard tableau found will be X/k. 

Similarly if M G M. is the integral encoding of some (A^')ieN G SST(A/k), then we have seen that 
M] = (A^ +1 )) - (A^) for all j, which together with A<°) = k implies \W = k + J2j<i Mj for I G N. 
By definition the sequence (A^^)j £ N so defined for a given k and M G M is a semistandard tableau 
if and only if A^ t— A^ +1 ' for all I G N (which implies that all A^ are partitions), in other words 
if and only if (k + J2 0< i M j) ^ (« + Ej<; M j) for al l 1 G N. The value of A (() ultimately becomes 
k + row(M), so the semistandard tableau found will have shape X/k if and only if row(M) = A — k. □ 
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1.2 Commuting cancellations 



Littlewood-Richardson tableaux are semistandard tableaux satisfying some additional conditions, 
and the Littlewood-Richardson rule expresses certain decomposition multiplicities by counting such 
tableaux (details, which are not essential for the current discussion, can be found in [vLee3]). In [vLee5, 
theorems 5.1 and 5.2], a generalised version of that rule is twice stated in terms of matrices, using respec- 
tively binary and integral encodings. A remarkable aspect of these formulations is that the additional 
conditions are independent of the tableau conditions that these matrices must also satisfy, and notably of 
the shape X/k for which they do so; moreover, the form of those additional conditions is quite similar to 
the tableau conditions, but with the roles of rows and columns interchanged. We shall therefore consider 
these conditions separately, and call them "Littlewood-Richardson conditions" . 

1.1.2. Definition. Let v/a be a skew shape. The set LR l2] (v/n) C M [2] is defined by M G LR [21 (zVm) 
if and only if row(M) = v — u, and u + J2j>i Mj G V for all I G N, and the set LR(u/u) C M is defined 
by M G LR(u/u) if and only if col(M) = v~- a, and (u + J2 i<k M i) ^ (P + E J<fe M i) for a11 k <EN. 

Thus for integral matrices, the Littlewood-Richardson conditions for a given skew shape are just the 
tableau conditions for the same shape, but applied to the transpose matrix. For binary matrices, the 
relation is as follows: if M is a finite rectangular binary matrix and M' is obtained from M by a quarter 
turn counterclockwise, then viewing M and M' as elements of M l2] by extension with zeroes, one has 
M G LR [21 (A/k) if and only if M' G Tabl [21 (A/K). Note that rotation by a quarter turn is not a well 
defined operation on M. [21 , but the matrices resulting from the rotation of different finite rectangles that 
contain all nonzero entries of M are all related by the insertion or removal of some initial null rows, and 
such changes do not affect membership of any set Tabl [21 (A/ft) (they just give a shift in the weight of the 
tableaux encoded by the matrices). 

1.2. Commuting cancellations. 

We can now concisely state the expressions mentioned above for the scalar product between two skew 
Schur functions, which were given in [vLee5]. What interests us here is not so much what these expressions 
compute, as the fact that one has different expressions for the same quantity. We shall therefore not recall 
the definition of this scalar product {s\/ K | s^/^), but just note that the theorems mentioned above 
express that value as #(Tabl [21 (A/«) n LR |21 (v/fi)) and as #(Tabl(A/«) flLR(i///i)), respectively (the 
two sets counted encode the same set of tableaux). Those theorems were derived via cancellation from 
equation [vLee5, (50)], which expresses the scalar product as an alternating sum over tableaux. That 
equation involves a symbol s(a, A), combinatorially defined for a G C and A G V with values in { — 1, 0, 1}. 
For our current purposes the following characterisation of this symbol will suffice: in case a is a partition 
one has e(a, A) = [a = A] , and in general if a, a' G C are related by (c^,c^ +1 ) = (o;j+i — l,aj + 1) 
for some i G N, and a'j = aj for all j {i,i + 1}, then e(a, A) + e(a',\) = for any A. Another 
pair of equations [vLee5, (55, 54)] has an opposite relation to equation [vLee5, (50)], as they contain an 
additional factor of the form s(a, A) in their summand, but they involve neither tableau conditions nor 
Littlewood-Richardson conditions. These different expressions, stated in the form of summations over all 
binary or integral matrices but whose range is effectively restricted by the use of the Iverson symbol, and 
ordered from the largest to the smallest effective range, are as follows. For the binary case they are 




)e{a + vow{M),v) 
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and for the integral case 
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1.3 Crystal operations for binary matrices 

The expressions in (l)-(2)-(3) as well as those in (4)-(5)-(6) are related to one another by successive 
cancellations: in each step one of the factors e(a, A) is replaced by an Iverson symbol that selects only 
terms for which the mentioned factor e(a, A) already had a value 1; this means that all contributions 
from terms for which that factor e(a, A) has been replaced by cancel out against each other. 

The symmetry between the tableau- and Littlewood-Richardson conditions allows us to achieve the 
cancellations form (1) to (3) and from (4) to (6) in an alternative way, handling the second factor of the 
summand first, so that halfway those cancellations one has 

(«a/«|Vm)= E £(« t + col(M),A t )[MGLR [2] W M )] (7) 

MeMW 

in the binary case, and in the integral case 

(sx/ K \a v/ll )= ^ e( K + row(M),A)[MGLR(^/ M )]. (8) 

MEM 

Indeed, the cancellation form (1) to (7) is performed just like the one from (1) to (2) would be for matrices 
rotated a quarter turn (and for \/k in place of v/ p), while the cancellation form (4) to (8) is performed 
just like the one from (4) to (5) would be for the transpose matrices. Slightly more care is needed to 
justify the second cancellation phase, since the Littlewood-Richardson condition in the second factor of 
the summand does not depend merely on row or column sums, as the unchanging first factor did in the 
first phase. In the integral case, the second cancellation phase can be seen to proceed like the cancellation 
from (5) to (6) with matrices transposed, but in the binary case the argument is analogous, but not quite 
symmetrical to the one used to go from (2) to (3). Of course, we already knew independently of this 
argument that the right hand sides of (7) and (8) describe the same values as those of (3) and (6). 

Although for the two factors of the summand of (1) or (4) we can thus apply cancellations to the 
summation in cither order, and when doing so each factor is in both cases replaced by the same Iverson 
symbol, the actual method as indicated in [vLee5] by which terms would be cancelled is not the same 
in both cases. This is so because in the double cancellations leading from (1) to (3) or from (4) to (6), 
whether passing via (2) respectively (5) or via (7) respectively (8), the first phase of cancellation has rather 
different characteristics than the second phase. The first phase is a Gessel-Viennot type cancellation; it 
is general (in that it operates on all terms of the initial summation) and relatively simple (it just needs to 
make sure that a matrix cancels against one with the same row- or column sums). By contrast the second 
phase is a Bender-Knuth type cancellation that only operates on terms that have survived the first phase 
(for matrices satisfying the pertinent tableau condition), and it has to be more careful, in order to assure 
that whenever such a term is cancelled it does so against a term that also survived the first phase. 

The question that motivated the current paper is whether it is possible to find an alternative way 
of defining the cancellations that has the same effect on the summations (so we only want to change 
the manner in which cancelling terms are paired up), but which has the property that the cancellation 
of terms failing one of the (tableau- or Littlewood-Richardson-) conditions proceeds in the same way, 
whether it is applied as the first or as the second phase. This requires the definition of each cancellation 
to be general (in case it is applied first), but also to respect the survival status for the other cancellation 
(in case it is applied second) . The notion of crystal operations on matrices described below will allow us to 
achieve this goal. We shall in fact see that for instance the cancellation that cancels terms not satisfying 
the Littlewood-Richardson condition for v/ fi is defined independently of the shape X/k occurring in the 
tableau condition; in fact it respects the validity of the tableau condition for all skew shapes at once. 

1.3. Crystal operations for binary matrices. 

Consider the cancellation of terms that fail the Littlewood-Richardson condition, cither going from (2) 
to (3), or from (1) to (7). Since the condition M e LR [21 (z^/^) involves partial sums of all columns 
of M to the right of a given one, this condition can be tested using a right-to-left pass over the columns 
of M, adding each column successively to composition that is initialised as fi, and checking whether that 
composition remains a partition. If it does throughout the entire pass, then there is nothing to do, since in 
particular the final value /i + row(M) will be a partition, so that e(n + row(M),v) — [ M € LR [21 (^//i) ] = 
[ (i + row(M) = v ] . If on the other hand the composition fails to be a partition at some point, then one 
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1.3 Crystal operations for binary matrices 

can conclude immediately that M £ LR [21 (y/^i), so the term for M cancels. Up to this point there is no 
difference between a Gessel-Viennot type cancellation and a Bender-Knuth type cancellation. 

Having determined that the term for M cancels, one must find a matrix M' whose term cancels 
against it. The following properties will hold in all cases. Firstly M' will be obtained from M by moving 
entries within individual columns, so that col(M) = col(M'). Secondly, the columns of M that had been 
inspected at the point where cancellation was detected will be unchanged in M', so that the term for M' 
is sure to cancel for the same reason as the one for M . Thirdly, a pair of adjacent rows is selected that 
is responsible for the cancellation; all moves take place between these rows and in columns that had 
not been inspected, with the effect of interchanging the sums of the entries in those columns between 
those two rows. In more detail, suppose (3 is the first composition that failed the test to be a partition, 
formed after including column I (so [3 — fi + J2j>i ^j)> then there is at least one index i for which 
[3 i+ i = fa + 1; one such i is chosen in a systematic way (for instance the minimal one) and all exchanges 
applied in forming M' will be between pairs of entries Mij, Mi+ij with j < I. As a result the partial 
row sums a = J2j<i and a' = wu l be related by (a-,a- +1 ) = (a, + i,aj) (the other parts 

are obviously unchanged), so that fi + row(M) = a + [3 and fi + row(M') = a' + [3 are related in a way 
that ensures e(/i + row(M), v) + e(/i + row(M'), v) = 0, so that the terms of M and M' may cancel out. 

Within this framework, there remains some freedom in constructing M', and here the Gessel-Viennot 
and Bender-Knuth types of cancellation differ. If our current cancellation occurs as the first phase, in 
other words if we are considering the cancellation from (1) to (7), then the fact that we have ensured 
e(ji + row(M), v) + e(fi + row(M'), v) = suffices for the cancellation of the terms of M and M 1 , and 
M' can simply be constructed by interchanging all pairs of bits (M ij 7 M i+ ij) with j < I, which is what 
the Gessel-Viennot type cancellation does (of course such exchanges only make any difference if the 
bits involved are unequal). If however our current cancellation occurs as the second phase (so we are 
considering the cancellation from (2) to (3)), then we must in addition make sure that M £ Tabl [21 (X/k) 
holds if and only if M' £ Tabl [2] (A/f«) does. This will not in general be the case for the exchange 
just described, which is why the Bender-Knuth type of cancellation limits the number of pairs of bits 
interchanged, taking into account the shape X/k for which the tableau condition must be preserved. The 
(easy) details of how this is done do not concern us here, but we note that among the pairs of unequal bits 
whose interchange is avoided, there are as many with their bit '1' in row i as there are with their bit '1' in 
row i + so that the relation between /? and (3' above is unaffected. The alternative construction given 
below similarly ensures that M £ Tabl [2] (A/«;) holds if and only if M' £ Tabl' 21 (A/k) does, but since it is 
defined independently of X/k, it works for all shapes at once, and it can be applied to any matrix, unlike 
the Bender-Knuth cancellation which is defined only for (encodings of) tableaux of shape X/k. 

Our fundamental definition will concern the interchange of a single pair of distinct adjacent bits in 
a binary matrix; this will be vertical adjacent pair in the discussion above, but for the cancellation of 
terms failing the tableau condition we shall also use the interchange of horizontally adjacent bits. Our 
definition gives a condition for allowing such an interchange, which is sufficiently strict that at most 
one interchange at a time can be authorised between a given pair of adjacent rows or columns and in a 
given direction (like moving a bit T' upwards, which of course also involves a bit '0' moving downwards). 
Multiple moves (up to some limit) of a bit in the same direction between the same rows or columns can 
be performed sequentially, because the matrix obtained after an interchange may permit the interchange 
of a pair of bits that was not allowed in the original matrix. 

1.3.1. Definition. Let M £ Ai l2] be a binary matrix. 

a. A vertically adjacent pair of bits (Mk,i, Mk+i.i) is called interchangeable in M if Mk,i ^ Mk+i,i, and 

if Y! 7 Jv M Kj > E l j r Ji> M k+ i,j for all V < I, while ^U+i ^ £5=l+i Mfe+1 ^ for a11 V ~ L 

b. A horizontally adjacent pair of bits (M^^, M/.j+i) is called interchangeable in M if M^ i ^ M^.i+i 

and if J2tk' M i,l < Etfc' M M+i for all k' < k, while EiLk+1 M H > M ^+i for a11 k ' ^ k - 

Applying a upward, downward, leftward, or rightward move to M means interchanging an interchange- 
able pair of bits, which is respectively of the form (Mk,i, Mk+i.i) = (0,1), (Mk,i, Mk+i.i) = (1,0), 
(M M , M M+1 ) = (0, 1), or (M M , M M+1 ) = (I, 0). 

These operations are inspired by crystal (or coplactic) operations, and we shall call them crystal 
operations on binary matrices. Indeed, horizontal moves correspond to coplactic operations (as defined 
in [vLee3, §3]) applied to the concatenation of the increasing words with weights given by the (nonzero) 
rows of M, from top to bottom; vertical moves correspond to coplactic operations on the concatenation 
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of increasing words with weights given by the columns of M, taken from right to left. Applied to the 
binary encoding of a semistandard tableau T, vertical moves correspond to coplactic operations on T. 

This definition has a symmetry with respect to rotation of matrices: if a pair of bits in a finite 
binary matrix is interchangeable, then the corresponding pair of bits in the matrix rotated a quarter turn 
will also be interchangeable. However the definition does not have a similar symmetry with respect to 
transposition of matrices, and this makes it a bit hard to memorise. As a mnemonic we draw the matrices 
(J J) and (J ) with a line between the pairs of bits that are not interchangeable (and they will not be 
interchangeable whenever they occur in a 2 x 2 submatrix of this form, since the conditions allowing 
interchange can only get more strict when a matrix is embedded in a larger one); the pairs not separated 
by a line are in fact interchangeable in the given 2x2 matrices: 



As a somewhat larger example, consider vertical moves in the binary matrix 

/l 00101100000 l\ 
M = 1 1 1 1 1 1 1 1 1 . (10) 
\0 10100111010 1/ 

The pair (*) at the top right is interchangeable, because in every initial part of the remainder of rows 
and 1, the pairs (°) are at least as numerous as the pairs (J). Since they are in fact always strictly 
more numerous, the pair (°) in column 1 is also interchangeable (the closest one comes to violating the 

second inequality in 1.3.1a is the equality J2j= 2 Mo,j = 3 = J2j=2 ^ij, an< ^ tnc nrst inequality poses 
no problems). None of the remaining pairs in rows and 1 are interchangeable however; for the pair 
in column 2 the first inequality in 1.3.1a fails for V = 1 since M^i = ^ M\ y \ — 1, and in fact this 
inequality continues to fail for V = 1 and all further columns (often there are other inequalities that fail 
as well, but one may check that for column 7 the mentioned inequality is the only one that fails). In 
rows 1 and 2, only the pair (J) in column 4 is interchangeable (while all inequalities are also satisfied for 
columns 5 and 12, these columns contain pairs of equal bits (°) and (J), which are never interchangeable). 
As an example of successive moves in the same direction, one may check that, in rows and 1, after 
interchanging the pair (°) in column 1, one may subsequently interchange similar pairs in columns 7, 8, 
10, and 11, in that order. 

Let us now show our claim that at most one move at a time is possible between any given pair of 
rows or columns and in any given direction. Consider the case of adjacent rows, i, i + 1 and suppose they 
contain two interchangeable vertically adjacent pairs of bits in columns jo < ji- Then one has two opposite 
inequalities for the range of intermediate columns, which implies that X^Lj +i = SjLjo+i ^i+i,j- 
One can also see that the interchangeable pair in column jo is (*) and the one in column j\ is (°) , since 
any other values would contradict definition 1.3.1. So there can be at most one downward move and at 
most one upward move that can be applied between rows i and i + 1, with the downward move being 
to the left of the upward move if both occur. Similarly, at most one leftward move and at most one 
rightward move can be applied between a given pair of adjacent columns, and if both can, the leftward 
move is in a row above of the rightward move. 

These uniqueness statements justify the following crucially important definition. In keeping with 
the usual notation for crystal operations, we use the letter e for raising operations and the letter / for 
lowering operations, but since we have a horizontal and a vertical variant of either one, we attach an 
arrow pointing in the direction in which the bit '1' moves. 

1.3.2. Definition, (binary raising and lowering operations) Let M £ Ai [2] . 

a. If M contains an interchangeable pair of bits in rows i and i + 1, then the matrix resulting from the 
interchange of these bits is denoted by e\(M) if the interchange is an upward move, or by //(M) 
if the interchange is a downward move. If for a given i G N the matrix M admits no upward 
or no downward move interchanging any pair bits in rows i and i + 1, then the expression e\(M) 
respectively // (M) is undefined. 

b. If M contains an interchangeable pair of bits in columns j and j + 1, then the matrix resulting from 
the interchange of these bits is denoted by ef (M) if the interchange is a leftward move, or by ff (M) 
if the interchange is a rightward move. If for a given j G N the matrix M admits no leftward or 
no rightward move interchanging any pair bits in columns j and i + 1, then the expression ef(M) 
respectively ff (M) is undchned. 
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Since an interchangeable pair of bits remains so after it has been interchanged, it follows that 
whenever e\(M) is defined then so is //(e}(M)), and it is equal to M. Similarly each of the expressions 
e-(//(M)), ,ej~(fj*(M)) and fp(ey(M)) is defined as soon as its inner application is, in which case 
it designates M. Our next concern will be characterising when expressions such as e\{M) are defined, 
and more generally determining the number of times each of the operations e\ , // , ej~ and /r* can be 
successively applied to a given matrix M, which we shall call the potential of M for these operations. 

1.3.3. Definition. For M G M [2] and i,j G N, the numbers n\{M),n\{M),nj{M),nj{M) G N are 
defined by 

n\{M) = max{^.> ; (M mj - M itj ) \l G N}, (11) 

n\{M) = max{X) j < I (M iiJ - - M i+hj ) | I G N}, (12) 

nr(M) = max{£ i<fc (M iJ+1 - M w ) | fc G N}, (13) 

n;(M) = max { £ l>fe (M,,, - M i;j+1 ) | fc G N }. (14) 



1.3.4. Proposition. For M G -M [2] and i,j G N, the numbers of times each of e\, //, ej~, and can 
be successively applied to M are respectively given by the numbers n\(M), n\{M), n~(M), and nj'(M). 
Moreover n\{M) - n\{M) = row(M)j - row(M) i+ i and nj{M) - nf {M) = col(M)j - col(M)j+i. 

Proof. Suppose M is an n x m binary matrix (so all entries outside that rectangle are zero) that admits 
an upward move interchanging a pair of bits (°) in column I of rows i,i + 1. Then it follows from 

Ejl7ii( M *+ij - Mij) > that n\{M) > Y,j>i( M i+i,i ~ M i,j) > 0- Conversely if nJ(M) > 0, then 
let I < m be the maximal index for which the maximal value of Ej>;(-^+i,j — 1S attained. One 

then verifies that M admits an upward move in column I of rows i, i + 1: the fact that the pair in that 
position is (°) follows from the maximality of I, and failure of one of the inequalities in 1.3.1a would 
respectively give a value I' < I for which a strictly larger sum is obtained, or a value /' + 1 > I for which 
a weakly larger sum is obtained, either of which contradicts the choice of I. 

The statement concerning e\ can now be proved by induction on n\{M). For n\(M) — we have 
just established that no upward moves in rows 1 are possible. So suppose n\(M) > and let M' 

be obtained from M by an upward move in column Iq. Then replacing M by M' decreases the sums 
^2j>l(Mi+i t j — Mij) by 2 for all I < lo, while those sums are unchanged for I > lo- The sums for I < Iq 
therefore become at most n\{M) — 2, while the sums for I > lo remain at most n\(M) ~ 1 (since lo was the 
maximal index for which the value n\ (M) is attained for M, as we have seen). Therefore the maximal sum 
for M' is attained for the index l + 1, and its value is n\{M') = Ej>z +i(-^i+i ,j ~ = n l(M) — 1; by 
induction e\ can be applied precisely that many times to M', and so it can be applied n\{M) times to M 
as claimed. The statements for ej", //, and fj* follow from the statement we just proved by considering 
the (finite) matrices obtained from M by turning it one, two, or three quarter turns. The statements in 
the final sentence of the proposition are clear if one realises that for instance Ej<z(-^i,j " -^i+i.j) an< ^ 
Ej>((-^i+i,j ~ Mij) differ by row(M)i — row(M) i+1 independently of I, so that their maxima n\(M) 
and n\(M) are attained for the same (set of) values of I, and also differ by row(M)i — row(M) i+1 . □ 

With respect to the possibility of successive moves between a pair of adjacent rows or columns, we 
can make a distinction between pairs whose interchange is forbidden in M but can be made possible after 
some other exchanges between those rows or columns, and pairs whose interchange will remain forbidden 
regardless of such exchanges. We have seen that when a move is followed by a move in the opposite 
direction, the latter undoes the effect of the former; it follows that if a given move can be made possible 
by first performing one or more moves between the same pair of rows or columns, then one may assume 
that all those moves are in the same direction. Moreover we have seen for instance that successive upward 
moves between two rows always occur from left to right; this implies that if a pair (®) in column I of 
rows i,i + 1 is not interchangeable due to a failure of some instance of the second inequality in 1.3.1a 
(which only involves columns j > I), then this circumstance will not be altered by any preceding upward 
moves between the same rows, and the move will therefore remain forbidden. On the other hand if the 
second inequality in 1.3.1a is satisfied for all /' > I, then the value of Ej>;(^+iJ ~ Mij) is larger 
than the one obtained by replacing the bound I by any I' > I; it may still be less than that overall 
maximum n\ (M), but that value can be lowered by successive upward moves between rows + which 
must necessarily occur in columns j < I, until the pair (J considered becomes interchangeable. 
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We may therefore conclude that, in the sense of repeated moves between two adjacent rows, failure 
of an instance of the first inequality in 1.3.1a gives a temporary obstruction for a candidate upward 
move, while failure of an instance of the second inequality gives a permanent obstruction. For candidate 
downward moves the situation is reversed. The following alternative description may be more intuitive. 
If one represents each pair (J by "(", each pair (J) by ")", and all remaining pairs by "— " (or any 
non-parenthesis symbol), then for all parentheses that match another one in the usual sense, the pairs 
in the corresponding columns are permanently blocked. The remaining unmatched parentheses have the 
structure ") •••)(••• (" of a sequence of right parentheses followed by a sequence of left parentheses (either 
of which might be an empty sequence). An upward move between these rows is possible in the column 
corresponding to the leftmost unmatched "(" if it exists, an a downward move between these rows is 
possible in the column corresponding to the rightmost unmatched ")" if it exists. In cither case the move 
replaces the parenthesis by an opposite one, and since it remains unmatched, we can continue with the 
same description for considering subsequent moves. In this view it is clear that all unmatched parentheses 
can be ultimately inverted, and that upward moves are forced to occur from left to right, and downward 
moves from right to left. For instance, in the 3 x 13 matrix given as an example after definition 1.3.1, 
the sequence of symbols for the two topmost rows is ") ( (— ( ))((—( (— " , and from this it is clear that 
one downward move is possible in column 0, or at most 5 successive upward moves in columns 1, 7, 8, 

10, and 11; for the bottommost rows we have the sequence " )— ) — ( )— " and only successive 

downward moves are possible, in columns 4 and 2. For moves between adjacent columns the whole picture 
described here must be rotated a quarter turn (clockwise or counterclockwise, this makes no difference). 

We now consider the relation of the definitions above to the tableau- and Littlewood-Richardson 
conditions on matrices. The first observation is that these conditions can be stated in terms of the 
potentials for raising (or for lowering) operations. 

1.3.5. Proposition. Let M G M l2] and let X/k and u./v be skew shapes. 

(1) M G Tabl [21 (A/«;) if and only if col(M) = A* - k* and n*f{M) < k) - for all j G N. 

(2) M G \M [2] {v/u) if and only if row(M) = v - u. and n\{M) < ft - ft+i for all i G N. 

The second parts of these conditions can also be stated in terms of the potentials of M for lowering 
operations, as nj(M) < X* — Aj +1 for all j G N, respectively as n^(M) < Vi — z/j+i for all ieN. 

Proof. In view of the expressions in definition 1.3.3, these statements are just reformulations of the 
parts of proposition 1.1.1 and definition 1.1.2 that apply to binary matrices. □ 

The next proposition shows that vertical and horizontal crystal operations on matrices respect the 
tableau conditions respectively the Littlewood-Richardson conditions for all skew shapes at once. 

1.3.6. Proposition. If binary matrices M, M' G Ai l2] are related by M' = e\{M) for some i G N, then 
ny(M) = ny(M') and nj{M) = nj{M') for all j G N. Consequently, the conditions M G Tabl [21 (A/ft) 
and M 1 G Tabl [21 (A/re) arc equivalent for any skew shape X/k. Similarly if M and M' are related 
by M' = ff{M) for some j G N, then n\{M) = n\{M') and n\{M) = n\{M') for all i G N, and 
M G LR [21 (i7^) 4=^ M' G LR [21 (iV^i) for any skew shape u/fi. 

Proof. It suffices to prove the statements about M' = e\(M), since those concerning M 1 — eJ~(M) 
will then follow by applying the former to matrices obtained by rotating M and M' a quarter turn 
counterclockwise. For the case considered it will moreover suffice to prove nj(M) = ny(M') for 
any j G N, since nj*(M) = nj*(M') will then follow from col(M) = col(M'), and the equivalence of 
M G Tabl [21 (A/K) and M' G Tabl [2) (A/K) will be a consequence of proposition 1.3.5. One may suppose 
that the pair of bits being interchanged to obtain M' from M is in column j or j + 1, since otherwise 
W~(M) — ny(M') is obvious from (13). Let (pk)ke~N be the sequence of partial sums for M of which 
n*"(M) is the maximum, in other words pt — J2i'<k(Mi> j+i— Mi',j), and let (p' k )k^ be the corresponding 
sequence for M'. Then the only index k for which pk ^ p' k is k = i + 1: one has p' i+1 = Pi + \ — 1 if 
the move occurred in column j, or p' i+1 = Pi+\ + 1 if it occurred in column j + 1. The only way in 
which this change could make n~(M) = m&XkPk differ from n^(M') = max^pj. is if k = i + 1 were 
the unique index for which pk = n^(M) (in the former case) or for which p' k = n^(M') (in the latter 
case) . That would in particular require that the indicated value be strictly larger than pu = p' k and than 
Pfe+2 = p'k+2-> so M or M' would have to contain a submatrix (° at the intersection of rows i, i + 1 and 
columns j, j + 1, while the other matrix would differ by the interchange of one of those two vertically 
adjacent pairs of bits. But we have seen that in such a submatrix neither of those two pairs of bits can be 
interchangeable, which excludes this possibility, and one therefore has n^{M) = ny(M') in all cases. □ 
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One can summarise the last two propositions as follows: Littlcwood-Richardson conditions can be 
stated in terms of the potentials for vertical moves, which moves preserve tableau conditions, while 
tableau conditions can be stated in terms of the potentials for horizontal moves, which moves preserve 
Littlcwood-Richardson conditions . 

We shall now outline the way in which crystal operations can be used to define cancellations either of 
terms for matrices not in Tabl [21 (A/«;) or of those not in LR [21 (y//i), in the summations of (1), (2), or (7). 
One starts by traversing each matrix M as before, searching for a violation of the condition in question, 
and of an index that witnesses it; this amounts to finding a raising operation e (i.e., some ej~ or e\) for 
which the potential of M is larger than allowed by the second part of proposition 1.3.5 (1) or (2). 

Now consider the set of matrices obtainable from M by a sequence of applications of e or of its inverse 
lowering operation /; these form a finite "ladder" in which the operation e moves up, and / moves down. 
Note that the potential for e increases as one descends the ladder. The condition of having a too large a 
potential for e determines a lower portion containing M of the ladder, for which all corresponding terms 
must be cancelled, and the witness chosen for such a cancellation will be the same one as chosen for M 
(there may also be terms cancelled in the remaining upper portion of the ladder, but their witnesses will 
be different). Now the negation of the expression of the form e(a, A), which is required for cancellation, 
can be obtained by reversing the indicated lower part of the ladder. Since a pair of matrices whose terms 
cancel are thus linked by a sequence of horizontal or vertical moves, their status for any Littlewood- 
Richardson respectively tableau condition (the kind for which one is not cancelling) will be the same, 
which allows this cancellation to be used as a second phase (starting from (7) or (2)). 

Let us fill in the details of the description above, for the cancellation of terms for matrices not in 
LR [2] [y I fjb) , in other words leading from (1) to (7) of from (2) to (3). As described at the beginning of 
this subsection, we start by finding the maximal index I such that the composition (3 = \i + J2j>i * s 
not a partition, and choosing an index i for which /3 i+ i = fa + 1; this implies that n\ (M) > fn — /Uj+i, so 
the potential of M for e = e\ exceeds the limit given in 1.3.5 (2). For convenience let us use the notation 
e d to stand for (e\) d when d > 0, for the identity when d = 0, and for (f-)~ d when d < 0; then the ladder 
mentioned above is { e d (M) \ —n\(M) < d < n\(M) }, and its lower part that must be cancelled because 
of a too large potential for e\ is { e d (M) \ -n\(M) <d< n\(M) - (^ - }. 

From the maximality of / it follows that M contains a pair (°) in column I of rows i,i + 1, and 
that this pair is not permanently blocked for upward moves in rows i,i + 1 (in other words, one has 
jyjLi+i Mij < Y^jLi+i Mi+ij for all m > I); indeed the pair will be interchanged in e d (M) when d 
satisfies n\(M) — (in — /i,+i) < d < n\(M), i.e., in the mentioned upper part of the ladder. So the 
lower part of the ladder is precisely the part in which that pair is not interchanged, and the matrices in 
this part will give rise to the same indices I and i as M to witness their cancellation. The expression 
for d such that the term for M cancels against the one for e d (M) can be found as follows. If M is at 
the bottom of the ladder (n\(M) = 0) then d has the value n\(M) — (m — m+i) — 1 that gives the 
topmost value of the bottom part of the ladder, and d decreases with the level nj(M) of M, so the 
expression is d — n\(M) — (fii — (M+i) — 1 — n\(M). Putting a = \i + row(M), this can also be written 
as d — cti+i — on — 1, by proposition 1.3.4. Since each application of e\ increases the sum of entries in 
row i while decreasing the sum in row i + 1, the value a' = \i + row(M') for the matrix M' = e d (M) 
satisfies (a-,a- +1 ) = (c*i+i — l,Q!i + 1) while its remaining components are unchanged from a, which 
ensures that e(a, v) + e(a', v) = 0. The fact that M and M' are related by vertical moves implies that 
M £ Tabl [21 (A/ft) <^=> M' E Tabl [21 (A/«;) for any skew shape X/k, so the terms for M and M' do indeed 
cancel, whether we are considering the passage from (1) to (7) or the one from (2) to (3). 

For the cancellations involved in passing from (1) to (2) and from (7) to (3) the description is similar, 
but rotated a quarter turn counterclockwise: the initial scan of the matrix is by rows from top to bottom, 
and the raising operations n\ are replaced by raising operations nj . 

The reader may have been wondering whether we have been going through all these details just to 
obtain more aesthetically pleasing descriptions of the reductions (1)— >(2)— >(3) and (1)— >(7)— >(3) (and 
maybe the reader even doubts whether that goal was actually obtained). But crystal operations turn out 
to be useful in other ways than just to define cancellations, and several such applications will be given 
below; those applications alone largely justify the definition of crystal operations. We have nevertheless 
chosen to introduce them by considering cancellations, because that provides a motivation for the precise 
form of their definition and for treating separate cases for binary and integral matrices; such motivation 
might otherwise not be evident. For our further applications it is of crucial importance that horizontal 
and vertical moves are compatible in a stronger sense than expressed in proposition 1.3.6. Not only 
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do moves in one direction leave invariant the potentials for all moves in perpendicular directions, they 
actually commute with those moves, as stated in the following lemma. 

1.3.7. Lemma, (binary commutation lemma) Let M,M',M" G M [2] be related by M' = e\(M) and 
M" = e^(M) for some i,j £ N; then e'~(M') = e\(M"). The same holds when e\ is replaced both times 
by ff and/or ej~ is replaced both times by fj*. 

Proof. Note that the expressions eJ~(M') and e\(M") are defined since nj(M') = n^(M) > and 
n\(M") = n\(M) > by proposition 1.3.6. The variants given in the second part of the lemma can be 
deduced from the initial statement either by rotation symmetry or by a suitable change of roles between 
the four matrices involved. So we shall focus on proving the initial statement. 

Suppose first that the pairs of bits interchanged in the moves c\:M i— > M' and Cj:M i— > M" are 
disjoint. In this case we shall argue that these pairs of bits are in the same position as the pairs of bits 
interchanged in the moves c\: M" *—> e](M") and ej: M' ^> e^(M'), respectively; then it will be obvious 
that e^(e\(M)) = e](eJ~(M)). To this end must show that the conditions in definition 1.3.1, which are 
satisfied in M for each of the two pairs of bits considered, remain valid after the other pair is interchanged. 
Since the values of one pair of bits is not affected by the interchange of the other pair, we only need to 
worry about the four inequalities in that definition. Depending on the relative positions of the two pairs, 
at most one of those inequalities can have an instance for which the values being compared change, but 
since we do not know which one, this does not help us much; nevertheless the four cases are quite similar, 
so we shall treat only the first one explictly. Each inequality, with its quantification, can be reformulated 
as stating that some maximum of partial sums does not exceed (actually it equals 0); for instance the 
first inequality is equivalent to 'max { X^=l' (-^fc+i ,j ~ ^k,j) < I' < I } < 0' (this condition applies for 
k = i if the move of e\ occurs in column I). That maximum of partial sums is of the same type as the 
one in one of the equations (11)— (14), but for a truncated matrix; in the cited case they are the partial 
sums of (11) but computed for M truncated to its columns j < I. Therefore the same reasoning as in 
the proof of proposition 1.3.6 shows that although one of the partial sums may change, their maximum 
remains the same, so that the pair of bits considered remains interchangeable. 

Now suppose that to the contrary the pairs of bits (°) and (0 1) being interchanged in M do overlap. 
Then after performing one interchange, the pair of bits in the position of the other pair can no longer be 
interchangeable, as its bits will have become equal. There is a unique 2x2 submatrix of M that contains 
the two overlapping pairs, and since it contains both a vertical and a horizontal interchangeable pair of 
bits, its value can be neither (° J) nor (J °) . Therefore it will contain either (° J) if the two pairs overlap 
in their bit '0' (at the top left), or (° °) if the two pairs overlap in their bit '1' (at the bottom right). 
In either case it is not hard to see that the overlapping bit, after having been interchanged horizontally 
or vertically, is again (in its new position) part of an interchangeable pair within the 2x2 submatrix, 
in the direction perpendicular to the move performed; the other bit of that pair is the one in the corner 
diametrically opposite to the old position of the overlapping bit in the submatrix considered (the bottom 
right corner in the former case and the top left corner in the latter case). This is so because comparing 
that new pair with the interchangeable pair that used to be in the remaining two squares of the 2x2 
submatrix, the only difference for each of the pertinent inequalities of definition 1.3.1 is the insertion or 
removal of a bit with the same value in each of the two sums being compared, which does not affect 
the result of the comparison. Therefore the succession of two raising operations, applied in either order, 
will transform the submatrix (" J) into (J J), or the submatrix (° °) into (J [[), as illustrated below. □ 
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1.4. Crystal operations for integral matrices. 

Motivated by the existence of cancellations (4)— »(5)— >(6) and (4)— »(8)— >(6), we shall now define opera- 
tions like those defined in the previous subsection, but for integral instead of binary matrices. Much of 
what will be done in this subsection is similar to what was done before, so we shall focus mainly on the 
differences with the situation for binary matrices. 
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1.4 Crystal operations for integral matrices 

A first difference is the fact that for integral matrices the operation of interchanging adjacent entries 
is too restrictive to achieve the desired kind of modifications. We shall therefore regard each matrix 
entry m as if it were a pile of m units, and the basic type of operation will consist of moving some units 
from some entry m > to a neighbouring entry, which amounts to decreasing the entry m and increasing 
the neighbouring entry by the same amount. We shall call this a transfer between the two entries; as in 
the binary case we shall impose conditions for such a transfer to be allowed. Another difference is the kind 
of symmetry implicitly present in equation (6) compared with (3), which in fact stems from the difference 
between the cases of binary and integral matrices in the relation of definition 1.1.2 to proposition 1.1.1, 
which we already observed following that definition. As a consequence the rules for allowing transfers 
will not be symmetric with respect to rotation by a quarter turn, but instead they will be symmetric 
with respect to transposition of the integral matrices and with respect to rotation by a half turn. 

This new type of symmetry somewhat simplifies the situation, but there is also a complicating factor, 
due to the fact that the tableau- and Littlewood-Richardson conditions are more involved for integral 
matrices than for binary ones. In the binary case it sufficed to construct a sequence of compositions by 
cumulating rows or columns, and to test each one individually for being a partition. But in the integral 
case one must test for each pair a, [3 of successive terms in the sequence whether a t— fa in other words 
whether /3/a is a horizontal strip. That test amounts to verifying fa + \ < on for all i 7 since cti < fa 
already follows from the circumstance that fa is obtained by adding a matrix entry to on. Thus if we 
focus on the inequalities involving the parts i and i + 1 of the compositions in the sequence, then instead 
of just checking that part i + 1 never exceeds part i of the same composition, one must test the stronger 
requirement that part i + 1 of the next partition in the sequence still does not exceed that (old) part i. 

This will mean for the analogues of definitions 1.3.1 and 1.3.3, that the final entries in partial sums 
in two adjacent rows or columns will not be in the same column or row, but in a diagonal position 
with respect to each another (always in the direction of the main diagonal). This also means that the 
conditions required to allow a transfer must take into account some of the units that are present in the 
matrix entries between which the transfer takes place, but which are not being transferred themselves (in 
the binary case no such units exist). Although the precise form of the following definition could be more 
or less deduced from the properties we seek, we shall just state it, and observe afterwards that it works. 

1.4.1. Definition. Let M e M, k, I G N, and a e Z - {0}. 

a. Suppose that J2j^i'(Mk+i,j+i — ^k,j) > max(a,0) for all V < I, or if I = that M k+1 ^ > a, and 

suppose Y^jZi (Mk^j— Mk+i.j+i) > max(— a,0) foralll' > I. Then we allow the entries (Mk,i, Mk+ij) 
to be replaced by (Mk,i + a, Mk+i.i — a); this is called an upward transfer of a units between rows 
k and k + 1 if a > 0, or a downward transfer of —a units between those rows if a < 0. 

b. Suppose that J2iZk< (Mi+i.i+i — Mij) > max(a,0) for all k' < k, or if k = that Mo^+i > a, and 
suppose ^iSf, 1 (Mij — M i+ u + i) > max(— a,0) forallk' > k. Then we allow the entries (Mk,i, Mkj+i) 
to be replaced by (Mk,i + a, Mkj+i — a); this is called a leftward transfer of a units between columns 
I and I + 1 if a > 0, or a rightward transfer of —a units between those columns if a < 0. 

Remarks. (1) The occurrence of the quantity a in the inequalities has the effect of cancelling its contri- 
bution to the entry from which it would be transferred. It follows that the transfer can always be followed 
by a transfer of a units in the opposite sense between the same entries, which reconstructs the original 
matrix. (2) The exceptional conditions Mo^+i > a and Mk+\,o > a compensate for the absence of any 
inequality where a occurs in the way just mentioned. They serve to exclude the introduction of negative 
entries by a transfer; note that for instance this is taken care of for upward moves with I > by the 
condition M^+i^i — Mk t i-i > a, and for downward moves by M^/ — Mk+i.i+i > —o- Hence the cases k = 
and I = arc not treated any differently from the others. (3) Wc could have restricted the definition 
to the cases a = 1 and a = — 1, since transfers with |a| > 1 can be realised by repeated transfers with 
\a\ = 1. The current formulation was chosen for readability, and because it permits a straightforward 
and interesting generalisation to matrices with non-negative real coefficients. 

We shall call these transfers crystal operations on integral matrices. It can be seen that horizontal 
transfers correspond to |o| successive coplactic operations on the word formed by concatenating weakly 
decreasing words whose weights are given by the (nonzero) rows of M, taken from top to bottom; vertical 
transfers correspond \a\ successive to coplactic operations on the word similarly formed by concatenating 
weakly decreasing words with weights given by the columns of M, taken from left to right. Horizontal 
transfers in the integral encoding of a semistandard tableau T correspond to coplactic operations on T. 
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1.4 Crystal operations for integral matrices 

Here is a small example of vertical transfers; for an example of horizontal transfers one can transpose 
the matrix. Consider vertical moves between the two nonzero rows of the integral matrix 

M= ( 2 lU 2 5 2 0j' (16) 

For the three leftmost columns no transfer is possible because Y^j=ii^o,j — M\j + i) < for I G {0, 1, 2}. 
For column 3 however an upward transfer of at most 2 units is possible, where the limit comes from the 
value J2j = i(Mij + i — M j) = 2 (and starting the sum at j = would give the same value). No downward 
transfer is possible in that column because Ej=3(-^o,j — M\j + i) = 0; however a downward move of at 
most 4 units is possible in column 7. If all 4 units are transferred downwards in that column it is replaced 
by (j?) , and no further downwards between these rows will be possible; however if 2 units are transferred 
upwards in column 3 one obtains the matrix 

M'=( l 2 1 5 3 1 2 4 0\ 

^2 1 1 2 2 5 2 y ' { ' 

and an upwards transfer in M' is still possible, namely in column where both units from row 1 can be 
transferred to row 0. 

We now consider the ordering of transfers between a given pair of adjacent rows or columns. Since 
any simultaneous transfer of more than one unit can be broken up into smaller transfers, we cannot claim 
that at most one transfer in a given direction between a given pair of rows or columns is possible, but 
it is important that the only choice concerns the amount transferred, not the pair of entries where the 
transfer takes places. Now if vertical transfers of a and a' units are possible between rows i and i + 
respectively in columns jo and ji with j < ji, then the inequalities in 1.4.1a give the inequalities 
min(o, 0) > X^=j (-^j+M+i — — max(a',0) of which the middle member must be 0, and which 

together with a, a' ^ therefore imply a > > a'. This shows that the transfer in column j is upwards 
and the one in column j\ is downwards. Thus successive upward transfers between the same pair of rows 
take place in columns whose index decreases weakly, and successive downward transfers between these 
rows take place in in columns whose index increases weakly; the same is true with the words "rows" 
and "columns" interchanged. Note that these directions are opposite to the ones in the binary case 
for transfers between rows, but they are the same as the ones in the binary case for transfers between 
columns. The following definition is now justified. 

1.4.2. Definition, (integral raising and lowering operations) Let M g M and i,j E N. 

a. If M admits an upward or downward transfer of a single unit between rows i and i+1, then the 
resulting matrix is denoted by e\ (M) respectively by // (M) . If M admits no upward or no downward 
transfer between rows i and i + 1, then the expression e\(M) respectively //(M) is undefined. 

b. If M admits a leftward or rightward transfer of a single unit between columns j and j + 1, then 
the resulting matrix is denoted by eJ~(M) respectively by fp(M). If M admits no leftward or 
no rightward transfer between columns j and j + 1, then the expression e'~(M) respectively fp(M) 
is undefined. 

If an upward, downward, leftward, or rightward transfer of a > 1 units between a given pair of 
rows i,i + 1 or columns j, j + 1 is defined, the the resulting matrix can also be obtained by a succession 
of a raising or lowering operations, as (e-) a (M), (fl) a (M), (e^) a (M) or (jj*) a (M), respectively. A 
succession of even more transfers between the same rows or columns and in the same direction may be 
possible, and the potentials for such transfers are given by the following expressions. 

1.4.3. Definition. For M e M and i,j e N, the numbers nl(M),n l i (M),ny(M),nj'(M) 6 N are 
dehned by 



n\(M) = max {M i+lfi + Y,j<i{Mi +w - M itj ) | I e N}, (18) 

n\{M) = max{X) J -> I (M i j - M i+W ) ieN}, (19) 

n-{M) = max{M 0j+ i + E l<fc ( M t +U+i - M i,j) ^N}, (20) 

nJ{M) = max{E i > fe (M ilj - M i+1J+1 ) \ k e N}. (21) 
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1.4 Crystal operations for integral matrices 



1.4.4. Proposition. For M E M and i,j G N, the numbers of times each of e\, //, , ej~ and fj° can 

be successively applied to M are given respectively by the numbers n\(M), n\{M), n*~(M), and nj*(M). 
Moreover n\{M) - n\(M) = row(M), - row(M) l+1 and nj{M) - nj{M) = co\{M)j - co\(M) j+1 . 

Proof. If M admits an upward transfer of a units between rows i and i + 1 and in column I, then 
n\{M) > M l+lfl + Ej<i( M »+i,j'+i - M i,j) > M i+i,o + a > if I > 0, while one has n\{M) > M i+lfi > 
a > in case I = 0, so nJ(Af) is nonzero either way. Similarly if M admits a downward transfer of 
a units between those rows then n\{M) > a > 0. Now suppose conversely that n\(M) > or that 
n\{M) > 0. In the former case let lo be the minimal index I for which the maximum in (18) is attained, 
and in the latter case let l\ be the maximal index I for which the maximum in (19) is attained; in 
either case let the maximum exceed all values at smaller respectively at larger indices by a. Then it is 
easily verified that an upward transfer of a units in column l , respectively a downward transfer of a 
units in column l\ 1 is possible between rows + We know that this is the only column in which a 
transfer in that direction between rows i,i + 1 is possible. Moreover the expressions in (18) and (19) 
have a constant difference row(M) i+ i — row(M)i for every /, so we may conclude that, whenever a 
transfer in either direction between M^i and Mi + u is possible, the index I realises the maxima defining 
n\{M) and n\(M) in both expressions. Then the fact that the transfer can be followed by an inverse 
transfer shows that an upward transfer of a units decreases n\(M) by a, and that a downward transfer 
of a units decreases n\(M) by a; a straightforward induction on n\{M) or n\(M) will complete the 
proof. The statements concerning nJ~(M) and nj{M) follow by transposition of the matrices involved. □ 

Like for binary case, we can describe the possibilities for transfers between a given pair of rows or 
columns using parentheses. We shall describe the case of rows i,i + 1, where the notion of matching 
parentheses is most easily visualised; for the case of adjacent columns one must apply transposition of 
the matrix (unlike the binary case, where rotation must be used) . Since successive transfers between two 
rows progress in the sense opposite to the binary case, namely from right to left for successive upward 
moves, we must invert the correspondence between parentheses and moves: an upward transfer of a unit 
will correspond to the transformation of ")" into "(". Each unit in row i is a potential candidate for an 
downward transfer, and therefore represented by "(", and each unit in row i + 1 is represented by ")". 
All these symbols are gathered basically from left to right to form a string of parentheses, but the crucial 
point is how to order the symbols coming from a same column j. The form of the summations in the 
definitions above makes clear that the rule must be to place the Mij symbols "("to the right of the 
Mj+ij symbols ")", so that these cannot match each other; rather the symbols from Mij may match 
those from Mi+ij+i. Now, as for the binary case, the units that may be transferred correspond to 
the unmatched parentheses, and the order in which they are transferred is such that the corresponding 
parentheses remain unmatched: upward moves transform unmatched symbols ")" into "(" from right to 
left, and downward moves transform unmatched symbols "(" into ")" from left to right. In the example 
given the string of parentheses is ))(|)((|)(|))))(((|))(((| (|)))))((|))(((G where we have inserted bars to 
separate the contributions from different columns, and underlined the maximal substrings with balanced 
parentheses; this makes clear that 4 successive upward transfers are possible in columns 3,3,0,0, or 
4 successive downward transfers, all in (the final nonzero) column 7. 

Note that any common number of units present in entries Mij and Mj+ij+i will correspond to 
matching parentheses both for transfers between rows i,i + 1 and between columns j, j + 1. Therefore 
no transfer between those rows or those columns will alter the value min(M ij j, Mj+ij+i) (but it can be 
altered by transfers in other pairs of rows or columns). In fact one may check that those instances of 
the inequalities in definition 1.4.1 whose summation is reduced to a single term forbid any such transfer 
involving Mjj or M i+lj+1 when that entry is strictly less than the other, and in case it is initially greater 
or equal than the other they forbid transfers that would make it become strictly less. 

The relation of the above definitions to tableau- and Littlewood-Richardson conditions for integral 
matrices is given by the following proposition, which like the one for the binary case is a direct translation 
of the pertinent parts of proposition 1.1.1 and definition 1.1.2. 

1.4.5. Proposition. Let M G Ai and let X/k and \ijv be skew shapes. 

(1) Me Tabl(A//t) if and only if row(M) = A - k and n\ (M) < - n t+1 for all i G N. 

(2) M G LR(jVm) ^ and on ly ^ col(M) = v - fj, and nf(M) < fij - fi j+1 for all j G N. 

The second parts of these conditions can also be stated in terms of the potentials of M for lowering 
operations, as n\{M) < \ — X i+1 for all i G N, respectively as nj'(M) < i/j — Vj+i for all j G N. □ 
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1.4 Crystal operations for integral matrices 

And as is the binary case, the potentials for transfers in a perpendicular direction are unchanged. 

1.4.6. Proposition. If integral matrices M, M' E M. are related by M' = e\(M) for some i g N, then 
n t-(M) = rij-(M') and ny(M) = nj{M') for all j E N. Consequentfy, the conditions M E LR(zVm) 
and M' E LR(^/^) are equivalent for any skew shape vj\i. Similarly if M and M' are related by 
M' = ff(M) for some j E N, then n\(M) = n\(M') and n\(M) = n\(M') for all i E N, and one has 
M E Tabl(A/re) ^> M' E Tabl(A/«;) for any skew shape X/k. 

Proof. Like for the binary case we may focus on proving nJ~(M) = ny(M') when M 1 — e\(M), 
and we may assume that the upward transfer involved in passing from M to M' occurs in column 
j or j + 1. Suppose first that it occurs in column j. Then the only partial sum in (20) that differs 
between M and M' is the one for k = i + 1, which has Mi + ij + \ — Mij as final term; it will decrease 
by 1 in M' . But this decrease will not affect the maximum taken in that equation unless it was strictly 
greater than the previous partial sum, for k = i, which means that M^j < M i+lj+1 ; however that 
would contradict the supposition that e\ involves Mij, so it does not happen. Suppose next that the 
upward transfer occurs in column j + 1. Then the only one of the values of which the maximum is 
taken in (20) that differs between M and M' is the one for k = i: it is Mij + \ if i = 0, and otherwise 
contains a partial sum with final term Mij+i — Mi-ij; it will increase by 1 in M' in either case. 
But that increase will not affect the maximum unless the value for k = i is at least as great as the 
one for k = i + 1, which means Mij > Mj+ij+i, but this would contradict the supposition that e\ 
decreases M i+ ij + i, so this does not happen either. Therefore n^(M) = n^(M') holds in all cases. □ 

Using crystal operations for integral matrices, it is possible to describe cancellations that will realise 
the transitions (4)^(5)— >(6) and (4)— >(8)— >(6), as we have done for binary matrices. We shall describe 
the cancellation of the terms not satisfying M E LR(v/[i), which realises (5)^(6) and (4)^(8); the 
cancellation of the terms not satisfying M E Tabl(«/A) is similar, transposing all operations. One 
starts traversing M by rows from top to bottom, searching for the first index k (if any) for which one 
has (/i + X^j<fe Mi) ^f- (/U + X)i<fe Mi), and then choosing the maximal witness j for this, i.e., with 
{u, + J2i<k < (m + J2i<k ■ This implies that nj~ (M) > u,j — , and similarly to the binary 

case the cancellation is along the "ladder" of matrices obtainable from M by operations e*j and fj*, 
reversing its lower portion where the potential for exceeds — /ij+i: the term for M is cancelled 
against the term for M' = (e^) d (M), where d = nj(M) - (p, 3 - + 1) - nj{M), and as before {ej) d 
means {ff)~ d if d < 0. 

One can also write d = ctj + i — otj — 1 where a = /i + col(M), from which it easily follows that 
the terms for M and M' do indeed cancel each other. The fact that the same values k and j will be 
found for M' (so that we have actually defined an involution on the set of cancelling terms) is most 
easily seen as follows, using the parenthesis description given above. Throughout the lower portion of 
the ladder, the unmatched parentheses ")" from the left, up to and including the one corresponding to 
the unit of Mkj+i that "causes" the inequality (fi + J2i<k < (fi + X^i<fe -^Oj+i) are unchanged, so 
M[ = Mi for all i < k, and the entry M' k - +1 remains large enough that the mentioned inequality still 
holds for M' . Since no coefficients have changed that could cause any index larger than j to become a 
witness for (n + ^2 i<k M-) -j- {p- + ^Z i<k the index j will still be the maximal witness of that relation 
for M'. (The change to the entry M k j could cause the index j — 1 to become, or cease to be, another 
witness of this relation for M', so it is important here that j be chosen as a maximal witness, unlike in 
the binary case where any systematic choice works equally well.) 

Again it will be important in applications to strengthen proposition 1.4.6, in a way that is entirely 
analogous to the binary case. The proof in the integral case will be slightly more complicated however. 

1.4.7. Lemma, (integral commutation lemma) Let M, M' , M" E M, be related by M' — e\(M) and 
M" = e^(M) for some i,jE N; then e*"(M') = e\(M"). The same holds when e\ is replaced both times 
by // and/or e*~ is replaced both times by fj*. 

Proof. As in the binary case, it suffices to prove the initial statement, and both members of the equation 
eJ~(M') = e\(M") are well defined by proposition 1.4.6. That equation obviously holds in those cases 
where the applications of ej~ and e\ occurring in it are realised by unit transfers between the same pairs 
of entries as in the application of these operations to M. This will first of all be the case when the pairs 
of entries involved in the transfers e\:M M' and ep M *—> M" are disjoint: in that case, for each 
inequality in definition 1.4.1 required for allowing one transfer, an argument similar to the one given in 
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the proof of proposition 1.4.6 shows that the minimum over all values I' or k' of its left hand member is 
unaffected by the other transfer. 

In the remaining case where the pairs of entries involved in the two transfers overlap, they lie inside 
the 2x2 square at the intersection of rows i,i + 1 and columns j,j + 1. Since we are considering 
upward and leftward transfers, there are only two possibilities: if Mij > M i+lj+1 , then both transfers 
are towards Mij, while if Mij < M i+ ij + i then both transfers are from M i+ i J+ i. In either case there is 
for each transfer only one of the inequalities in definition 1.4.1 whose left hand member is affected by the 
other transfer, namely the one involving just Mij and Mi + \j + \. In case both transfers are towards Mij, 
the inequalities in question both read Mij — Mj+i^+i > max(— 1,0) = 0, and they continue to hold 
when Mij is increased in M 1 and in M" . In case both transfers are from Mj+i^+i, the inequalities both 
read Mj+ij+i — Mij > max(l, 0) = 1, and here the other transfer may invalidate the inequality: indeed 
the inequality will continue to hold, when Mj+ij+i is decreased by 1 in M' and in M", if and only if 
Mi+ij+i > Mij + 2. Therefore the only case where the unit transfers involved in e^(M') and in e\(M") 



are not in the same squares as in e*"(M) and in e-(M), respectively, is when Mj+ij+i = Mjj + 1. In 
that case an elementary verification, comparing the required inequalities in M' and M" to those that are 
known to hold in M, shows that those transfers are from M t ' - +1 to M[ ^ respectively from M" +1 j to M"j; 
then ej{M') — e\(M") still holds, as illustrated by the display below of the relevant 2x2 squares. 
Instead of that verification, one may also argue that the transfer involved in e^(M') or in e\{M"), which 
as we saw cannot occur in the same place as the corresponding transfer in M, cannot occur in a pair of 
squares disjoint from those of the perpendicular transfer preceding it either (such a transfer in disjoint 
squares cannot validate a relevant inequality any more than it could invalidate one in the first case 
considered above); a transfer towards the entry at (i, j) is therefore the only possibility that remains. □ 

+ 1 y \ ( e T fx y + 1 
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By a simple induction, the commutation lemma implies that {ej) m ((ej) n (M)) = (e\) n ((ej) m (M)) 
whenever (eJ) n (M) and (ej) m (M) are defined (this is true in the binary case as well). It is interesting 
to consider how this commutation is realised if one assumes that the operations applied to M can each 
be performed by a single transfer. Much of the proof above can be applied to such transfers as well, but 
the final case is more complicated. One can see that transfers remain in the same position and are by 
the same amount, unless those in (e\) n (M) and (e~~) m (M) are both from the entry Mk+ij+i and one 
has Mk y i < Mfc+i^+i < Mk.i + n + m; then in fact Mkj + max(n, m) < Mk+i.i+i by the assumption that 
(eJ)"(M) and (e*~) m (M) are defined by single transfers. In this case one can put d = Mk+i.i+i — Mk.i 
and decompose {ej) m — {ej) a o (ej~) b and (e-) n = (e-) a o {e\) c where a = m + n — d and b = m — a, 
c = n — a; then the commutation takes place as in the diagram below, which shows the evolution of the 
2x2 square at the intersection of rows i,i + 1 and columns j, j + 1. One sees that only in the top- left 
commuting square the transfers realising the same operation are between different pairs of entries. 
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We see that when considering non-unit transfers between pairs of adjacent columns and rows, a transfer 
that can be realised in a single step can become a transfer broken into several steps after the commutation. 
The opposite can also happen, and it suffices to consider the inverse operations in the diagram to obtain 
an example of that situation. 
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2 Crystal graphs 

§2. Crystal graphs. 

We shall now discuss some basic properties of crystal operations on matrices, starting with those that 
involve only horizontal or only vertical operations. So for the moment we assume that one of the binary 
and the integral cases is chosen, as well as one of the two orientations horizontal and vertical; then raising 
and lowering operations can be indicated simply by ej and /j, respectively. Doing so we essentially find 
the structure of crystal graphs of type A n , so the results in this section are generally known in some form. 
Our goal here is to give an brief overview of facts, formulated for matrices, and of terminology used. 

A crystal graph of matrices is a set closed under the application of operations e, and /, for i G N 
whenever these are defined, with matrices considered as vertices, linked by directed edges defined and 
labelled by those operations. A connected crystal graph will be called a crystal. Since for any nonzero 
matrix there is always some lowering operation that can be applied, such graphs are necessarily infinite. 
It will in some cases be more convenient to fix some n > and to restrict our attention to matrices with 
at most n rows (in case of vertical operations) or at most n columns (in case of horizontal operations), in 
which case the set of crystal operations is limited to and fi for i G [n — 1]; then crystals will be finite. 

Crystals have a complicated structure, but a number of properties can be established directly. 
Whether considering finite or infinite sets of crystal operations, it is easy to see that the crystals that 
occur do not depend, up to isomorphism, on whether one uses horizontal or vertical operations. In the 
case of integral matrices, transposition defines an isomorphism between the two structures on the set of 
all such matrices. In the binary case, rotation of matrices by a quarter turn counterclockwise defines an 
isomorphism from crystal graphs defined by vertical moves to crystal graphs defined by horizontal moves, 
but such a rotation can only be defined for a subset of binary matrices at a time. For any crystal for 
vertical moves, there is some to G N so that all its nonzero matrix entries are contained in the leftmost to 
columns, and an appropriate rotation maps these columns to the topmost m rows in reverse order. 

Using suitable rotations by a half turn, one can similarly show that when crystal operations ei and fi 
are limited to i e [n — 1] , the structure obtained from a crystal graph by interchanging the labels and 
fn-2-i is again a crystal graph of the same type. (It can be shown using other known properties that 
each crystal thus gives rise to a crystal isomorphic to itself, but that is not obvious at this point.) 

It is also true that crystals defined using binary matrices are isomorphic to those defined using integral 
matrices. Here the correspondence is rather less obvious than in the above cases, so we shall state it 
more formally. The statement basically follows from the fact that coplactic operations on semistandard 
tableaux can be defined either using the Kanji reading order (down columns, from right to left), or 
using the Semitic reading order (backwards by rows, from top to bottom); see [vLee3, proposition 3.2.1]. 
Translating the operations on tableaux in terms of matrices, using binary respectively integral encodings, 
gives crystal operations on matrices, and leads to the following formulation. 

2.1. Proposition. Let T G SST(A/k) be a semistandard tableau, and let its binary and integral encod- 
ings be M G M l2] and N G M, respectively. Then for to G N, e^M) is dchned if and only if e^(N) 
is defined, and if so, (M ) and e~ (N) are respectively the binary and integral encodings of one same 
tableau T' G SST(A/k). The same is true when and are replaced by and f^, respectively. 

This proposition immediately implies that the correspondence of M and N as encodings of the 
same semistandard tableau of shape X/k can be extended to a correspondence between the crystal for 
vertical moves that contains M and the crystal for horizontal transfers that contains N, so that these 
two structures are isomorphic. Note that a tableau may satisfy the tableau condition for many different 
shapes (indeed any tableau satisfies infinitely many tableau conditions, though not all of them are equally 
interesting), so for a crystal of one type, one obtains isomorphisms with many different crystals of the 
other type. This indicates that it is quite common to have distinct but isomorphic crystals of matrices; 
that observation could also have been made after commutation lemmas of the previous section. 

When we consider in a crystal graph only the edges labelled ej and fi for a single index i, then as we 
have seen, each vertex M is part of a finite ladder. The cancellations we considered in the previous section 
were defined, after selecting an appropriate ladder, by reversal of a lower portion of a ladder. It is natural 
to consider also the simpler operation of reversing each complete ladder for a fixed index i. This defines 
an involution on the set of all matrices, which in the case of vertical operations preserves all column sums 
while interchanging the row sums row(M)i and row(M)i + i, and in case of horizontal operations preserves 
all row sums while interchanging col(M)i and col(M)i + i. Now if one takes any crystal graph for which 
the multivariate generating series of its matrices by row sums (for vertical operations) or by column sums 
(for horizontal operations) exists (i.e., its coefficients remain finite), then that series will be invariant 
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2 Crystal graphs 



under all permutations of the variables; if the monomials of the series are of bounded degree, this means 
that the series is a symmetric function. (While the crystal graph of all binary or integral matrices is too 
large to have such a series, each individual crystal does indeed define a symmetric function.) 

This is similar to the situation for Bender-Knuth involutions, which can be used not only to de- 
fine cancellations, but also to show the symmetry of the weight generating series of skew semistandard 
tableaux, in other words of skew Schur functions. Indeed, the crystal graphs on the sets Tabl' 21 (A/k) or 
Tabl(A/«;) for some skew shape X/k are simple examples of graphs to which the symmetry result obtained 
by ladder reversal can be applied, and here the generating series is the skew Schur function s A / k LYn] in 
both cases. Although this result can also be obtained by using the simpler Bender-Knuth involutions, 
there is at least one sense in which it is "better" to use the ladder reversing involutions, namely that 
they actually define a symmetric group action of the underlying set of any crystal graph: the ladder 
reversing involutions satisfy the Coxeter relations of type A, whereas the Bender-Knuth involutions do 
not. (On the other hand ladder reversal does not respect the potentials for crystal operations at other 
indices, so the action is not by graph isomorphisms of any kind.) This symmetric group action means 
that one not only knows, for any crystal graph for vertical operations and any pair of compositions a, a' 
whose parts are related by some permutation, that there are as many matrices M in the crystal graph 
with row(M) =aas with row(M) = a' , but that one also has a natural bijection between those sets of 
matrices. 

We shall now state the above matters more formally. We first define the involutions generating the 
action, which come in four flavours: binary or integral, and horizontal or vertical. However, like for the 
crystal operations, there is no explicit distinction in the notation between the binary and the integral 
case, so only two symbols are used, o\ for the "vertical" involutions and for the "horizontal" ones. 

2.2. Definition. The symbols a] and cry for i, j € N denote involutions both on Ai l2] and on M., given 
in either case by the following expressions: 

(t\{M) = (e]) s (M) where S = n\{M) - n\{M) = row(M) m - row(M),, (24) 
C7~(M) = (ey) s (M) where S = nj{M) - nj{M) = co\(M) 3+1 - col(M)j. (25) 

and where (e\) s and {ey) s are to be interpreted as (//)~ 5 and {fy)~ s , respectively, when S < 0, and as 
the identity when S = 0. 

It is immediately obvious that one has the following identities, where Si denotes the transposition of 
i and i + 1, acting as usual on compositions by permuting their parts: 

row( aj(M)) = s t • row(M) row(a"(M)) = row(M) 

t and . (26) 

col(<7?(M)) = col(M) col(<7~(M)) = 8j ■ col(M) 

Moreover, when row(M)i = row(Af)i+i one has aj(M) — M, and similarly one has cry(M) — M when 
col(M)j = col(M)j + i; this contrasts with the Bender-Knuth involutions, which may send one tableau 
whose weight already has equal values at the positions being interchanged to another such tableau. It 
follows from the commutation lemmas that a] and a J commute for all i,j 6 N. 

2.3. Lemma. The operations u\ satisfy the Coxeter relations (crj) 2 = 1, <j\<j\ = <jj.aj when \k — i\ > 2, 
and ajaj +1 aj = crj +1 crjcrj +1 . Similar relations hold for the operations ay . 

Via the correspondence of crystal operations to coplactic operations on semistandard tableaux, this 
lemma is essentially equivalent to [Loth, theorem 5.6.3], but we shall give a proof at the end of our paper. 
The following theorem summarises the main conclusions that follow from the relations given. 

2.4. Theorem. There are two commuting actions of the group Soo, denoted by ir l (M) and n"(M), 
on each of M [2] and M; they are defined in both cases by s\{M) — <j\{M) and s~(M) = cr 4 "(M) for a 
generator Si of Soo (the transposition (i, Any operation -k 1 permutes the parts of row(M) by ir 
while leaving col(M) invariant, and if ir fixes row(M) then it 1 fixes M; similar statements hold for n". 

Proof. The only point that is not immediate is the final claim that tt 1 fixes M, which only follows directly 
from the observed properties of a] and ay in case the stabiliser of row(M) is generated by a subset of 
the generators Si (i.e., any equalities among the parts of row(M) occur in consecutive ranges). But that 
stabiliser is always conjugate to such a (parabolic) subgroup, which allows reduction to the given case. □ 
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We continue with some simpler but important observations. While every matrix admits some lowering 
operation, both vertically and horizontally, there are matrices that permit no raising operations in a given 
direction at all. Indeed it follows from propositions 1.3.5 and 1.4.5 that a M G M [2] or N G M admits no 
upwards operations if and only if Ms LR [2] (v / '(0)) where v — col(M) G V respectively N G Tabl(A/(0)) 
where A = col(M) G V, while the matrix admits no leftwards operations if and only if M G Tabl [21 (A/(0)) 
respectively N G LR(z//(0)), where A = row(M) and v = row(TV). The given values of A and v arc 
the only ones for which the respective conditions could hold, so these properties state that the fact that 
no raising operations are possible means that the matrix encodes, according to the case considered, a 
semistandard tableau of straight shape, or some (skew) Littlewood- Richardson tableau. 

Fixing one direction as before, if we start with a given matrix, and repeat the step of applying some 
raising operation that can be applied, until no such operation exists, then we end up with a matrix of 
the type just described. The process cannot go on indefinitely since at every step we move a unit to a 
row or column with a lower number, which decreases an obvious statistic (the sum over all units of their 
row respectively column number). The matrix found at the end of the process will be said to be obtained 
from the original matrix by exhausting raising operations (of the type considered) . We shall show later 
that the final matrix found does not depend on the choices made during the exhaustion process (in other 
words, the process is confluent). This is an important and quite nontrivial statement; it implies that every 
crystal contains a unique vertex in which no raising operations are possible. For the moment however it 
suffices to know that some final matrix can always be reached in this way. 

Similar properties are known to hold for jeu de taquin on skew semistandard tableaux, which can 
be used to "rectify" a skew tableau to one of straight shape by a nondeterministic procedure, but whose 
final result is independent of the choices made. In fact there is a close relation between that game and 
crystal operations on matrices encoding the tableaux. In fact it is known, see [vLee3, theorem 3.4.2], 
that jeu de taquin slides applied to a tableau T correspond to sequences of coplactic operations applied 
to a its "companion tableau" , which is a tableau whose integral encoding is the transpose of the integral 
encoding of T. Since coplactic operations on tableaux correspond to horizontal crystal operations on 
their integral encodings, this means that jeu de taquin slides correspond to sequences of vertical crystal 
operations. A correspondence exists also for binary encodings. More precisely, one has the following. 

2.5. Proposition. Let T G SST(A/k) and let T G SST(A'/k') be obtained from T by an inward jeu de 
taquin slide starting at the corner (io,jo) of k and ending in the corner (zi,ji) of A. The their respective 
binary encodings M, M' are related by M' = e^_ 1 (- • • (e^ +1 (e^(M)) • • ■)), and their respective integral 
encodings N, N' are related by N' = et i _ 1 (- • • (e? o+1 (eJ o (Nj) ■■■)). 

Proof. In the description of a jeu de taquin slide in terms of a skew tableau with an empty square into 
which entries as successively slid, one can at each point associate both a binary and an integral encoding 
matrix, which count the entries in columns respectively rows in the usual fashion (ignoring the empty 
square). It is a straightforward verification that under this correspondence each leftward slide of a tableau 
entry results in applying a leftward move in corresponding pair of columns of the binary matrix while 
leaving the integral matrix unchanged, and that each upward slide of a tableau entry results in applying 
an upward unit transfer in corresponding pair of rows of the integral matrix while leaving the binary 
matrix unchanged. Indeed, the inequalities in definition 1.3.1b reflect the fact that the entries above and 
below a horizontally sliding entry, in the same pair of columns, satisfy weak increase along rows, while the 
inequalities in definition 1.4.1a reflect the fact that the entries to the left and right of a vertically sliding 
entry, and in the same rows, satisfy strict increase down columns. The proposition immediately follows. □ 

Thus one can see that if a skew tableau T has rectification P, then the binary encoding of P can be 
obtained from that of T by exhausting leftward moves, while the integral encoding of P can be obtained 
from that of T by exhausting upward transfers. However this reasoning does not yet imply that these are 
the only possible results of such exhaustion of raising operations, even if one admits the uniqueness of the 
rectification of T, since jeu de taquin slides correspond only to specific composition of raising operations, 
and there might exist transitions that can be realised by raising operations but not by jeu de taquin. 

Since the raising operations on matrices in the direction perpendicular to the one just considered 
correspond to coplactic operations on tableaux, one easily recognises a relation between the commutation 
lemmas of the previous section and the commutation of coplactic operations with jeu de taquin, as stated 
for instance in [vLee3, theorem 3.3.1]. In the next section we shall elaborate on this relation, developing 
a theory for matrices that is similar (but more symmetric) than the one for semistandard tableaux. 
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§3. Double crystals. 

We shall now return to considering at the same time vertical and horizontal operations. A set of all 
matrices that can be obtained from a given one by using both vertical and horizontal crystal operations, 
provided with the structure of labelled directed graph defined by those operations, will be called a double 
crystal. Due to the commutation lemmas, double crystals will be much easier to study than ordinary 
crystals. We shall see that they can be used to define in a very natural way correspondences quite similar 
to the RSK correspondence and to its dual correspondence that was (also) defined in [Knu] , respectively. 
In fact, the correspondence that is most naturally obtained for integral matrices is the Burge variant of 
the RSK correspondence, described in [Fult, A. 4.1]. 

3.1. Knuth correspondences for binary and integral matrices. 

We have just seen that by exhausting either vertical or horizontal raising operations, one can transform any 
integral or binary matrix into one that encodes, according to the case considered, a straight semistandard 
tableau or a Littlcwood-Richardson tableau. Now by exhausting both vertical and horizontal raising 
operations, one can obviously transform any integral or binary matrix into one that encodes a Littlewood- 
Richardson tableau of straight shape. It is well known that the possibilities for such tableaux are extremely 
limited: in their display, any row i can contain only entries i, whence such tableaux are completely 
determined by either their shape or their weight, which moreover are the same partition A. One easily 
sees that the binary and integral encodings of such tableaux are the matrices described by the following 
definition. 

3.1.1. Definition. Let A E V. We denote by Diag [21 (A) G M l2] the binary matrix E [A]]) JjGN 
whose bits 1' occur precisely in the shape of the Young diagram [A]. We denote by Diag(A) E M the 
integral matrix = ([i = j] \i)i,j£N, which is diagonal with the parts of A as diagonal entries. 

One immediately sees that this binary matrix satisfies row(Diag [2] (A)) = A and col(Diag [21 (A)) = A', 
while this integral matrix satisfies row(Diag(A)) = col(Diag(A)) = A. From the consideration of tableaux 
one may deduce that these matrices are precisely the ones for which all vertical and horizontal raising 
operations are exhausted, but it is easy and worth while to prove this directly. 

3.1.2. Lemma. A binary matrix M E M [2] satisfies Vi:nJ(M) = and Vj:nJ~(M) = if and only if 
M = Diag [21 (A) for some A E V. Similarly an integral matrix A E Ai [2] satisfies Vi:nJ(A) = and 
Vj: nj-(N) = if and only if N = Diag(A) for some XeV. 

Proof. One easily sees that no raising operations are possible for Diag [2I (A) or for Diag(A), since for 
Diag [21 (A) there are not even candidates for upward of leftward moves, while for A — Diag(A) and all i 
the units of N i+ i^ + i are blocked by those of A, ; j for upward or leftward transfers. This takes care of the 
"if" parts. For the "only if" parts we shall apply induction on the sum of all matrix entries, the case of 
null matrices being obvious (with A = (0)). 

Suppose that Mi:n\{M) — and Vj:n^(M) = 0, and let (k,l) be the indices of the bit Mfe ; ; = 1 
for which / is maximal, and k is minimal for that /. If k > one would have nI_ 1 (M) > 0, so k = 0. 
But then also Moj = 1 for all j < I, since otherwise the largest j violating this would give nj{M) > 0. 
Let M' = (Mi + ij)i je N be the matrix obtained by removing the topmost row Mo. The hypothesis 
Uj-(M) = implies nf(M') = for j < I since (M ,j, M ,j+i) = (1, 1); one also has ny(M') = for 
j > I, for lack of any bits M' t = 1. Therefore M' satisfies the hypotheses of the lemma, and by 
induction one has M' = Diag [2) (A') for some A' E V, with obviously X' < I. Then M = Diag [21 (A), where 
A is the partition obtained by prefixing / to the parts of A' (formally: Ao = / and Aj+i = A^ for i E N). 

Now suppose that Vi:nj(N) — and \/j:ny(N) — (the integral case). The conditions n\{N) = 
and Uj~(N) = respectively imply that A i+ i !0 = and A 0j +i = 0. But one cannot have row(A) = 0, 
since then for the smallest index i such that row(A)i + i > one would have n](N) = row(A)i + i > 0. 
Therefore Ao t o > and we can apply induction to the matrix A' = (Aj+ij+i^jgN obtained by removing 
the initial row and column from A, since clearly n-(A') = nJ +1 (A) = and n^(N') = nJ +1 (N) = 0. 
Then A' = Diag(A') for some A' E V, and Ao.o > A follows either from n 1 (N) = or from 
Uq'(N) = 0. So one has A = Diag(A), where A E V is obtained by prefixing A 0j o to the parts of A'. □ 

Now that we have characterised the matrices at which the process of applying raising operations can 
stop, and with the commutation lemmas at our disposal, the general structure of double crystals can be 
easily analysed, as follows. 
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3. 1 Knuth correspondences for binary and integral matrices 

Fix a matrix M G M l2] or M G A4. Suppose that R r is some sequence of raising operations e\ for 
varying i G N that exhausts such operations when applied to M, in other words such that the resulting 
matrix i? T (M) satisfies Vz: n\(W {M)) = 0. Similarly let R*~ be some sequence of raising operations 
e'j for varying j G N that exhausts such operations when applied to M; then the resulting matrix 
R^(M) satisfies \/j:n^(R^(M)) = 0. By lemma 1.3.7 or lemma 1.4.7, the expressions i?^(i? T (M)) and 
P T (P^(M)) are both defined and they designate the same matrix, which we shall call TV. By proposition 
1.3.6 or 1.4.6, one has Vz:nJ(iV) = and Vj:rij~(N) — 0, and so lemma 3.1.2 states that N = Diag [21 (A) 
or N — Diag(A) (according to the case considered) for some A G V . Since R*~ preserves row sums one has 
row(P T (M)) = row(N), which composition equals A both in the binary and in the integral case; similarly 
i? T preserves column sums, so one has col(P^(M)) = co\(N), which equals A in the integral case and A* 
in the binary case. Thus either of the values row(P T (M)) or col(P*~(M)) completely determines A. 

Now let S* T and be other sequences of upward respectively leftward crystal operations that can be 
applied to M and that exhaust such operations. Then the matrices 5^(i? T (M)) and S r (R*~ (M)) satisfy 
the conditions of lemma 3.1.2, and since their row sums respectively their column sums are the same 
as those of N, they must both be equal to N. From R*~ we can form the "inverse" operation R~" by 
composing the lowering operations /~* corresponding to the operations ej~ used, in the opposite order. 
Applying R~ to the equation R^(W(M)) = N = S T (P^(M)) = i?^(S* T (M )), one finds W{M) = S T (M), 
and one similarly finds R"(M) = S"(M). So the matrices obtained from M by exhausting vertical or 
horizontal crystal operations are uniquely defined. The matrix N obtained by exhausting both types of 
operations is also uniquely determined by M, and we shall call N the normal form of M. 

Moreover from the knowledge of i? T (M) and R^(M) one can uniquely reconstruct M. To state this 
more precisely, let P and Q be matrices of the same type (binary or integral) that satisfy Mi: n\(P) = 
and \/j:n,j~(Q) = (which implies that row(P) and col(Q) are partitions), and that satisfy moreover 
row(P) = co^Q)' in the binary case, or row(P) = col(Q) in the integral case. Let R^(P) and P T (Q) 
be obtained by from P by exhausting leftward crystal operations, respectively from Q by exhausting 
upwards crystal operations (as before R*~ and P T denote particular sequences of operations that realise 
these transformations). Both matrices satisfy the conditions of lemma 3.1.2, and by the hypothesis on 
column and row sums, the partition A in the conclusion of that lemma is the same in both cases; therefore 
R"~(P) = P T (Q). As before we can form inverse operations R~" and R l of R*~ and of i? T , respectively. 
We successively apply these inverse operations to the equation R^(P) = i? T (Q), and by commutation of 
R~* and R l (which follows that of R*~ and i? T ) one finds R l (P) = R~"(Q). In the members of this final 
equation we have found a matrix M for which i? T (M) = P and R^(M) = Q. 

To show uniqueness of M, suppose that M' is another matrix, and that 5 T and are other 
sequences of upwards respectively leftwards crystal operations for which P = S T (M') and Q = S^(M'). 
Then (P) is the normal form of P, and therefore equal to R*~ (P). From this equality and P — S r (M') 
one deduces by commutation that S r (S^(M')) = S r (R~(M')), and since S^(M') = Q = R~(M) this 
gives 5" T (i?^(M)) = 5 T (i?^(M')). But 5" T and R^ arc invcrtiblc operations, so it follows that M = M' . 

The above reasoning shows the extraordinary usefulness of the commutation lemmas. But despite 
the relative ease with which the conclusions are reached, one should not forget that each application of a 
crystal operation has to be justified; for instance in the final part we could talk about S~{P) only because 
P = S r (M') and by assumption S^(M') is defined. Also it should not be inferred from S*~(P) = J2*~(P) 
that and R*~ are equivalent sequences of crystal operations, or even that they can be applied to the 
same set of matrices. We now state more formally the statement that follows from the argument given. 

3.1.3. Theorem, (binary and integral decomposition theorem) 

(1) There is a bijective correspondence between binary matrices M G A4 l2] on one hand, and pairs 
P,Q G M m of such matrices satisfying Mi : n\ (P) — and Vj:n^(Q) — and row(P) = col(Q) ' on the 
other hand, which is determined by the condition that P can be obtained from M by a finite sequence 
of operations chosen from {ej | i G N } and Q can be obtained from M by applying a finite sequence 
of operations chosen from {ej \ j G N }. In particular col(P) = col(M) and row(Q) = row(M). 

(2) There is a bijective correspondence between integral matrices M G M on one hand, and pairs 
P,Q G Ai of such matrices satisfying Vi: n\(P) = and Vj: n*~(Q) = and row(P) = col(Q) on the 
other hand, which is determined by the condition that P can be obtained from M by a finite sequence 
of operations chosen from { e ■ | i G N } and Q can be obtained from M by applying a finite sequence 
of operations chosen from { e*j~ \ j G N }. In particular col(P) = col(M) and row(Q) = row(M). 

These correspondences can be explicitly computed by the constructions given above. □ 
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3.2 Relation with jeu de taquin and with Robinson's correspondence 

3.1.4. Corollary. The crystal defined by any binary or integral matrix, and by one type (vertical or 
horizontal) of crystal operations, contains a unique vertex at which the potentials for all raising operations 
are zero. Moreover, the crystal is determined up to isomorphism by the potentials for the lowering 
operations at that vertex. 

Proof. We have seen that exhausting upward or leftward crystal operations in all possible ways from a 
given initial matrix always leads to the same final matrix (the matrix P respectively Q of theorem 3.1.3), 
where the potentials for all raising operations are zero. Then whenever two matrices M, M' are linked 
by a crystal operation in the chosen direction, this final matrix will be the same whether M or M' is 
taken as point of departure. Consequently, every matrix in the crystal considered leads to the same final 
matrix, which is the unique matrix in the crystal where the potentials for all raising operations are zero. 

For the final part of the claim, consider the common normal form N — Diag [21 (A) or N = Diag(A) 
of each of M, P, and Q in theorem 3.1.3. Since horizontal crystal operations induce isomorphisms of 
crystals defined by vertical crystal operations (due to their commutation) and vice versa, the crystal 
for vertical operations containing P is isomorphic to the one for the same operations containing N, 
and similarly for the crystal for horizontal operations containing Q; it will suffice to show that A can 
be deduced from the sequences of potentials for the appropriate lowering operations at P respectively 
from those at Q. Now nj(P) = n iC^0 = ^* ~~ f° r au * m both the binary and the integral case, 
while nj{Q) = nj*(N) for all j, which equals Aj — Aj +1 in the binary case, or Xj — Aj+i in the integral 
case. Since A can be reconstructed from the sequence of differences its consecutive parts or of those 
of A' (for instance = J2i>k^i f° r au k £ N, where di = Aj — Aj+i), this completes the proof. □ 

Part (2) of the theorem is very similar to the statement of bijectivity of the RSK correspondence 
(see for instance [Stan, theorem 7.11.5]), since the matrices P and Q* both lie in Tabl(A/(0)) where A is 
the partition row(P) = col(Q), whence they encode straight semistandard tableaux of equal shape A; the 
statements col(P) = col(M) and row(Q) = row(M) give the usual relation between the weights of these 
tableaux and the column and row sums of M. However, the construction of the bijection is quite different, 
and as we shall see, the bijection itself is does not correspond directly to the RSK correspondence either. 
Yet it is appropriate to call the bijection of part (2) a "Knuth correspondence" . This is a generic term for 
a large class of generalisations of the RSK correspondence that was used in [Fom] and in [vLee4] , without 
being formally defined. It actually refers to the result of the type of construction defined in [Fom] using 
growth diagrams; that construction is rather different from the one above, but as we shall see below, the 
bijection considered does allow an alternative construction using growth diagrams. 

As for the bijection of Part 1, it has some resemblance to dual RSK correspondence, the second 
bijection that was defined together with the usual RSK correspondence in [Knu]. That correspondence 
is a bijection between binary matrices and pairs of straight tableaux that are less symmetric than those 
associated to integral matrices: one can define them to be both semistandard but of transpose shapes 
(as in the original paper), or one can define one of them to be semistandard and the other transpose 
semistandard, in which case their shapes will be equal (as in [Stan, §7.14]); in fact at the end of the 
original construction one of the tableaux is transposed to make both semistandard. In our bijection 
however, while Q is the binary encoding of a straight semistandard tableau of shape A = col(Q)', it is 
not entirely natural to associate a straight semistandard or transpose semistandard tableau to P. One 
could rotate P a quarter turn counterclockwise to get the binary encoding of a semistandard tableau of 
shape row(P)' = A', but that matrix, and the weight of the tableau it encodes, depend on the number 
of columns of the matrix that is rotated, which should be large enough to contain all bits T' of P. The 
type of straight tableau most naturally associated to P is a reverse transpose semistandard tableau of 
shape row(P) = A, namely the tableau given by the sequence of partitions (A^); £ n where A^ = X^>2 Pj- 

3.2. Relation with jeu de taquin and with Robinson's correspondence. 

Let us give an example of the bijections of theorem 3.1.3. We could take arbitrary matrices M in both 
cases, but it will be instructive to use the binary and integral encodings of a same semistandard tableau T. 
We shall take 
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3.2 Relation with jeu de taquin and with Robinson's correspondence 
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= M, 



where the composite raising operations are for instance R*~ = ef e§ eg" ej" ej~ e$ ej~ ejr e^ e 1 ^ eg" and 
R) = (el) 3 (e\f (e;y (e\Y {elf (etf (e{f (4) 3 (ej) 3 (e\f (e^f (e\f (4) 4 (ej) 4 e\ e\ ej e\ ej; as usual 
operations are applied from right to left, and we have systematically chosen to climb the ladder with the 
lowest possible index first, and to climb it completely (so without switching ladders when a lower index 
ladder becomes available halfway up). As an illustration of theorem 3.1.3(2), taking for M the integral 
encoding of T, one has the diagram 
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M, 



where the composite raising operations are for instance (using the same rule to choose ladders) given by 
R~ = (eD 3 i^f i^) 7 (eZ) 7 (e^) 7 {e^ f (ef) 3 fe-) 3 (e^f (^) 4 K) 2 (^ f (^Y i^ f ef e$ efr ef e» 
and W=4 (e^ e\ (e\f e\ (e\f (ej) 2 . 

Comparing (28) and (29), one may observe that in both cases the normal form TV is parametrised by 
the same partition A = (8, 8, 5, 3, 1). This can be explained using proposition 2.1, which implies that the 
same composite raising operations can be taken for R r in the binary case as for R*~ in the integral case, 
up to the replacement of e\ by e"~ ; therefore the binary matrix P and the integral matrix Q both encode 
the same tableau L, of the shape (9, 8, 5, 5, 3)/(4, 1) of M. The raising operations in the given direction 
are exhausted, so this is a Littlewood-Richardson tableau, namely 
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(30) 



The partition wt(i) = (8, 8, 5, 3, 1) is equal to row(P) in the binary case and to col(<5) in the integral 
case, and in both cases it is therefore also equal to the partition A = row(N) that parametrises N. One 
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3.3 Identifying the factors in the decomposition 



can also read off A from the matrix Q in the binary case, or from the matrix P in the integral case, namely 
as A = col(Q)' respectively as A = row(P). Those matrices are in fact the binary and integral encodings 
of the rectification 5* of T, which has shape A, since exhausting inward jeu de taquin slides applied to T 
translates into particular sequences that exhaust leftward respectively upward crystal operations when 
applied to M, by proposition 2.5. That rectification is easily computed to be 
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(31) 



and its binary and integral encodings are indeed the binary matrix Q above and the integral matrix P, 
respectively. The matrices L and S are the ones associated to T under Robinson's bijection, as is described 
in [vLee3, corollary 3.3.4]; in the notation of that reference one has (L, S) = 1Z(T). 

Thus in a certain sense both bijections of theorem 3.1.3 model Robinson's bijection. However that 
point of view depends on interpreting matrices as encodings of tableaux, which might not always be 
natural. In any case one should be aware that one matrix can simultaneously encode tableaux of many 
different shapes. This point of view also gives somewhat incomplete information about double crystals, 
because even though the exhaustion of appropriate raising operations can be realised using jeu de taquin, 
it is not easy to do the same with individual raising operations. 

3.3. Identifying the factors in the decomposition. 

Since any matrix can easily be seen to encode a skew semistandard tableau of an appropriately chosen 
shape, the relation established above gives a complete characterisation of the bijections in theorem 3.1.3 
that is independent of crystal operations. However, it describes the matrices P and Q in rather different 
terms, which seems unnatural given the characterisation of theorem itself. As we have observed above, 
the theorem in fact suggests a close relationship with the well known RSK correspondence (and its dual) , 
in which the two straight tableaux associated to a matrix play a much more comparable role. We shall 
now proceed to characterise the bijections of our theorem in the language of those correspondences. 

It is well known that the "P-tablcau" of the pair that is associated under the RSK correspondence 
to an integral matrix M can be computed using rectification by jeu de taquin; since the matrix P of 
3.1.3(2) can also be determined using rectification, this provides a useful starting point to relate the two 
correspondences. In more detail, the mentioned P-tableau can be found by rectifying a semistandard 
tableau H whose shape is a horizontal strip, and whose "rows", taken from left to right, successively 
contain words (each of which is necessarily weakly increasing) whose weights given by the rows of M, 
from top to bottom. This observation was first made, for the case of Schensted correspondence, in [Schu.2], 
see also [Stan, Al.2.6]; the generalisation to the RSK correspondence follows directly from its definition. 
For instance, the initial skew tableau corresponding to the matrix 
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M = 2 2 1 would be H = 
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Note that by filling the horizontal strip from left to right, the vertical order of rows is interchanged: the 
topmost row of M determines the bottom row of H. Therefore the relation between M and H is not 
that of a matrix encoding a tableau; rather the integral encoding of H is the matrix formed by reversing 
the order of the rows of M. As for the matrix P of 3.1.3(2), we have seen above that it is the integral 
encoding of the rectification S of any tableau T of which M is the integral encoding. There are many 
possibilities for T, all of which differ only by horizontal shifts of their rows and which therefore have the 
same rectification. It is in particular always possible to take for T a tableau whose shape is a horizontal 
strip. Comparing with what we said for the RSK correspondence, and using the fact that both for that 
correspondence and for the one of theorem 3.1.3(2), transposing M results in interchanging P and Q (or 
actually P and Q l for our theorem), we arrive at the following description. 

3.3.1. Proposition. Let M € M, which we view as a some sufficiently large Gnite rectangular matrix. 
The corresponding matrices P,Q 6 M of theorem 3.1.3(2) can be found as follows: P is the integral 
encoding of the P-symbol under the RSK correspondence of the matrix obtained from M by reversing 
the order of its rows, and Q t is the integral encoding of the Q-symbol under the RSK correspondence 
of the matrix obtained from M by reversing the order of its columns. □ 
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3.3 Identifying the factors in the decomposition 

Note that reversal of the order of rows or columns is not a well defined operation on AA, but the 
resulting ambiguity about the matrices of which the P-symbol and the Q-symbol are to be taken is 
harmless, since it only affects the other symbol associated to each matrix, the one that is not used in the 
proposition. For instance the matrix obtained by the reversal of the order of rows is only determined up 
to a vertical shift, but this only induces a shift in the entries of the associated Q-symbol, which is unused. 

The description of the proposition can be reformulated as follows: P and Q l are the integral encodings 
of the insertion tableau for the traversal of M in reverse Semitic reading order (by rows, read from left to 
right, ordered from bottom to top), row-inserting at each entry Mij copies of the number j, respectively 
of the recording tableau for a similar insertion, but of negated values and traversing M in the opposite 
(Semitic) order. The latter tableau can also be described either as the Schiitzenberger dual of the 
recording tableau for the initial insertion process, or as the insertion tableau for a traversal of M in Kanji 
reading order, where at each entry Mij copies of the number i are row-inserted, but in all these cases 
this involves two separate applications of an RSK or Schiitzenberger algorithm. 

However the pair of tableaux of which P and Q t are encodings is in fact the pair associated to M 
under the "Burge correspondence" . This is a variation of the RSK correspondence that was first defined 
(for a somewhat specialised case) in [Bur], and which is discussed in detail in [Fult, A. 4.1] and in [vLee4, 
§3.2]. To be more specific, this gives us the following description. 

3.3.2. Proposition. The matrices P, Q e M corresponding to M G M under the bijection of theo- 
rem 3.1.3(2) are respectively the integral encoding of the P-symbol and the transpose of the integral 
encoding of the Q-symbol that are associated to M under the Burge correspondence. This P-symbol can 
be obtained by traversing M in the Semitic reading order (by rows from right to left) and successively 
column inserting, into an initially empty tableau, for each position (i,j) encountered Mij copies of the 
number j; the Q-symbol is the recording tableau (A^)i e N for this procedure, where is the shape of 
the tableau under construction before the traversal of row i of M . 

The sequence of numbers that is being inserted is also the Semitic reading of any tableau of which 
M is the integral encoding; for instance for the matrix M in (29) this is the sequence 5, 6, 6; 2, 3, 4, 6, 6; 
1, 1, 2, 4, 5; 0, 1, 3, 4, 6, 6, 6; 0, 2, 4, 5, 5, read from right to left. (For column insertion it is best to imagine 
the entries being inserted as coming from the right, so writing this sequence backwards is actually quite 
natural; for instance, in the notation of [Fult] the insertion tableau may be written as 5 — > (6 — > (6 — > 
(2 — > • • • — > (2 — > (4 — > (5 — > (5)))) • • •))), where we added redundant parentheses to stress the right- 
associativity of the operator '—>'.) Since column insertion can be simulated by jeu de taquin just like 
row insertion can, with as only difference that the entries inserted are initially arranged into a tableau 
according to the Semitic reading order, it is clear that the tableau obtained by column inserting the 
sequence given is equal to the rectification S of the tableau L given in (31). This argument justifies the 
"P"-part of the proposition above; to justify its "Q"-part, the easiest argument is to recall that the Burge 
correspondence enjoys the same symmetry property as the RSK correspondence, and that the bijection 
of theorem 3.1.3(2) has a similar symmetry 

We illustrate this computation for the integral matrix M of (29). We are performing column insertion, 
and the display below should be read from right to left; it shows the state of the insertion process after 
each row of M is processed. The entries being inserted are written above the arrows between the tableaux; 
they too should be read from right to left to get their order of insertion. 

| | 1 | 2 | 4 | 4 | 5 | 5| | | | 1 | 2 | 4 | 4 | 5 | 5| 

1 1 3 4 5 6 1 6 1 6 I 5.6.6 1 3 4 5 §j6jM 2,3.4,6.6 

2 2 1 6761 * 2 2 4 1 6 | < 

3 e I e I 1 3 1 6 1 

5] 

Indeed the final tableau coincides with S of (31) that is encoded by P of (29), and the sequence of shapes 
(0) v- (5) v- (7, 5) v- (8, 7, 2) v- (8, 8, 4, 2) v- (8, 8, 5, 3, 1) defines the semistandard Young tableau 
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encoded by Q t ; it is the companion Young tableau of the Littlewood- Richardson tableau L of (30). 
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3.3 Identifying the factors in the decomposition 



For the bijection of 3.1.3(1), one can use a similar argument to identify the matrix Q, which is the 
binary encoding of the rectification of any tableau T encoded by the binary tableau M. There always 
exists such a tableau T whose shape is a vertical strip, and as in the case of a horizontal strip, the 
rectification can then be computed as the insertion tableau either for the column insertion of the entries 
taken from top right to bottom left, of for the row insertion of the entries taken in the opposite order. 
In either case the shape of the tableau grows by a vertical strip after the insertion of the entries for 
each complete column of T, which defines a transpose semistandard recording tableau; since its vertical 
strips are most naturally indexed by the numbers of the columns of T traversed, the recording tableau 
for the right-to-left column insertion is actually a reverse tableau. As we have observed above, the 
binary matrix P of 3.1.3(1) corresponds to a reverse transpose semistandard tableau, and this is in fact 
the recording tableau of the top right to bottom left column insertion of the entries of T (for the row 
insertion of the entries in the opposite order, it is the Schiitzenberger dual of the recording tableau). One 
thus has the following description of the bijection between M and (P, Q) of 3.1.3(1). 

3.3.3. Proposition. The matrices P,Q € Ai l2] corresponding to M s Ai l2] under the bijection of 
theorem 3.1.3(1) can be computed as follows. The matrix Q is the binary encoding of the semistandard 
Young tableau S obtained by traversing the matrix M in the Kanji reading order (down columns, from 
right to left), and successively column inserting, into an initially empty tableau, for each bit Mij equal to 1 
an entry i. The matrix P determines the reverse transpose semistandard Young tableau R = (X^)jeN 
that is the recording tableau for this insertion process: X^ = J2i>j Pi JS ^ e shape of the tableau under 
construction after the the traversal of column j of M. 

We illustrate this computation for the binary matrix M of our running example. As before we display 
successive stages of column insertion from right to left; the strictly increasing sequence of numbers inserted 
at each rightward arrow is written from top to bottom, as it is a column of the tableau T of (27). 
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Indeed the final tableau coincides with the tableau S of (31) of which the matrix Q of (28) is the binary 
encoding, and the reverse transpose semistandard Young tableau R given by the sequence of shapes 
(8, 8, 5, 3, 1) -, (7, 7, 4, 3, 1) -, (6, 6, 3, 2, 1) -r (5, 5, 2, 2) -, (5, 4, 1, 1) (4, 3) -r (3, 2) -r (2, 1) -r (1) -r 
(0) corresponds to the binary matrix P of (28). This tableau can be displayed by filling for each i the 
squares of the vertical strip [X^ / \( l+1>> ] with entries i, giving 
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(35) 



We shall finally compare the bijection of 3.1.3(1) with the dual RSK correspondence. One should 
first of all realise that, unlike for the ordinary RSK and Burge correspondences, there is no universal 
agreement of which pair of tableaux corresponds to a given binary matrix (versions differ by transposition 
of one or both tableaux). We shall use the form in [Stan], where the entries coming from a conventional 
reading of M are row-inserted into a transpose semistandard Young tableau, and the recording tableau 
is semistandard. Note that the following statement directly relates to the tableaux coming from the dual 
RSK correspondence applied to M (even if a Schiitzenberger dual is involved), whereas for the integer 
correspondence we could not use the tableaux corresponding to M itself under the RSK correspondence. 

3.3.4. Proposition. The matrices P,Q £ Ai l2] corresponding to M £ Ai l2] under the bijection of 
theorem 3.1.3(1) are related as follows to the pair (R*,S) that is associated to M under the dual RSK 
correspondence as described in [Stan, theorem 7.14.2]: the transpose semistandard Young tableau R* is 
the Schiitzenberger dual of the reverse transpose semistandard Young tableau R determined by P, and 
S is the semistandard Young tableau encoded by Q. 
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3.4 Relation with decomposition of pictures 

3.4. Relation with decomposition of pictures. 

We have established relations between the bijections of theorem 3.1.3 and known correspondences defined 
by various insertion procedures, but the descriptions could use some clarification. We have justified the 
correspondences with insertion tableaux via their relation to jeu de taquin, but the correspondences with 
the recording tableau are only justified using symmetry properties for the integral case, and only by an 
example in the binary case. We shall now present a correspondence that explains both factors in a more 
symmetric way, and allows to transparently justify the descriptions of the previous subsection. 

If an integral matrix M satisfies M G Tabl(A/fc), it encodes a semistandard tableau T of shape X/k, 
and if it also satisfies M G LR(v/u), then M' encodes a semistandard tableau T of shape v/u, which is 
a companion tableau of T. The pair (T, T) can be represented in a symmetric manner by the notion of a 
picture X/k — > v/u, which designates a particular kind of bijection between skew diagrams [X/k] — > [v / u], 
in a way that we shall now describe. 

In fact two different flavours of pictures can be defined, both of which have the property that the set 
of pictures X/k — > v/u is in bijection with Tabl(A/«;) C\LR(v/u), and both of which give rise to essentially 
the same theory. Our reference for the theory of pictures will be [vLee2] , but we shall define the opposite 
flavour of pictures than was used there, since that simplifies the correspondences that are of interest here. 
This means that we must adapt the results we cite to our situation; in fact the form of these results 
changes little between the two flavours of pictures, except for explicit references to the Schensted row 
insertion procedure (for instance in theorem 3.2.1 and lemma 3.3.2 of [vLee2]), which must be replaced 
by column insertion. The pictures used in the present paper are defined as follows (the difference with 
the other flavour of pictures is the use of rather than its opposite partial ordering). 

3.4.1. Definitions. 

a. The partial orderings '< N ' and '</, ' on N 2 are defined by < x (k, I) if and only if i < k and 
j < I, respectively by <y (k, I) if and only if i < k and j > I. 

b. For skew shapes X/k and v/u, a picture X/k — ► v/u is a bijection f: [X/k] — ► [v/u] between their 
diagrams, such that s < N t implies f(s) <y f(t) and f(s) < x f(t) implies s <y t for all s,t G [X/k]. 
The set of pictures X/k — > v/u is denoted by Pic(A/ k, v/u). 

c. For a picture f G Pic(X/k,v/u), the matrices Int(/) G Mx-K.y-ti and Bin(/) G M [ * ] _ are 
respectively given by Int(/)j, fe = #{ (j, I) £ N 2 G [X/k], (k, I) G [v/u], f(i,j) = (k, I) } and by 
Bin(/) fcJ =#{(*, G N 2 | (i,j) G [X/K],(k,l) G [v/u],f(i,j) = (k,l)}. 

The expressions in part c can be described in words as follows: Int(/)j ; fe counts the number of squares 
of row i that the picture / maps to row k, while Bin(f)kj counts the number of squares of column j 
that / maps to row k. 

From the form of part b it is clear that the inverse of a picture X/k — > v/u is a picture v/u — ► X/k. 
The point of part c of the definition is the following fact. 

3.4.2. Proposition. For any pair of skew shapes X/k and v/u, the operations Int and Bin define 
bijections from Pic(A//c, v/u) respectively to Tabl(A/«) n ~LR(v/u) and to Tabl [21 (A/«;) n LR [21 (^/^i). 

Proof. For a given matrix M, a picture / with Int(/) = M or Bin(/) = M is straightforwardly 
constructed, defining images of the squares of [X/k] in the Semitic (in the integral case) or the Kanji (in 
the binary case) order. The conditions of 3.4.1b allow only one candidate image for each square, and one 
checks that the tableau and Littlewood-Richardson conditions for M assure that this is a valid choice. □ 
When v/u is a horizontal strip, the conditions M G LR(v/u) and M G LR l2] (v/n) just require 
col(M) = v — u respectively row(M) = v — u, so in this case Pic(A//t, v/u) is in bijection with the set 
of tableaux in SST(A/k) of weight v — u. Moreover, if / so corresponds to a tableau T, then Int(/) and 
Bin(/) are respectively the integral and binary encodings of T. Thus we find our encoding of semistandard 
tableaux by integral or binary matrices as a special case of representing pictures by matrices using Int(/) 
or Bin(/). Similarly, if v/u is arbitrary but X/k is a horizontal strip, then the pictures in Pic(A/fc, v/u) 
are in bijection, via their inverse picture, with tableaux in SST(v/u) of weight A — k. In the latter case 
Int(/) = Int(/ -1 )' G LR(v/u) is the transpose of the integral encoding of the tableau thus corresponding 
to /; however there is in this case no similar interpretation for Bin(/) (which is in fact a rather sparse 
matrix, having column sums of at most 1). On the other hand, if X/k is a vertical strip rather than a 
horizontal one, then it is most practical to consider Bin(/) instead of Int(/): in this case Pic(A/«;, v/u) 
is in bijection with the set of reverse transpose semistandard tableaux of shape v/u and weight A* — k 1 , 
and their correspondence to the matrix Bin(/) E LR l2] (v/u) is exactly the one that we used before. 
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3.4 Relation with decomposition of pictures 



A central result concerning pictures is the decomposition theorem [vLee2, theorem 3.2.1], which 
generalises the Robinson-Schensted and RSK correspondences; it gives a bijection between Pic(A/«, v / \i) 
and the union over all 7r £ V of Cartesian products Pic(7r, v/[i) xPic(A//c, tt). This bijection can be defined 
by a version of the Schensted algorithm using (for our notion of pictures) column insertion, translated 
into the language of pictures, as described in [vLee2, lemma 3.3.2]. It follows that for specific shapes X/k 
and v/n, this decomposition specialises under the appropriate correspondences to the column insertion 
Robinson-Schensted correspondence, to the Burge correspondence, and to the version of the dual RSK 
correspondence using column insertion described in our proposition 3.3.3. For the Robinson-Schensted 
correspondence one takes X/k and u/n to be shapes that are both horizontal and vertical strips, so that 
Pic(A/K, v/ii) consists of all bijections [X/k] — > [v/fi\, and that Pic(7r, v/n) and Pic(A/«, n) are both in 
bijection with the set of standard Young tableaux of shape n. For the other two cases, one takes v / 'fi to 
be a horizontal strip, and X/k to be a horizontal respectively a vertical strip. Then the operations Int 
and Bin define bijections Pic(A/«, v / \i) — ► M.\- K , v -p respectively Pic(A/K, v / \i) — > M. l * ] _ \t_ K t, and the 
sets Pic(7r, v/n) and Pic(A//t, n) are in bijection with sets of semistandard Young tableaux of shape tt 
as described above. Thus these latter two special cases of the decomposition of pictures give rise, under 
application of Bin and Int, precisely to the bijections that occur in our theorem 3.1.3. 

The description of the decomposition of pictures using the column insertion Schensted algorithm 
is rather different from the definition of our bijections in theorem 3.1.3. However, the decomposition of 
pictures has an alternative description that is given in [vLee2, theorem 5.4.1]. This description is in terms 
of two forms of jeu de taquin that are defined for pictures, one that changes the shape X/k of the domain, 
and another that changes the shape u/fi of the image of the picture. The two types of operations commute 
[vLee2, theorem 5.3.1], and a decomposition of pictures based on these operations follows in a manner quite 
similar to the way we obtained theorem 3.1.3. From the fact that for pictures Schensted column insertion 
can be emulated by jeu de taquin (as it can for tableaux), it follows that the two descriptions define the 
same decomposition of pictures. Under the operation Int, each domain slide of a picture corresponds 
to a sequence of upward crystal operations on integer matrices, just like in proposition 2.5, and each 
image slide corresponds to a sequence of leftward crystal operations; hence the decomposition of pictures 
translates directly into the decomposition of integral matrices of theorem 3.1.3(2). For the operation Bin, 
the situation is similar, except that domain slides correspond to leftward crystal operations and image 
slides correspond to upward crystal operations. Thus the decomposition of pictures also translates into 
the decomposition of binary matrices of theorem 3.1.3(2). 

With these links between the decomposition of matrices and that of pictures, and the fact that the 
latter can be computed using Schensted column insertion, the descriptions of the previous subsection 
that involve column insertion are easy to understand. For instance, the facts that in proposition 3.3.2 
the integral matrix M is traversed in Semitic order while in proposition 3.3.3 the binary matrix M is 
traversed in Kanji ordering both reflect the fact that the column insertions used to compute the first 
factor of the decomposition of a picture traverse its domain from top-right to bottom-left (in fact any 
order will do that never makes a decreasing step for '</>'). This is so because M a .j3 is modelled by 
Pic(A/K, v/\i) where X/k and v/[i are horizontal strips with A — k = a and v — fi = (3, and with the 
rows of the matrix corresponding to those of X/k, while p is modelled by Pic(A/«;, v/fi) where X/k 
is a vertical strip and u/fj, a horizontal strip, with A* — k* = (3 and v — \i = a, and with the columns 
of the matrix corresponding to those of X/k. In both cases top-right to bottom-left traversal of [X/k] 
corresponds to the indicated traversal of the matrix. The description of the recording tableaux is also a 
direct translation of the way recording tableaux are defined for pictures. For propositions 3.3.1 and 3.3.4, 
which use row insertion, the more complicated formulations reflect a less direct correspondence to the 
decomposition of pictures; they can however be understood by using the symmetry properties of pictures 
that are discussed in [vLee2, §4]; we shall not elaborate the details here. 
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§4. Growth diagrams. 

After what was said in the last subsection, one might be tempted to conclude that crystal operations on 
matrices and the correspondences of theorem 3.1.3 are merely shadows of similar structures for pictures 
under the maps 'Int' and 'Bin', and that as such they do not add much insight, although they are somewhat 
simpler to present than pictures and the operations of jeu de taquin on them. In this section we shall 
however establish a link with the descriptions as those in [Roby] and [Fom] of Knuth correspondences 
using growth diagrams, which link does not exist in the theory of pictures. In doing so we build a bridge 
between the previously unrelated approaches to Knuth correspondences via jeu de taquin on one hand, 
and via growth diagrams on the other hand (the approach via jeu de taquin often also uses invariants, 
an in [Schii2] or [Fult]; we shall discuss these below in §5). 

4.1. Submatrices and implicit shapes. 

For a discussion of the notion of growth diagrams we refer to [vLee4, §§2,3]. In general, a growth diagram 
associates "shapes" to the elements of the grid N x N, and it also associates "entries" to squares of the 
grid; in our case the shapes will be partitions, and the entries are our usual matrix entries, lying either 
in [2] = {0, 1} or in N. Formally a square of the grid is an element of N x N, just like a grid point is, but 
the "square" (i, j) is supposed to have as corners the "grid points" in {i, i + 1} x {j, j + 1}. We shall build 
a growth diagram from a given matrix M by associating to the grid point (k, I) a shape determined by 
the submatrix (^i,j)ie[fe],jem 5 of which that grid point is the lower right hand corner. The shape will in 
fact be the one that parametrises the normal form of that submatrix, which by our usual interpretation 
of finite matrices is identified with the matrix obtained from M by clearing the entries outside [k] x [I] ; 
we shall call this the implicit shape of the submatrix. The formal definitions are as follows. 

4.1.1. Definitions. Let M e Ai l2] or M E M, and let N be the normal form of M, obtainable from it 
by exhausting the operations chosen from { e\ | i € N } U { ej~ | j £ N }. We dehne 

a. 7r(M) = row(iV), which we call the implicit shape of M; 

b. Mi x j = ([ E I x J] Mij)i.j e N for intervals I, J C N, which is the matrix of the same type 
(binary or integral) as M and obtained from if by replacing all entries M^j with I x J by 0. 

The term implicit shape refers to the fact that it cannot be easily read off directly from M, although 
in §5 we shall give a characterisation that does not require modifying M. Once either of the factors P, Q 
associated to M under the decomposition of theorem 3.1.3 has been determined, the implicit shape of M 
can however be read off, for instance as ir(M) = row(P); in other words it is the shape of the tableaux 
corresponding to M under the Burgc correspondence or the column insertion dual RSK correspondence. 
If M is the encoding of a skew semistandard tableau T, then n(M) is the shape of the rectification of T. 

Growth diagrams should satisfy a local rule, which we called "shape datum" in [vLee4], that relates 
the shapes at the four corners of any grid square and the matrix entry in that square. This rule is 
such that when the shapes at (i + and + 1) arc arbitrarily fixed, it defines a bijection between 
the possible values for the pair consisting of the shape at (i, j) and the matrix entry on one hand, and the 
possible values for the shape at + + 1) on the other hand. When the shape datum is prescribed, this 
permits deducing the entire growth diagram from a partial specification of it, namely from the shapes 
along any lattice path in the grid from bottom left to top right, together with all matrix entries to the 
bottom right of that path. In our current situation this will eventually provide alternative methods to 
compute implicit shapes, but at this point we just have the matrix of shapes (tt(M^ x \j]))%,j^N, which is 
defined without any reference to a shape datum. So we must first prove that it indeed defines a growth 
diagram, and determine the shape datum for which it does so. 



4.2. Growth diagrams for integral matrices. 

The key to proving that (tt(M^ x [j] ))i,jeN defines a growth diagram will be the study of a particular way of 
reducing M to its normal form. Consider for the moment integral matrices, and take the matrix M of (29). 
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4-2 Growth diagrams for integral matrices 

One sees that the submatrix indexed by [3] x [6] has been reduced to the normal form Diag((8, 4, 2)). 
We shall first argue that this is indeed the normal form of the original submatrix M[ 3 ] x [ 6 ] , and that this 
reduction is no accident (a similar reduction can be achieved for any integral matrix). 

4.2.1. Proposition. Fix M G M, and k,l G N>o- Let M' G M be the matrix obtained from M by 
exhausting the operations in { ej \ i G [k — 1] } U { ej~ | j G [I — 1] }, and let N G M be the normal form 
of M[ fe ] x [;] , which is obtained from it by exhausting the same set of crystal operations. Then x ^ = N. 

Proof. It is clear that the normalisation of M[ fe ] x ^ cannot involve any raising operations other than the 
kind being applied to M, as M[ fe ] x [z] contains no units that could be transferred by such raising operations. 
It is also clear that the operations applied to M cannot transfer units across the boundaries of the rectangle 
[k] x [I]: they cither transfer units within that rectangle or within its complement. If we extract, from a 
sequence that exhausts these operations applied to M, those that involve a transfer within [k] x [I], then the 
extracted sequence can be applied to M\u\ x m , and each operation will involve the same transfer as it did 
in the original sequence; this is so because any valid transfer in a matrix remains valid in any submatrix 
containing the squares involved in the transfer. Therefore the extracted sequence reduces M[j,] x [z] to 
Mr' fe j x j Z j, and it remains to show that the latter matrix is a normal from, which then necessarily equals N. 

This means that we have to prove that n](M^ x ^) — nj~ {M^ x ^ ) = for any i G [k — 1] and j G [I — 1], 
while we know that n](M') = ny(M') = 0. Now if U\V denotes an integral matrix vertically split into 
submatrices U,V, then it follows from (18) that n\{U\V) > n\(U), so n\{U\V) = implies n\(U) = 0. 
For a horizontal division into submatrices, (20) similarly gives nf (|) > rij-(S), so that n^(|) = 
implies ny(S) = 0. These relations allow us to conclude n\ (M^ x ^) = n j~(M^ x ^) — as desired. □ 

Note that while the final part of this argument amounts to showing that if any raising operation 
can be applied to a submatrix of the form M^] x [i], then the same operation could also be applied to 
the full matrix M, this does not mean that these operations would necessarily involve the same transfer. 
Indeed, if one has n\{U\V) > n\(U) > 0, then the first application(s) of e\ to U\V will be between entries 
of V, and only once n](U\V) and n]{U) have become equal will the remaining applications of e\ transfer 
entries within U . The same is true for leftward moves with respect to a horizontal split. In the example 
given one sees that the exhaustion of operations in ej^ involves leftward transfers below as well as above 
the horizontal line; the former precede the latter in case of transfers between the same pair of columns. 

Let us now consider the relation between the four shapes ^(M^xyj) for (i, j) G {k, k+ 1} x {1,1+ 1} 
and the matrix entry M k .i- All of these values only depend in the submatrix -Mj fe+1 ] x [z+i], and none of 
them changes when one applies operations in { e\ \ i G [k — 1] } U { ej~ \ j G [I — 1] }. Therefore let A be 
the matrix obtained from M[ fe+1 ] x [ ;+1 ] by exhausting these operations, and to simplify the notation put 
A = 7r(M[ fe ] xW ), /j, = !r(M Wx[l+1 ]), v = 7r(M [fe+1]xW ), and n = 7r(M [fe+1 ] x [j + i]), which arc also the 
implicit shapes of the corresponding submatrices of A (here we abandon our previous conventions for the 
use of A, ... ,k, and adopt those of [vLee4]). By proposition 4.2.1 we have ^4[fe]x[;] = Diag(A). 

The normalisation of -A[fc+i] x [z] will only involve upward transfers, since the leftward operations in 
{ ej~ | j G [I — 1] } remain exhausted throughout. This normalisation transfers any units present in (the 
bottom) row k of ^4.[fe+i] x [i] upwards, until they arrive on the main diagonal, where together with the 
units that were already there in ^4[fe] x [;] = Diag(A) they form the matrix Diag(^). It follows that row k 
of j4[fc + i]x[z] must have been equal to the composition v — A (and in particular that it contained no units 
strictly above the main diagonal). A similar argument shows that the final column I of j4[ fc ] x [;+i] is equal 
to (j, — A. The matrix A is therefore completely determined by the values of A, /i, v, and Mk y V- one has 
A [k]x[i] = Diag(A), Ai.i = in - A; for i G [k], A kJ = v. } - Xj for j G [I], and A k j = M k j. This establishes 
one part of requirement for families of shapes (7r(M[ i ] x [j]))i,jeN to define a growth diagram, namely that 
A, /i, v, and M kt i together determine k; it also provides a procedure to compute this direction of the 
shape datum (namely construct the matrix A and find its implicit shape). It is easy to see that one has 
M = H + \ u \ ~ 1^1 + as required by the definition of a shape datum [vLee4, definition 2.2.1]. 

Let us apply the procedure for computing k in the case of the matrix M of (29) that was used in 
the example above, for k = 3 and / = 6; so our goal is to compute k — 7r(M[ 4 ] x [7]). From the last matrix 
in (36) one reads off A = (8, 4, 2), fj, - A = (0, 3, 0) so that fi = (8, 7, 2), v - A = (0, 1, 1, 1, 0, 0) so that 
v = (8, 5, 3, 1), and = 2. We may now determine the following diagram, from which it can be read 
off that k — (8, 8, 4, 2) (unlike the previously drawn diagrams we now let the arrows between the matrices 
point in the direction opposite to the crystal operations that are being applied; this is done so that the 
rectangle [i] x [j] for which the submatrix has been normalised to 7r(A^ x yj) grows from left to right and 
from top to bottom). 
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Diag(fc). 



To complete our argument we shall also have to show that conversely /i, v, and k determine A 
and Mk t i, but before we do so we must specify which relations should hold among the shapes for the 
above construction to apply. First of all, each implicit shape tt(M^ x ^) can have at most mm(i,j) nonzero 
parts. Apart from that and A C fi and ACv (which must hold since neither v — A nor fi — A have negative 
entries), there is only one requirement for A, /i, and v, namely that the matrix A constructed from the 
pieces Diag(A), /j, — A, \i — A and Mjy as indicated must satisfy the hypothesis of initial exhaustion, in 



other words that n\(A) 



'(A) = for all i G [k — 1] and j G [/ — !]• This amounts to the requirement 



Aj+i — Aj + (/i — A)i+i < or equivalcntly /ii+i < A^ for i G [fc — 1], and similarly Vj+\ < Aj for j £ [Z — 1]. 
Together with the limit on the number of nonzero parts, this implies that A, /i, and v satisfy the relations 



A n and A v- v. Conversely, for any \,fj,,f G V satisfying /i m in(fe,j+i) = ^minffc+i.n = 0, and the 
relations A v- fi, and A v- v, the construction of A (with an arbitrary entry Ak,i G N) provides an 
example where A, /i, and v occur as the implicit shapes of appropriate submatrices. 



Now assume that /j,,u,k € "P are given with fj, 



min(fe,;+l) — u min(k+l,l) 



= 0, 



and We 



want to find ^4 as above, giving rise to these shapes. We know its normal form Diag(fc), as well as the 
normal forms Diag(((i) and Diag(^) of -<4.[fc] x [i+i] ano ^ ^[fe+i]x[;]; respectively. Exhausting the operations 
{ ej~ | j G [I] } applied to A must give the (k + 1) x (I + 1) matrix B whose restriction to [k] x [I + 1] 
equals Diag(/i) and whose row k equals k — [i. The normal form Diag(^c) of B can be written as R 1 '(B) 
for some sequence i? T of upward crystal operations in { e\ \ i G [k] }. Then i? T can also be applied to A, 
exhausting upward transfers and resulting in the (k + l)x(l + l) matrix C whose restriction to [k + 1] x [I] 
equals Diag(/i) and whose column I equals n — v. Then inverting the operations in i? T , it is clear that 
A = R l (C), and from this matrix one can extract the values A = 7r(^4[ fe ] x ^]) and Mk,i = Akj- 

As an example of the computation in the inverse direction, let partitions k = (13,9,9,5), /i = 
(11,9,8), and v = (10,9,8,2) be given, and take k — 3 and I = 4 (the minimal possible values, 
given jU and v). We form a commutative diagram like the one drawn above, but starting with the 
matrices B, C, and Diag(«) below, which correspond to the initial data. One easily determines com- 
posite raising operations transforming B or C into Diag(fc), for instance i? T = (ej) 2 (e\) 2 (e^) 3 and 
R*~ = (e^ ) 3 (e^) 3 (e^) 4 (e^) 7 (the exponents are just initial partial sums of the parts of the compositions 

- v = (3, 0, 1,3)). Then one computes A by inverting one of these composite op- 
R l (C) where R l = (fe) 3 (f±) 2 (/ 1 ) 2 , and from A one then reads off A = (9, 9, 6) 



k — [i= (2, 0, 1, 5) and k - 
erations, for instance A = 
and M 3j4 = 3. 
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Diag(«). 



As is clear from the examples, the only aspect of the reduction A — > Diag(fv) that is not completely 
straightforward is the transformation during the first reduction phase of the row or column that is 
outside the submatrix being normalised in that phase. Consider for instance applying upward transfers 
to reduce the submatrix A[ fe+1 ] x [;j to its normal form Diag(i^). Here column I, which initially describes 
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the composition a = ji — A and the entry Mk t i, is transformed into a composition a' that will determine 
k = v + a'. The operations e\ can be applied by weakly decreasing indices i. For fixed i, the effect on 
rows i, i + 1 of exhausting the operation e\ is a transformation of the following form (with (i = v — A): 

/o •■■ o a, o o ••• o «A (eir /A) ••• A-i ^ ••• m' \ 

^/3 ••• A-i fti/itiO-Omj V ... ■■■ a{ + J' 1 j 

where = min(m, Aj — fj+i), and m' = on + m — a' i+1 . The number m is determined by the previous 
raising operations, except for the initial case i = k — 1, where m = M^j. Since the diagonal entry i/j 
is unchanged by further applications of e\_ 11 and Vi > Aj > fj+i + aij +1 , there is no need to reapply 
later. Thus we have established a simple iterative procedure for determining the transformation of the 
final column of the matrix, which will then determine n. 

We rewrite this procedure into an algorithm for directly computing k from A, ji, v, and M^j, by 
taking into account a = fi — A and k = v + a'; in this procure the roles of \x and of v will be symmetrical. 
Adding v i+i to the equation for a' i+1 we obtain m+i = min(m + fj+i, Aj). The amount in' — on added 
to on is equal torn- ct' i+1 = m + — Kj+i; we shall use this "carry" c in our algorithm rather than m. 
To that end we replace m in our description by the sum of the original entry ctj+i = Hi+t — Aj+i and 
the value c left by the previous step; initially we take that expression with c = Mjy. Then we get 
= min(d, Aj) where d = /Uj+i — Aj+i + c + I'j+i, and d — will replace c. The index z + 1 now 
occurs much more often than i, and we shall therefore shift the index one up, so that i + 1 becomes i. 
Finally, the last part of n that can be nonzero is Kmi^fc.f), so we can start at i = min(fc, /) instead of at 

1 = k and gain in efficiency and symmetry. After these transformations we get the following algorithm: 

c := M M ; 

for i from min(fc, I) down to 1 

do d := Hi + Vi — Aj + c ; Kj := min(d, Aj_i) ; c := d — Ki od ; 
k := Mo - A + c + z/ . 

If we execute the algorithm for the values of the example above A = (8, 4, 2), /i = (8, 7, 2), v = (8, 5, 3, 1), 
k = 3, I = 6, and = 2 then we find for i = 2: <i := 3, K3 = 2, c := 1; for i = 1: := 4, «2 = 4, 
c := 0; for i = 0: d := 8, Ki = 8, c := 0, and finally k = 8, so that indeed one finds k — (8, 8, 4, 2). 

Note that the algorithm ensures that m+i < Aj for all i, which means that A k will always hold; 
this is a stronger condition than follows just from the known relations A t— fi k and A v k. 
In fact the correspondence between the four partitions and the entry can be characterised by the 
following properties: (1) there exist numbers Cj G N for < i < min(fc, I) such that putting c_i = one 
has Cj + Hi + Vi = Aj + Kj + Cj_i for all those i; (2) «j+i < Aj for all i, and Cj > implies «j+i = Aj; 
(3) c m ; n (fc.() = Mfe^. For i > 0, a will be the value of c in the algorithm at the start of the loop indexed 
by i, while Co is the final value of c. One may eliminate the dependence of this statement on (k,l) by 
having i range over N in property (1) and replacing property (3) by limj^oo Cj = M/~j (in this case 
the sequence Cj becomes stationary for i > min(fc,/)). This characterisation is useful for deducing the 
following algorithm inverse to the one above, which will compute A and Mjy from /x, v, and k. 

c := ; 

for i from to min(fc, I) 

do d := Hi + v % — K « — c ; Aj := max(rf, Kj+i) ; c 
M M := c. 

Executing this with as above k = (13,9,9,5), fi = (H>9>8), ^ = (10,9,8,2), fc = 3 and Z = 4 gives for 

2 = 0: d := 8, A = 9, c := 1; for i = 1: d := 8, Ai = 9, c := 1; for i = 2: d := 6, A 2 = 6, c := 0; for i = 3: 
d := —3, A 3 = 0, c := 3, and finally M^a — 3, which together with A = (9, 9, 6) is the result found before. 

It should come as no surprise that the shape datum we have just described is in fact the shape datum 
for the Burge correspondence, which was given in [vLee4, §3.2]; after all we know from proposition 3.3.2 
that the implicit shape computed for M G M by the growth diagram for our shape datum is the shape of 
the tableaux associated to M under the Burge correspondence. Although none of the descriptions given 
in that paper coincide exactly with the one given in (40), it is not hard to see their equivalence. 

Those descriptions compute [n] by making a bottom-left to top- right pass over the squares that could 
be added to [A] (namely the topmost squares of the columns of the complement of that diagram), using 
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a "stock of squares" whose size may grow and shrink. Initially the stock contains M^j squares, and at 
each square s visited, a square is added to the stock if s £ [a*/ A] and if s £ [v/X] (if both conditions hold, 
two squares are added), and then a square from the stock is contributed to [k] unless the stock is empty 
(so if s lies either in \p/X] or [v/X] it will certainly be in [k], and if it lies in both there will in addition 
be a net increase of one square in the stock; if it lies in neither and the stock was not empty then s will 
still be in [n] while the size of the stock is decremented). Algorithm (40) can be seen to compute this 
same diagram [k] a whole row at the time. The value d = jJLi + v% — A« + c = Aj + c + (p — A)j + (y — A)j 
represents the maximal length row i of [«] can have, obtained by adding to Aj the size of the stock when 
arriving at this row and the sizes of rows i of [///A] and of [v/X]. Either that length does not exceed the 
limit Aj_i, in which case the one sets Ki = d and the stock becomes empty, or d does exceed Aj_i, in 
which case Ki is set to that limit and the size of the stack becomes d — Aj_i = d — Ki. 

It is interesting to see how this description of the shape datum for the Burge correspondence arises 
naturally from the rules governing the normalisation process of integral matrices. Summarising, we have: 

4.2.2. Theorem. If for any M £ M. one associates to each grid point (i,j) £ N 2 the shape tt(M^ x ^), 
then one obtains a growth diagram for the shape datum of the Burge correspondence. □ 

Having determined the shape datum involved, we can now determine the family of shapes n(M^ x yj) 
directly from the matrix M, without actually performing any crystal operations: it suffices to use the fact 
that Tr(Mu] x ij-i) = (0) when i = or j = 0, and to construct the growth diagram for the shape datum. 
Diagram 1 illustrates this for the integral matrix M of our running example. 
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Diagram 1. Growth diagram with integral matrix and shapes {n(M^ x ^))ij e ^i. 

From this growth diagram one not only can read off n(M) in the bottom right corner, but also the 
matrices P, Q associated to M under the bijection of 3.1.3(2). To see that, observe that the operation 
M i ► P of exhausting upward transfers commutes with the operation M i— ► M Nx rji of discarding all 
columns beyond a certain point; this is just the fact observed above that n\(U\V) = implies n\ (U) = 0. 
Therefore the matrix fWx[j] can be obtained from M Nx [j] by exhausting upward transfers, and one has 
row(Pisr x [j]) = 7r(-PNx[j]) = 7t(Mn x [j]). Putting X^ = row(pNx[j]) for all j defines the semistandard 
Young tableau S = (A^)j £ N of which P is the integral encoding, and we see that S coincides with the 
sequence (7r(M Nx r J i))j e N that can be read off along the lower border of the growth diagram. Similarly 
Q l is the integral encoding of the semistandard Young tableau L = (n(M^ X N))ieN that can be read 
off along the right hand border. Indeed in the example of growth diagram 1, one finds S of (31) and L 
of (33), whose integral encodings are the matrices P and Q t of (29). 

This argument confirms (and provides an alternative proof of) proposition 3.3.2, which can be 
interpreted as giving one particular way in which the growth diagram can be computed: column insertion 
from right to left of the letters of the weakly increasing word w whose weight is given by the row Mj gives 
the transition from the semistandard Young tableau (n(M^ x yj))j e ^, formed by the shapes along row i 
of the grid, to the tableau (7r(M[ i+1 ] x [ J -]))j eN along the next row of the grid. Indeed one may check that 
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the successive stages of the insertion process in (32) correspond to the shapes along the successive rows 
of the growth diagram above. (It may be a bit confusing that whereas the letters of w are inserted from 
right to left, the growth diagram is extended from left to right; however, completing the square (i,j) in 
the growth diagram does not correspond to the insertion of the letters j of w, but rather to determining 
the effect on entries < j of the insertion of w, relative to the effect of that insertion on entries < j.) 

4.2.3. Summary. Given M E A4, the shapes tt(M^ x [^) can be determined as follows for all i,j E N 
(as there are only finitely many distinct matrices Mui x ui, the computation is effectively finite). Starting 
with 7r(JUf[i] x ij-] ) = (0) whenever i = or j = 0, one completes a growth diagram by repeatedly computing 
shapes k — 7r(M[ i+1 ] x [ J+1 ]) from previously computed shapes A = tt(M^ x ^), fj, = 7r(M[ fe ] x y +1 ]), and 
v = 7r(M[fc + i] x [j]), and the entry Mk,i, using the algorithm (40). From this growth diagram one can 
read off the semistandard Young tableaux S = (ir(Misix[j]))j€'isi &nd L = (7r(M[ i ] X N))i £ N associated 
to M under the Burge correspondence, as well as the matrices P, Q associated to M under the bijection 
of theorem 3.1.3(2): P is the integral encoding of S and Q f is the integral encoding of L. One can 
reconstruct M from (P, Q) or from (S, L) by reversing the computation, using algorithm (41 ) to compute 
A and Mk t i given \i, v, and k. 

4.3. Growth diagrams for binary matrices. 

We shall now consider the family of shapes n(M^ x y^) for binary matrices. With respect to the integral 
case there is a complication, because one cannot in general reduce a submatrix M^ x ^ to its normal 
form using raising operations acting on the entire matrix. For instance after exhausting e$, C\ , and ej 
for the binary matrix M of (28), its submatrix M[ 2 ] x [3] is reduced to J °) which is not a normal form; 
in isolation the bottom right entry could be moved up by ej, but in the whole matrix it is blocked by the 
parts of rows and 1 further to the right. 

The difficulty comes from the reversal of left and right in the rule n\(U\V) > n\(U) that we used 
for integral matrices: for binary matrices one instead has n\(U\V) > n\(V) (for horizontal crystal 
operations there is no such difference with the integral case: the relation n^(^) > rij~(S) still holds 
for binary matrices). So we cannot exhaust upward moves in the submatrix U; however we have the 
rule n\(U\V) > n\(U) that allows us to eventually exhaust any downward crystal operation // on U 
(after a finite number of applications of // that move bits in V^). One cannot exhaust all downward 
moves on a matrix, but if fc > is such that M = M[ fc ] xN , then it will suffice to exhaust the operations 
// for i E [fc — 1]. If moreover leftward moves are also exhausted, then the matrix obtained from M 
by a quarter turn clockwise, mapping the topmost fc rows to the leftmost fc columns, will satisfy the 
conditions of lemma 3.1.2; it follows that M is obtained by the inverse rotation from Diag [21 (A'), where 
A* = col(M), and in particular that M is determined by its implicit shape A together with fc. It is in fact 
an upside-down version of Diag [2] (A), in other words the region occupied by its bits '1' is the French style 
display* of the diagram of A, in which the parts of A correspond to the rows taken from bottom to top. 

4.3.1. Definition. Let A E V and fc E N with A fe = 0. We denote by NF(A, fc) the binary matrix 
([ (fc — 1 — i, j) E [A] ] )ijgN, which is called the French normal form (relative to k) parametrised by A. 

Note that col(NF(A, fc)) = A', while row(NF(A, fc)) = (Afc_i, . . . , Ao); we shall write this composition 
obtained by reversing the fc initial parts of A as revfe(A). Now by exhausting the crystal operations in 
{ // | i E [fc — 1] } U { ej~ | j E [I — 1] }, the submatrix M^ x ^ of any binary matrix can be reduced to the 
form NF(A, fc), and this will be for A = 7r(M[ fe ] x [;]). We call this process French normalisation of the 
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For instance for the binary matrix M of (28) and (fc,^) = (6,4), this exhaustion of 



crystal operations gives: 
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One might have imagined that the occurrence of Young diagrams inside matrices, as in lemma 3.1.2, would 
give a coup de grace to the French style of displaying diagrams. Mais non. 
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We see the submatrix Mrgi x r4] has been transformed into NF((4, 4, 2, 1), 6), whence we may conclude 
that 7t(M[ 6 ] x [ 4 ]) = (4,4,2,1). As in the integral case, we shall now compare the results of the French 
normalisation of submatrices M[j] x for (i,j) G {k, k + 1} x {I, l + l}. Let A be obtained from M[ fe+1 ] x [/ +1 ] 
by exhausting the operations in { // | i G [k — 1] }U{eJ~ | j G [I — 1] }, and as before put A = n(M^ x ^), 
/i = n(M [k]x[l+1] ), v = 7r(M [fe+1]x[i] ), and k = n(M [k+1]x[l+1] ). Then A [k]x[l] = NF(A, k), and since 
French normalisation of A[fe] x [;+i] an d of -4[fc+i] x [Z] oru y involves horizontal respectively vertical moves, 
and it transforms these submatrices into NF(/i, k) respectively into NF(^, k + 1), column I of A[fej x [; +1 ] 
is equal to the composition revk(^ — A), while row k of A[ fe+1 ] x [;] is equal to v l — A*. In particular one 
always has A •<— /j, and A v- v. In the example displayed above, where A = (4, 4, 2, 1) and the 7x5 
matrix A can be read off from the matrix on the right, one finds fj, = A + (1, 1, 0, 0, 1, 0) = (5, 5, 2, 1, 1) 
and v = (A* + (0, 1, 1, 1))* = (4, 4, 4, 2). 

As in the case of integral matrices we see that A, and therefore k = n(A), is determined by the shapes 
A, /i, v, and the bit M k ^\ moreover for any such data the construction provides an example M = A of 
a matrix with the specified implicit shapes of submatrices and bit M k ,i, provided that one has A it, 
A t— v, and [i k = v\ = 0. By the same procedure as for integral matrices one can also reconstruct A 
and Mkj from /i, v 1 and k. Thus for any /i,i/£? with n k = v\ = one has a bijective correspondence 
between (A, M k ,i) on one hand and k on the other, where these values are restricted by the conditions 
A /j,, A t— v, Mkj G [2], v 1 — k and /j, k. Following the terminology of [vLee4, definition 6.1], these 
bijective correspondences for all fi, v G V together define an (asymmetric) shape datum for [V, t— , [2]). 
By assigning n(M^ x yj) to each grid point one obtains a growth diagram for this asymmetric shape 

datum. Note that the shapes associated to grid points that are horizontal neighbours differ by a vertical 
strip, while those associated to vertical neighbours differ by a horizontal strip (we have no choice); in 
[vLee4] an opposite convention was used (its shape datum was for (V, t— , , [2]), not for (V, t— , [2])). 

Before considering this shape datum more closely, let us compute an example. The example of (42) 
for (k,l) = (6,4) shows that during the French normalisation of A^+i^^ all bits '1' either remain in 
place or move one place down, while during the French normalisation of A[j.] x [; +1 ] some bits '1' may move 
several places to the left (the bit at (1, 4) moves all the way to (1, 0)); however this example is too simple 
to illustrate well what can happen in general. Therefore consider instead A = (5, 3, 2), /i = (6, 3, 3, 1, 1), 
v = (6, 3, 2, 2), and Mfc ; ; = 1 with (fc, I) = (5, 6). For the corresponding matrix A one has: 

/0 

. _ 1 1 
1110 
111110 
Vl 1 1 

(f( 5] )\ 

/0 



1 1 
1 1 
1110 

Vi i i i i i 

for k = (7, 4, 3, 2, 1, 1). The parts of this computation that are not easy to predict are the following ones. 
Firstly there is, for the exhaustion of downward moves applied to A, the transformation of column I of A 
(which represents vevk{n ~ ^) followed by M ktl ) into rev fc+1 (K — u); here some bits T move down a row 
(like in the example M 2 ,6; this corresponds to a square of [/x/A] for which [n/v] contains a corresponding 
square in the same row), while some bits '1' remain in the same row (in the example bits Mo. 6, M\^, M^, 
and M5,6; this corresponds to a square of [/x/A] (or Mk,i = 1) for which [k/v] contains a corresponding 
square in the next row (respectively in row 0)). Secondly there is, for the exhaustion of leftward moves 
applied to A, the transformation of row k of A (which represents v t — A* followed by Mk,i) into k* — /i*; 
here the bits '1' may exchange with zero or more bits '0' initially situated directly to their left (in the 
example bit M 5j5 interchanges successively with two bits '0' to its left while M 5 , , M 5i i, and M 5i6 stay 
in place; such a bit corresponds to a square of [f/X] (or Mk,i = 1), and the corresponding square of [n/(i\ 
lies in a column weakly to the left of it (respectively in a column < I)). 
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Giving details of cither of these transformations will give a description of how to compute k. However, 
although these descriptions follow by a straightforward application of the definitions, they are rather 
technical, and both are not stated in the same terms (the former is primarily by rows of the shapes 
involved, the latter primarily by columns), which makes it hard to see that they give the same result 
for k (as we know they must). Therefore we shall instead describe this shape datum using the approach 
that was detailed in [vLee4, §6]; we have to adapt that description to the interchange of '-5—' and lj — ' in 
the relations between shapes at horizontal and vertical neighbours in the grid. 

Given fi, v £ V for which there exists at least one A £ V with fi —r A v, or equivalently at least 
one it £ P with n k —r v, the possibilities for such A and n have a rather simple description: in 
each case the diagrams of these partitions are the disjoint union of a fixed set of "obligatory" squares 
and an arbitrary subset of a fixed set of "optional" squares. In case of A, the Young diagram [u] n [v] 
is the union of all obligatory and optional squares (and therefore the maximal choice for [A]), and its 
subset S of optional squares consists of those that are both at the end of a row of \p] and at the end 
of a column of [v], in equation: S = { £ N 2 | (i, j) — (vj — l,m — 1) }. In the case of k, the set of 
obligatory squares is the Young diagram [fj] U [v] , and the set T of optional squares consists of those that 
are both just beyond the end of a column of [fj] and just beyond the end of a row of [v\, in other words 
T = { ihj) £ N 2 I {hi) = (/"j) y i) }■ It was shown in [vLee4] that #T = #5 + 1, which is in agreement 
with the fact that there are twice as many possibilities for k as for A; moreover, there is an injective 
map matching every s € S with the first t £ T in a row strictly below s (or equivalently in a column 
weakly to the left of s), which leaves the topmost square t of T unmatched. Now the correspondence 
between (A, M fcj ;) and k is determined by the fact that for every optional square s £ S absent from [A] 
the corresponding optional square t £ T is present in [k], and vice versa, while one has to £ M if an d 
only if Mk : i = 1. We shall call this shape datum the "row-insertion shape datum for (V, [2])". 

4.3.2. Theorem. If for any M 6 M [2] one associates to each grid point (k , I) £ N 2 the shape tt(M^ x ^), 
then one obtains a growth diagram for the row-insertion shape datum for (V, t— , [2]). 

We omit the proof; the assiduous reader may check that the rules governing French normalisation of 
submatrices do indeed lead to this shape datum. Whether one starts normalising A^+i]x[i] or ^4[fc]xR+i]i 
the description obtained will not directly give a correspondence between S and T, but rather one between 
the sets of rows meeting [fi/X] and [n/v], respectively between the sets of columns meeting [v/\] and 
On one hand this means that a large number of steps may intervene between encountering the bit 
for some s £ S and the one for the corresponding t £ T, and on the other hand the occurrence of such 
squares may be hidden by the fact that a square in S shares a row or a column with a different square 
in T. As a consequence, cases need to be distinguished that play no role in the above description. 

Now that the shape datum is known, it can be used to build up the growth diagram without actually 
performing (French) normalisation submatrices; again one starts by associating the empty partition (0) 
to every grid point with i = or j = 0. The construction for the binary matrix M of (28) is given 
in diagram 2, with the shapes associated to grid points given by their diagram a la frangaise. 

The shapes along the right border define the semistandard tableau S of (31) whose binary encoding 
is the matrix Q of (28), for essentially the same reasons as in the integral case: one has col(<5m X N) = 
7r(Q[j] xN )' = n(M[q xN y for all i £ N. There is however no such direct relationship between matrix P 
of (28) and the shapes along the lower border; this is not surprising since those shapes define a transpose 
semistandard tableau, whereas P corresponds to the reverse transpose semistandard tableau R of (35). 
Indeed, although like in the integral case one still has 7r(P Nx [j]) = 7r(M Nx y]) for all j, these shapes are 
not given by row(P Nx [j])> since upward moves need not be exhausted in Pnx[j]- To read off R from a 
growth diagram, one should instead normalise the complementary submatrices MNx(N-[j])> as we shall 
do in the next subsection. On the other hand, our growth diagram can be interpreted as recording the 
insertion process of proposition 3.3.4, which gives as insertion tableau the Schutzenberger dual R* of R: 

| 1 | 2 | 3 | 4 | 5 | 6 | 7 
1 2 3 4 6 1 7| 8 

1 2 3 1 4 1 

1 2 1 5 | 
jj 

Indeed, since the shapes in our diagram are governed by the row-insertion shape datum for (V, ^— , [2]), 
the passage from the transpose semistandard Young tableau defined by the shapes along row i to those 
along row i + 1 is given by the row insertion (according to the rules for transpose tableaux) of the letters 
of the strictly increasing word of weight M i} from left to right. This proves: 
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Diagram 2. Growth diagram with binary matrix and shapes (n(M^ x ^))i t j^ . 

4.3.3. Proposition. For M e M m , the transpose semistandarcl tableau R* — (7r(M Nx [ J ]))jgN and the 
semistandard tableau S = (7r(M[ i ] xN )) ie N form the pair associated to M in [Stan, theorem 7.14.2]. □ 

From M one can obtain a matrix P whose rows give the weights of those of R* , but in reverse order, 
by exhausting downward crystal operations, up to some limit. Exhausting { // | i e [6] } in our example 
gives 

/0 0000000 0\ 
000000000 
000010000 
1 1 1 
1 1 
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1 1 



111110 111 

Vi 1111111 o/ 



which matches 



R* 
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(44) 



4.4. Other types of growth diagrams. 

As was suggested above, one can form a growth diagram from which both tableaux R and S encoded by 
the matrices P and Q corresponding to M G Ai l2] under the bijection of theorem 3.1.3(1) can be directly 
read off, by computing the family of shapes (7r(M[ i ] x ( N _y])))j i j e N rather than (7r(M[j] x [j])), ; j e N. The 
shapes defining R, a weakly decreasing sequence in which adjacent shapes differ by vertical strips, can 
then be found along the bottom border, while the shapes defining S are now found along the left border of 
the growth diagram. (In the direction from left to right one might prefer calling this a shrinkage diagram 
rather than a growth diagram; we shall call it a horizontally reversed growth diagram.) That S will be 
found is obvious, since we are considering the implicit shapes of the same sequence of matrices; that 
fact that R can be read off as indicated, follows from the fact that upward moves are exhausted in the 
submatrices Pnx(n-[j]) f° r au h whence one has row(P Nx ( N _ui)) = 7r(P Nx (N-[j])) = T (%x(N-[j]))- 
We illustrate this in diagram 3 for the binary matrix M of (28). 

Calling this a horizontally reversed growth diagram means that for any pair of indices k, I £ N, if we 
put A = 7r(M[ fc ] x ( N _[ i+1 ])), n = 7r(M[ fe ] x ( N _[j])), v = 7r(M[ fe+1 ] x ( N _[j +1 ])), and k = 7r(M[ fe+1 ] x ( N _[j])), 
there should be a relation between these shapes and the matrix entry M/-j that determines for fixed u., v 
a bijection between the pair (X,Mk,i) and the shape k. Indeed, if we determine the implicit shapes of 
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Diagram 3. Horizontally reversed growth diagram with binary matrix and shapes (7r(M[ i ] x (is r _[j])))i,jeN- 



submatrices by exhausting the pertinent raising operations (as in computing ordinary rather than French 
normal forms), then we have a situation in which everything is rotated a quarter turn clockwise with 
respect to that of the previous subsection. Taking our usual binary matrix M of (28) and (fc, /) = (5, 2), 
we can read off the shapes A = 7r(M[ 5 ] x ( N _[3])) = (3,2,1), /i = 7r(M[ 5 ] x ( N _[ 2 ])) = (4,2,2,1), and 
v = 7r(M[6] x (N-[3])) = (5,2,2) from diagram 3; from this and M 5y2 = we can compute k = (6,3,2, 1) 
as follows (the two leftmost columns that are of no interest here are represented by dots): 
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The rotation of the matrices induces the transposition of their implicit shapes and those of their subma- 
trices; therefore the relation here between A', v 1 , and k 1 (in that order) will be the same as was found 
between A, fi, v, and k in the situation of the previous subsection. It follows that one still has A /i k 
and A v- v k (which is also obvious from the computation considered), but the shape datum that 
determines k will be the coZwrnn-insertion shape datum for (V, [2]); its definition differs from the 

row-insertion shape datum by the fact that another injective map S — > T is used, one that maps each 
s E S to the first t G T in a column strictly to its right (and in a row weakly above it) . 
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4.4.1. Proposition. The bijection of theorem 3.1.3(1) between M € Ai l2] and the pair (P,Q), 
or the pair of tableaux (i?, S) encoded by them as in proposition 3.3.3, can be computed using a 
horizontally reversed growth diagram for the column-insertion shape datum for (V, , t— , [2]), reading 
off R = (7r(M Nx ( N _y]))j e N along the lower border and S = (n(M^ xnWn along the left border. □ 

We may interpret proposition 3.3.3 as a way to determine such horizontally reversed growth diagrams: 
the transition from the semistandard Young tableau formed by the shapes along column j + 1 of the grid 
and the one formed by those along column j is given by column insertion in strictly increasing order of 
the letters of a word w whose weight is given by the column M|. This may be verified by comparing the 
intermediate tableaux in (34) with the vertical sequences of shapes in the display above. Computing the 
shape k at the bottom left corner of a square (i, j) amounts to finding the shape occupied by letters < i 
after column-insertion of w, given the shape /i occupied by letters < i after that insertion, the horizontal 
strip v/\ filled with letters i before the insertion, and the number M j; j e {0, 1} of letters i in w. 

We have seen how n(M) can be computed using either an ordinary or a horizontally reversed growth 
diagram, which were determined by studying the exhaustion of appropriate crystal operations on sub- 
matrices from the top left respectively from the top right, and the two situations are the same up to 
a rotation of everything by a quarter turn. It should by now be clear that two more variations to this 
construction can be obtained by further rotations. Thus we obtain computations using a vertically and 
a doubly reversed growth diagram, which are respectively governed by the column-insertion an by the 
row-insertion shape datum for (V, , ^— , [2]). The relations this gives between the results of applying 
the corresponding variations of the dual RSK algorithm to M are both straightforward and in principle 
known (see for instance [Fult]), so we shall not give more details here. However, for integral matrices 
similarly defined variant growth diagrams are not all related by simple symmetries, as we shall now see. 

For integral matrices M £ M. we have seen that any upper left submatrix M[j] x ui can be reduced 
inside M to its normal form by exhausting the applicable upward and leftward crystal operations. If one 
views M as a finite rectangular matrix, then one can similarly bring any submatrix containing one of the 
three remaining corners into a standard form by exhausting an appropriate set of crystal operations, but 
the result is not the same in all cases. In the case of a lower right submatrix the situation is rotated a 
half turn with respect to that of an upper left submatrix, so exhaustion of all downward and rightward 
crystal operations between pairs of rows respectively columns of the submatrix will end with a rotated 
version of Diag(A): all units are on the diagonal passing through the bottom right corner, their numbers 
weakly increasing towards that corner. This rotated version has the same implicit shape A as Diag(A), as 
can be verified by any of the many ways we have seen to compute implicit shapes. Therefore the family 
of shapes (M(isf_[j]) x (isr_[j]))j 1 j£N defines a doubly reversed growth diagram in, which for all k, I G N the 
four shapes A = 7r(M( N _[ fe+1 ]) x ( N _[ ;+1] )), fi = 7r(M (N _[ fe+1 ] )x ( N _[ ; ])), v = 7r(M( N _[ fe ]) x(N _[ i+1 ])), and 
k = 7r(-M(N-[fc])x(N-[m) at the corners of a grid square and the entry M)~j of that square are related by 
the shape datum of the Burge correspondence, that can be computed by algorithms (40) and (41). 

For submatrices containing the top right or bottom left corner of M however, the situation is not 
symmetric to those of the top left and bottom right corners; on the other hand these two remaining cases 
are related by a rotation by a half turn, so we shall focus (arbitrarily) on the bottom left submatrices. 
A normalised form for the submatrix M(N-[fc]) x [;] that can be obtained by crystal operations acting 
on the entire matrix, is the one where upward and rightward crystal operations between pairs of rows 
respectively columns of the submatrix are exhausted, in other words the operations in { e\ i E N — [k] } U 
{ ff I j £ [I — 1] }• The result is not however that all units of the submatrix end up on one same diagonal 
or anti-diagonal, as we illustrate for our usual integral matrix M and its submatrix (k,l) — (1,6): 
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One sees that entries Mij with e (N — [k]) x [I] have become constant along diagonals, and the 

implicit shape of M( N _ m can be read off from the row sums of the submatrix after the transformation: 
in (45) on has 7r(M( N _[ 1 ]) x [ 6 ]) = (5,5,2,1). This must be so, for if one would isolate M( N _[ fe ]) x [;] and 
exhaust leftward crystal operations on it, then the result would on one hand be a shifted version of 
Diag(A) for A = 7r(M( N _[fe]) x [;]), an d on the other hand it would have unchanged row sums. Whereas 
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the Xi units forming such a row sum are all grouped together in one entry of Diag(A), they are distributed 
in the corresponding row of M^-[k])x[l] according to the partitioning of the squares of row i of [A] by 
the lengths of their columns: the unit corresponding to a square lies on the diagonal near the top right 
of (N — [k]) x [I] that has the same length as the column of the square. This is easily checked, since the 
entry M k+i j_i lies at the end of the row of Mrn-[k])x[i] with sum Aj and also of a diagonal of length i + 
and its value is Mk+%,i-i = A, — Aj+i, which is the number of columns of length i + 1 oi the diagram [A]. 
Thus the submatrix shows how that diagram can be sliced vertically into parts for each column length; 
for instance in the example the diagram [(5, 5, 2, 1)] appears as cut up into two single columns, of length 
4 and 3, and a group of the remaining 3 columns of length 2, while there are columns of length 1. We 
now define more formally the "sliced form" matrix that thus corresponds to A, for given k, I. 

4.4.2. Definition. Let A G V and k, I G N with A; = 0. We denote by SF(A, fc, I) the integral matrix M 
defined by Mij — [i > fc] [j < I] (A TO _i — A m ) for i, j G N, where m = i — k — j + l. 



M, 



(N-[fc])x[Z], 



4.4.3. Proposition. Let k,l G N with I > be fixed; then M G M satisEcs M 
n\{M) = for i > k, and nj*(M) — for j G [I — 1], if and only if there exists A G V with A; = such 
that M — SF(A, fc, /). 

Proof. It is clear that any M — SF(A, fc, I) with A; = satisfies the given conditions: for that case 
all the contributions to the sums in (18) and (21) are < 0. Conversely suppose that M satisfies those 
conditions. If its entries were not constant along all diagonals within (N — [fc]) x [I], then on the lowest 
diagonal with non-constant entries one could find a pair M^j ^ M i+1 j +1 with k < i and j + l < I; 
if < Mi + \j + \ one would have n\{M) > while if Mjj > M i+1 j + i one would have nj'(M) > 0, 
either of which contradicts the assumption. Therefore the entries of M are constant along such diagonals; 
moreover n\(M) = implies Mj+^o = for all i > fc, so the entries are zero on all diagonals below the one 
containing Mk,o- Then putting Ai = X^e[z] Mk+ij one checks that A G V, Xi — 0, and M — SF(A, fc, I). □ 

Now consider an integral matrix A = ^4(N-[fc])x[i+i] that satisfies n\(A) = nj~(A) = for alH > fc + 1 
and j < I — 1, so that in particular its submatrix A N _r fe+1 i) x r;i is of the form just described. Let us try 
to find the relation between the entry Ak,i just above and to the right of that submatrix, and the shapes 
at the grid points forming the corners of its square: A = 7r(j4(isf_[/s+i])x[i])i A 4 = 7r (^4(N-[fe+i])x[i+i]), 
v = 7r(j4( N _[/c])x[i])i and n — 7r(j4( N _[fc]) x [i+i])- The matrix after the transformation in (45) provides an 
example with (fc, I) = (0, 6), Mk y i = 0, and A = (5, 5, 2, 1): the computation below shows that in this case 
H = (8, 5, 5, 2), v = (8, 5, 3, 1, 1), and k = (8, 8, 5, 3, 1). 
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(46) 



We shall now analyse in detail what happens during this procedure. The description will be heavy 
going at times, but it eventually leads to the shape datum for the RSK correspondence; as this forms the 
most direct link of our constructions to that correspondence, we insist on presenting it, with its inevitable 
technical details. The reader may choose to either bear with us, or take our word for it and skip to (51). 

The parts of the transformations above that respectively reduce A^-[k])x[l] to SF(v, fc, I) by upward 
transfers (the part to the left of the vertical line for the left upward arrow), and ^4(N-[fc+i])x[i+i] to 
SF( / u, k + 1, 1 + 1) by rightward transfers (below the horizontal line for the bottom arrow), are less 
transparent that the corresponding parts in (38), but we can describe the former as follows (the latter is 
quite similar). From nj'(A) = it follows that A^j < Ak+ij+i for j G [I — 1]. Therefore, if one starts 
by exhausting n\ applied to A, the entry Akj will block as many units of Ak+ij+i; any further units 
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of Ak+ij+i however will not be permanently blocked, and they will eventually (when all transfers further 
to the right are completed) be transferred one row up to the entry Akj+i- As a result of this exhaustion, 
every entry Ak+ij+i will have become equal to the initial value of A k j, and A k+ i fi will have become 0. 
We illustrate this for the example above, showing just the first 6 columns of rows 0, 1: 



1 1 3\ (ej)* f 1 2 2 3 
001130/ ' I 1 1 



(47) 



After this exhaustion, the situation in the first / columns of rows k + 1, k + 2 is the same as the original 
situation in rows k, k + 1, but shifted one place to the right so that the final entries disappear, and with 
entries shifted in at the left. One can therefore continue by exhausting n k+1 , which will involve shifted 
copies of the same moves (except those that took place in the final column), and then similarly exhaust 
n k+2 , ■ ■ ■ , n l+i-i: after this, the entries in A(N-[fe])x[i] wm have become constant along diagonals, so that 
submatrix is then equal to SF(^, k, I). Thus we see that row k of the submatrix will already have obtained 
its final value after the initial exhaustion of el, and v can already be read off from this row. 

Now for any j > 0, the initial value of Ak+i,i-j represents the number Aj_i — Xj of columns of length j 
of [A], and the final value of the entry A k j-j above it represents the number of columns of that length 
of [v\. Then the units that are transferred upwards from A k+ \j_j to A k ^_j represent columns whose 
length is j in [A] and remains so in [u] , while the units that already were in A k j_j initially (and remain 
there), and whose number equals the number of units that remain in A k+l i _^_ 1 s ) if j > 1, represent 
columns whose length grows from j — 1 in [A] to j in [v] . The latter columns are the ones that for which 
row j — 1 of [v/X] contains a square. From this one deduces that vjX is a horizontal strip, and that the 
initial value of the reverse (A k ,i-i, ■ ■ ■ , A k ^) of the topmost row of ^4(N-[fc])x[i] gives the value of v — X. 
In the example above one has v - X = (8,5,3, 1, 1, 0) - (5, 5, 2, 1, 0, 0) = (3,0,1,0,1,0) = (A , 5 , . . . , A ,o)- 

We might have reached the same conclusion with somewhat less effort by considering column sums, 
which are unaffected by upward transfers, but the details will be needed anyway to understand the 
full transformation of A by upward moves, notably the changes to column We now know that this 
column initially contains, below the entry A k j, a vertical copy of the composition /i — A, and after the 
transformation it will contain a similar copy of k — v (but shifted one place less, starting at square (k, I)). 
When we consider the exhaustion of some e k+i , we can focus exclusively on columns I — 1 and I, since 
what we have seen above implies that at that time one has A k+i j < A k+i+ ij + i for all j 6 [I — 1], so that 
the presence of columns to the left of column I — 1 will not block any upward transfers in column I. 

The first 2x2 submatrix to consider is ^4{fc,fe+i}x{i-i,;} whose entries are initially (i.e., at the start 
of the exhaustion of e k ) the following: A k ,i-i is the final entry v n — A of the composition rcv;(^ — A) 
that is in the first I columns of row k; A k j is arbitrary; A k +i,i-i is the top-rightmost entry Ao — Ai of the 
submatrix ^4(N-[fe+i])x[i] = SF(A, k + 1, /); ^4fc+i,; is the top entry /xo — Ao of the vertical copy of fi — X in 
column I. At the start of the exhaustion of e k+i , the values in the 2x2 submatrix j4.{fe+» 1 fe+»+i}x{i-i,i} 
are similarly: A k +i,i-i = v% — A, is a copy of the initial value of A k j-i-i; A k+i j is some value a left after 
the exhaustion of e k , i _ 1 (which we shall determine presently); A k+ i + u-i is the final entry Ai — Aj+i of 
row k + 1 + i of SF(A, k + 1, 1); and finally A k +i+i,i = Hi — Aj. The upward moves in column / are then 
given by 

'^o-Ao A k j V^ )max( " r" 0, °V - A A kt i+max(jio-vo,0)\ (4g) 



A - Ai ho^XqI \X Q -X 1 min(^ , ^o) - A 



for the transfer between rows k and k + 1 (which will be followed by a transfer in column I — 1 changing 

tVQ — Vl 



it from (^°_^°) into (^°_^) , but that part was already known), and similarly in rows k + i and k + i + 1: 



Vi — Xi a \ ( e fc+i) maxf ^ "*' 0) / Vi-Xi a + max( / u i - u t , 0) 

Xi - X i+1 Hi- Xi J I Aj - A i+ i min(/i i , vA - X t 



(49) 



It is now clear that a = min(/i i _ 1 , Vi-i) — Aj_i, and one easily checks that for all i > 0, after completion 
of the exhaustion of e k+i and of the subsequent exhaustion of e£ +i+1 , the submatrix A^ k+i k+i+1 y x ^_ 1 ^ 
will have become 

Vi — Vi+\ a + max(^ii — fi, 0) 

Vi+\ - Vi+2 min(z/j, pti) - Xi + max(/i i+ i - 0) 



(50) 



Although max(^ i+ i — i/j + i,0) = max(/Xj + i, — v i+ \ units have been added to the bottom right 
coefficient by the exhaustion of e k+i+1 , we do not expect that this re-enables applications of e k+i (since 
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such applications were not required for reducing the submatrix ^4(N-[fc])x[/] to SF(i/, k,l) cither); indeed 
it follows from Ai > max(/ij + i, v-i+i) that one has v; t — Vi+i > min^j, ^) — Aj + max(/ij + i, fj+i) — fj+i, 
whence further upward transfers are blocked. We may conclude that after the mentioned sequence of 
operations we will have determined the new value of column I representing the composition k — v. 

In spite of the complications of the reduction, the resulting relation between A, fi, u, k, and Ak.i is 
surprisingly simple, as there is no "carry over" effect like there was in (39): the value of a has no effect 
on the bottom-right entry in (50) that determines /tj + i — z^+i. Therefore we can describe the parts of k 
by the following explicit formulae rather than by an algorithm: 



«o = A k j +max(p ,is ) 

= nrm(z/j,/ij) - A, + max(/z i+ i, v i+1 ) 



for all i G N. 



(51) 



Comparison with [vLee4, equation (15)] shows that this is the shape datum for the RSK correspondence, 
which reflects row-insertion on semistandard tableaux. 

One can give the following interpretation of (51) in traditional terms. If one imagines, during RSK 
insertion, a semistandard Young tableau S in whose display the entries < I fill the diagram [A] while 
the entries I fill A], and if some insertion of letters < I into £ makes the shape filled by entries < I 
grow to [u], then upon subsequent row-insertion of Mjy letters I, these will be placed in row starting 
at position max(^ , ^o) (i- c -, after the unbumped entries I of S, if any of them remain, or else after the 
entries < I present at that point) making the length of that row max(/io, ^o) + ^k,i- That insertion will 
not bump any letter I, but the previous insertion of letters < I will bump min(^, /i^) — Ai letters I from 
row i of S, and they will be placed in row i + 1 starting at position max(/i i+ i, ^,+i), making its length 
max(/ij_|_i , Vi+i) + min(i^, /Ltj) — A^. Note in particular that it is the units that are not transferred upwards 
from the entry at (i + k + 1, 1) in (49) that represent letters I being bumped from row i to row i + 

Now that we have determined the shape datum, the family of shapes (M( N _[i]) X [j])ijeN can be 
computed using a (vertically reversed) growth diagram, just like we have done before for other families 
of implicit shapes of submatrices. This time in fact we obtain a classical RSK growth diagram, growing 
from the bottom left corner of the matrix. We illustrate this in diagram 4 for our usual matrix M. 
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Diagram 4. Vertically reversed growth diagram with integral matrix and shapes 
(7r(M( N _[i]) x [j]))ijeN, computed using the RSK shape datum. 

4.4.4. Proposition. If for any M G A4 one associates to each grid point G N 2 the shape 

7r(M/isf_fj|) x ui), then one obtains a vertically reversed growth diagram for the RSK shape datum. □ 

Along the top border of diagram 4 we find the semistandard Young tableau S of (31); this is no 
surprise since we already knew that S = (7r(M Nx y]))j e N, which was computed along the bottom border 
in diagram 1. Along the right border we find a reverse semistandard Young tableau (7r(M( N _[ i ]) xN )) iG N 
that we did not encounter before. However, the way in which S, whose integral encoding is the matrix P 
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of (29), was found here, corresponds to the description of theorem 3.3.1; according to comments following 
that theorem, the recording tableau of that insertion is the Schiitzenberger dual of the semistandard 
Young tableau L encoded by the matrix Q l . Indeed the Schiitzenberger dual of 
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which is the reverse semistandard Young tableau that can be read off along the right border of diagram 4. 



4.5. The Schiitzenberger involution. 

We have encountered the Schiitzenberger involution in two similar situations, for binary and for integral 
matrices. In the former case it defined a bijection mapping the reverse transpose semistandard Young 
tableau R = (7r(MN X (N_[j])).j e N ^° ^ e transpose semistandard Young tableau R* = (7r(MNx[j]))jeN; 
in the latter case it gave a bijection transforming the semistandard Young tableau L — (7r(M[j] xN ))jgN 
into the reverse semistandard Young tableau L* = (ir(M(p ! -[i])x'N))iEN- The fact that the members of 
these pairs of tableaux mutually determine each other, independently of M, may seem surprising at first 
glance, since at corresponding positions in their sequence they contain implicit shapes of complementary 
submatrices of M, which shapes are completely unrelated. On the other hand, these pairs of tableaux 
correspond to matrices obtained from M by exhausting crystal operations in opposite directions. We 
have seen that R and R* correspond to the matrices P and P obtained from the binary matrix M by 
exhausting upward respectively downward operations (cf. propositions 3.3.3 and 4.3.3). In the integral 
case L corresponds to the matrix Q obtained from M by exhausting leftward operations, and one easily 
sees that L* similarly corresponds to a matrix Q obtained from M by exhausting rightward operations 
(up to some sufficiently large limit); in the case of the example this is 



Q 



/0 1 1 3 0\ 
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Vo 3/ 



(53) 



whose columns from right to left give the weights of the rows of the display of L* from top to bottom. 
Since for each of the pairs of matrices (P, P) and (Q, Q) the members mutually determine each other 
(by exhaustion of appropriate crystal operations), it is clear that the same must be true for the pairs of 
tableaux (R, R*) and (L,L*). Exhaustion in the remaining pairs of opposite directions also give pairs of 
matrices that correspond to tableaux related by the Schiitzenberger involution. 

We shall not study the Schiitzenberger involution in detail here, and just observe some facts without 
formal proof. Its original definition in [Schiil] for standard Young tableaux can be generalised, using 
appropriate notions of standardisation, to semistandard or to transpose semistandard tableaux; doing so 
one must map tableaux of weight a either to tableaux of weight rev; (a) for some sufficiently large I, or 
to reversed tableaux of weight a (while keeping the semistandard or transpose semistandard attribute). 
Above we have implicitly assumed the latter, more elegant, form of these generalisations. 

The original definition of the Schiitzenberger involution is by an algorithm that repeatedly applies a 
deflation operation (also called evacuation) to a standard Young tableau S, which removes a square from 
the diagram of its shape at every step, until the shape becomes (0); the resulting decreasing sequence 
of shapes defines the Schiitzenberger dual S* of S. It is natural to index the sequence defining the 
tableau deflated i times by N — [i]. If one then writes = (A^ J ')j G N-[i]j then the definition 
of the deflation procedure amounts to a simple rule for determining \( l+1 >j +1 ) once A^ +lj \ A^ J \ and 
are known. Thus S = (A^) )™ =0 determines S* = (A^'™)) ™ =0 via a two-parameter family of 
shapes governed by a local rule, in much the same way as a growth diagram, but without matrix entries. 

The generalisations to semistandard or to transpose semistandard tableaux can be formulated in 
terms of deflation by an entire horizontal respectively vertical strip at a time, and give rise to similar 
families of shapes (A^'^)^;,^, in which neighbouring shapes now differ by horizontal strips respectively 
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by vertical strips. Here the rule for determining \( l+1 >i +1 ) given A^ +lj \ \( l,: >\ and \( l >j +1 ) can be 
formulated as follows, using the notion of tableau switching for semistandard tableaux given in [vLee3, 
definition 2.2.4] (the transpose semistandard case is similar): applying tableau switching to the single-strip 
skew semistandard tableaux T = (\( l+1 >ti ■<— A^)) and 7\ — (A^'- 7 ' t— A^ ,:,+1 )) gives similar tableaux 
(T{,Tq) = X(T ,Ti) where Tq = (A^ +1 ^ v- and T{ = (\( i + 1 '3+ 1 '> v- A^' +1 )). An analysis 

of what tableau switching means in this particular case leads to the following alternative formulation: 
the two skew semistandard tableaux of two horizontal strips each, \( %+1 >i) ^— \( l >fi t— \( 1 > J+1 ") and 
,_ \(. l + 1 -J+ 1 ) ,_ A^ J+1 \ are related to one another by a Bender-Knuth involution (the unique 
applicable one). Interestingly Bender-Knuth involutions, which we preferred to replace by ladder reversal 
involutions for establishing the symmetry of skew Schur functions to get more pleasant properties, thus 
still have their proper place in the theory. 

The way in which we have encountered the Schiitzenberger involution in fact relates more naturally 
to a different way to compute it. The matrix P obtained from M € M. by exhausting upward crystal 
operations as in (29) encodes a unique semistandard Young tableau S, of shape 7r/(0) where n = n(M). 
If k > is so large that M = M[ fc ] xN , then the matrix P obtained from M by exhausting downward 
crystal operations // for i E [k — 1] determines the Schiitzenberger dual S* = {fJ-^)jeN OI S, namely 
by /j,tf> = 7r(P Nx (N-[j])) = revfc(row(Pisr x (isr_rj])))- But if one fixes / > 7To, and p 6 V is such that 
[p] = W\ x [']> then P is also the integral encoding of a unique semistandard tableau of shape p/n*, where 
7r° = p — rev fc (7r), namely of S"* = (M^^)jeN = (p — T ^' v k(P'^))je'N- This is the tableau whose display 
is obtained from that of S* by the rotation by a half turn that maps the rectangle [p] to itself. 

The semistandard tableaux S and S*° are encoded by integral matrices P and P that can be obtained 
from one another by exhausting upward respectively downward crystal operations, so it follows from 
proposition 2.5 that S is the rectification by jeu de taquin of S*°, and since rotation by a half turn 
transforms inward jeu de taquin slides into outward slides, 5**° is similarly the "bottom-right rectification" 
of S within the rectangular shape p. In our example one may check, for (k, I) = (5, 8), that 
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and that the rectification of S*° is the tableau S of (31); thus S <-> S* can be computed via rectification. 
The same relation exists between L and L*° in (52), and between R in (35) and obtained from R* 
in (44); in the latter case a version of jeu de taquin for transpose semistandard tableaux is used. 

This alternative method of computing the involution was first presented by Schiitzenberger in [Schii2] . 
It is not so easy to deduce it from the original definition of the Schiitzenberger involution, but the opposite 
is straightforward, once it is established that the rectification of any semistandard tableau T of shape pj A* 
is of shape A. Then the shape pffl in S* = (p^)j^ can be determined as the shape of the rectification 
of the skew tableau obtained by removing j initial horizontal strips from S, and these shapes can be 
found by successively deflating S by one horizontal strip at a time. By construction that rectification is 
given by the sequence of shapes (A^'^)^ , where (^'^)j,l;j<l<n is the semistandard tableau switching 
family used to compute S* = (A (j '™))™ =0 starting from S = (A ( ° ! ^)™ =0 . 

Now if P is the integral encoding of S, then Pnx(n-[j]) is the integral encoding of for all j G N, 
and the integral encoding of the rectificationof can be found by exhausting upwards crystal operations 
applied to Pnx(n-[]])- This gives have the following interpretation for the shapes in the semistandard 
tableau switching family: A^'^ = n(M-yi x ^_y^), which agrees with S* = (7r(Pisr x (N-[j])))jeN- Since ver- 
tical crystal operations do not change the implicit shape, (7r(M Nx ([;]_[j])))j,; ;: ,<i is a semistandard tableau 
switching family for any integral matrix M. The same is true for (7r(M(ru_ui) X ]y))i,fc;i<fc; these statements 
hold in the binary case as well, with the adjustment that the former family is transpose semistandard. 
However, in contrast with the situation for the various growth diagrams we considered, we have not been 
able to find a simple direct argument why a local rule for these families of shapes should exist, without 
referring to the entire family; for instance one showing that knowledge of the shapes 7r(M Nx ([;]_y +1 ])), 
7r(M Nx ([;]_y])), and 7r(M Nx ([; + i]_[j])) suffices to determine 7r(M Nx Q; +1 ]_y +1 ])). In principle one could 
establish the local rule for that family using those for (Mru x m)fe ; ; e isr and for (Mru x ( N _rj-i))fc ;:) - e isr, an d 
induction on k, similarly to the proof of [vLeel, theorem 4.1.1], but this is not very attractive. 
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§5. Characterisation of the implicit shape. 

The implicit shape of a matrix M characterises by definition the double crystal in which M occurs. We 
have by now seen many different ways to compute this shape, but these involve either transforming M 
into a more special matrix within its double crystal (by exhausting some operations), or constructing 
a growth diagram from M. It would be convenient to have a characterisation of ir(M) in terms of M 
itself, which would provide a quantity invariant on double crystals and distinguishing any two distinct 
double crystals. In this section we shall present such a characterisation of n(M), which will however 
usually not provide a practical way to compute it; the methods given above can therefore be interpreted 
as algorithms to compute the invariant that we shall describe below. 

It is well known that the length of the longest increasing and decreasing sequences extracted from 
a permutation a (represented as a word) are given by the lengths of the top row and leftmost column 
of the tableaux associated to a under the Schensted correspondence; extending this result, it was shown 
in [Grel] that the entire common shape of those tableaux can be characterised in terms of the permutation 
matrix (which is viewed as describing a partially ordered set). A version of Greene's poset invariant can 
also be defined for integral and binary matrices, with a similar interpretation as the shape of the tableaux 
associated under the RSK correspondence and its dual. Therefore we shall not waste any time trying 
to derive the invariant we are looking for, and simply present the proper formulation. (Note that the 
word invariant is used in two different senses here: the partition defined by Greene depends only on the 
abstract finite poset defined by the permutation, but its extension to matrices is invariant under crystal 
operations, even though these do in general alter the associated poset.) 

5.1. Definitions. Let M g M [2] be a binary matrix and N g M an integral matrix, and m g N. 

a. M is called a chain if M k j = 1 implies Mij = for all i < k and j > I, and it is called an antichain 
if Mk,i = 1 implies Mjj = for all i < k and j < I. 

b. N is called a chain if N k ,i > implies Nj = for all i < k and j < I, and it is called an antichain 
if Nk,l < 1 f° r &11 k, I g N, while N k ,i = 1 implies Nj = for all (i,j) ^ (k, I) with i < k and j > I. 

c. An m-chain or an m-antichain in M is a sequence (M^ fe ')j, e [ m ] of m binary matrices, all chains 
respectively antichains, such that M — X)fee[m] + M' for some M' g M l2] . 

d. An m-chain or an m-antichain in N is a sequence (N^) ke [ m ] of m integral matrices, all chains 
respectively antichains, such that N — J2ke[m] + ^' f° r some N' E A4. 

e. The m-length of a binary or integral matrix A is the maximum of X)fee[m] l-^ fe ^l> taken over all 
m-chains (A^) ke [ m ] in A, where \A\ — X^jeN^j- ^ e m-breadth of A is the maximum of 
Sfee[m] l-^^l) taicen over all m-antichains {A^) ke ^ n ^ in A. 

The terminology is derived from the partially ordered sets associated as follows to the matrices. 
For M g M l2] an element x^j is associated to every bit Mjj = 1, and Xij < x k ,i holds if and only if 
i < k and j < I. For N g M. a chain of Nij elements is associated to every entry Nij, and for every 
element x of the chain for Nij and every element y of the chain for another entry N k ,i one has x < y if 
and only if i > k and j < I. In this setting it is convenient to use some more related terminology. For 
m-chain or m-antichain (A^) ke [ m ], either binary or integral, we shall call X)fee[m] \ its size. If A, B 
are binary or integral matrices, the matrix (mm(Aij,Bij))ij e jq will be called their intersection A n B; 
in case B = A n B we shall write A D B and say that A contains B. In particular A contains the sum 
of any m-(anti)chain in A. Then the intersection of a chain and an antichain has size at most 1, and it 
follows that the sum of an m-chain cannot contain any antichain of size > m, and that the sum of an 
m-antichain cannot contain any chain of size > m. The converse is also true, namely that any matrix 
that does not contain any antichains of size > m is the sum of an m-chain, and one that does not contain 
any chains of size > m is the sum of an m-antichain, but this is somewhat less obvious. 

The definition is formulated so as to emphasise the similarities between the different cases, but there 
are notable differences as well. In the binary case chains and antichains can be interpreted as sets of 
squares where the matrix has bits T', and the sums involved in m-(anti)chains obviously must correspond 
to disjoint unions of such sets. In the integral case antichains can also be seen as sets of squares, but the 
antichains in an m-antichain need not be disjoint, or even distinct. On the other hand one could require 
in the definition of m-length for the integral case that the nonzero entries in the chains of the m-chain 
involved occur in disjoint squares: the indicated maximum would not decrease due to this restriction. 

It is not hard to check both in the binary and in the integral case that a matrix M is a chain 
if and only if ir(M) has at most one non-zero part, in other words if n(M)i = 0, and that M is an 
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antichain if and only if all parts of 7r(M) arc < 1, in other words if 7r(M)J = 0. In general it is quite 
hard to determine, for a given M and m, an m-(anti)chain in M of maximal size. In the special cases 
M = Diag [21 (A) and M — Diag(A) however, this is easy. For M = Diag [21 (A), the rows of M define chains 
and the columns antichains, and the m-chain or m-antichain formed by the m longest rows respectively 
columns has the maximal possible size, since its intersection with each individual column respectively 
row has maximal size. For M — Diag(A), each chain contained in M has at most one nonzero (diagonal) 
entry, so an m-chain of maximal size contained in M is formed by extracting the m largest diagonal 
entries: {M^ x ^) ie ^ m y In this case the size of an m-antichain in M is at most ^ ieN rnin(m, M^j), and 
this size is obtained by taking the m-antichain (l(Aj))j £ [ TO ], where 1(7) denotes the I x I identity matrix 
([* = 3 G [I]] )ij'GN- Note that in both cases the m-length of M is J2ie[ m ] anc ^ lis ^-breadth is 
5Zie[m] ^i- ^ ne mam result of this subsection is that this is true in general, where A = ir(M). 

5.2. Theorem. For any binary or integral matrix M, and any m G N, the m-length of M is equal to 
J2ie[m] an d the m-breadth of M is equal to J2ie[ m ] where A = ir(M). 

Since we know the theorem holds for normal forms, it is equivalent to the statement that for every m 
the m-length and m-breadth of a matrix are invariant under all crystal operations, and this is the way in 
which we shall prove it. Note that one consequence of the theorem is that the sequence of all m-lengths 
of M determines that of all m-breadths and vice versa. It is known that this statement is true in general 
for finite posets, see [Grc2, Theorem 1.6], but a proof for that case is rather more difficult than for the 
situation of our theorem. 

There are numerous ways to approach a proof of the theorem, many of which start by reducing to 
a simplified situation. It is fairly easy to see that the argument by which the theorem holds for normal 
forms can be adapted to show that is holds whenever either upward or leftward moves are exhausted, 
so that a proof for operations in one direction would suffice (but this facilitates only the binary case). 
The result is also known to hold for permutation matrices, see [Grel, Theorem 3.1], and it is possible 
to enlarge matrices, pulling their units apart so as to obtain a permutation matrix, while not changing 
their implicit shape or m-length and m-breadth; this would provide a proof. One can define a notion of 
chains in the display of a semistandard tableau T, and [Schu2, theoreme 3.1] states a variant of Greene's 
result for that case, where the shape is that of the rectification of T; it can be seen that the m-length 
and m-breadth so associated to T are the same as those of its integral encoding, so this would provide a 
more direct proof for the integral case (but since the chains of T do not correspond to those of its binary 
encoding, more work would be needed for that case). Finally one could try, as is done in the proofs of 
the results cited, to decompose the changes to the poset structure into minimal modifications involving 
only three elements, and prove that these do not effect m-length or m-breadth. 

All of these methods however require arguments to justify the reduction, which are not related directly 
to the definitions of m-lcngth or m-brcadth, or of crystal operations; they also depend on nontrivial results 
about the Robinson-Schensted correspondence. We believe that for understanding why the theorem is 
true, it is most instructive to use a head-on approach, showing directly that the definition of crystal 
operations prevents them from changing any m-length or m-breadth; in fact this is not too hard to do, 
although the situation is more complicated than in the 3- unit case. By using only the definition of crystal 
operations, this approach could even serve as a starting point for developing their theory (but we believe 
this would not would make that development more transparent). We also choose a "positive" point of 
view, characterising matrices that are the sum of an m-(anti)chain systematically by the existence of such 
a decomposition, rather than by the non-existence of any antichain (respectively chain) of size m + 1. 
The proof we give will consist of a series of lemmas that each prove a part of the mentioned invariance. 

The number of cases to be treated in principle is considerable, since we need to show both in the 
binary and the integral case that neither the m-length nor the m-breadth can increase nor decrease 
under crystal operations in any of the four possible directions. We reduce this number using the usual 
symmetry considerations to four distinct cases, which can moreover be regarded as two pairs of analogous 
cases. Both in the binary and integral cases, rotation by a half turn leaves both m-length and m-breadth 
invariant, so for each of these quantities and any pair of opposite directions it will suffice to show that the 
quantity will not decrease under a crystal operation in one direction: by rotation symmetry the inverse 
crystal operation will then not decrease the quantity either, so that it is invariant under both operations. 
In the binary case, rotation of matrices by a quarter turn transforms m-chains into m-antichains, which 
permits us to consider only the m-length; in the integral case we must consider both m-length and 
m-breadth, but transposition symmetry will allow us to consider transfers in one direction only. 
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The binary and the integral case are sufficiently different to merit separate treatment. We shall start 
with the somewhat simpler binary case, where as said we shall consider only the m-chains. Much if our 
discussion will be about locating bits Mjj = 1; we shall call this a unit at position (i,j), and we shall 
view the individual chains of an m-chain in M as sets of such units. 

5.3. Lemma. Let a binary matrix M G Ai l2] contain an m-chain of size n for some m, n G N, and let 
M' = e[(M) for some k G N. Then M' also contains an m-chain of size n. 

Proof. Let Mk,i = and Mk+i.i = 1 be the bits interchanged by the application of e\. One may assume 
that the unit at (k + 1,1) occurs in one of the chains C of the m-chain, for otherwise one can simply keep 
the same TO-chain. If, after removing the unit at (k + 1, 1) from C, any of the chains of the m-chain would 
remain a chain when a unit at (k, I) is added to it, then an m-chain of size n in M' can be formed in that 
way; we shall show that after possibly changing the choice of the m-chain of size n in M, we can always 
arrange to be in that situation. Let the units of C to the left of (k + 1, 1) in row k + 1 be at positions 
{ (k + | i G [r] } where r G N and jo < ■ ■ ■ < j r -i < I- If r = (i.e., if there are no such units) then 
after removing the unit at (k + 1, 1) we can add a unit at (k, I) to C itself, and we are done. 

Otherwise definition 1.3.1a gives us X^=j Mk.j > 2j=j Mk+i.j > T, so there are at least r different 
indices j < I such that Mk.j — 1; let the r largest such indices be j' < ■ ■ ■ < j' r _\. From the fact 
that 53j=ij > Ylj^ji Mk+i,j for i G [r], it follows that ji < j[ for all i. None of the units 

M/. ji = 1 can occur in C because of its unit at (k + l,jo)- If none of them occur in any of the 
other chains of the m-chain either, then one can modify C by removing the units at {k + and 
adding units at (k,jl), for i G [r]; this reduces us to the case r = already treated. If at least one 
of the units M^ji = 1 does occur in a chain of the m-chain, let j' s be the maximal index for which 
this happens, and let the unit at (k,j' s ) occur in the chain C . If C contains any units in row k to 
the right of the one at (k,j' s ) then this must be in a column > I, and a unit at (k + 1,1) can be 
added to C . If not, then we replace C, C as follows by a new pair of chains D,D' of the same total 
size: the chain D contains the units of C in columns < j' s , the units at {k,j^) for s < i < r and 
the units of C in columns > I, while D contains the units of C in columns < j s and the units of C in 
columns > j' s . After this change D takes over the role of C, and we are in the case r = 0, so we are done. □ 

This situation of the final case of the proof is schematically represented in figure 1, in which black 
circles represent units, the white circle is the bit Mk,i = 0, the red and blue lines indicate the chains 
C and C", with dotted parts that do not occur in D with the unit at (k + 1, 1) replaced by one at (k, I), 
nor in D'; the green lines indicate the new parts of the latter chains and the yellow dotted lines show the 
correspondence (k + 1, ji) <-> {k,j[). 




Figure 1. Schematic view of an partial exchange of chains for a vertical move in a binary matrix. 

5.4. Lemma. Let a binary matrix M G A4 l2] contain an m-chain of size n for some m, n G N, and let 
M' = ej~(M) for some I G N. Then M' also contains an m-chain of size n. 

Before stating the somewhat technical formal proof, we shall give an informal reasoning leading to 
it, ordering the cases that can occur by increasing complication. Let = and Mk^i+i — 1 be the 
bits interchanged by ep. First consider the relatively simple case m = 1. If the unit at (k, I + 1) does not 
occur in the given (l-)chain C, then C will still be a chain in M' . If the unit does occur in C, we can 
just replace it in C by the unit that was moved to the left, unless C already contains a unit in column I, 
necessarily at some (k',l) with k! < k. The inequality X^=fc' -^U — Y^i=k' ^h.i+i from definition 1.3.1b 
shows that there exists a unit at (fc", I + 1) with k' < k" < k; then by replacing in C the unit at (fc, I + 1) 
by the one at (k" , I + 1), one obtains a chain of size n in M that is also in M' . 

In the case m > 1 there is at most one chain C in the m-chain with a unit at (fc, I + 1), and if so 
one can apply the reasoning above to C, in order to find another chain to replace it in the m-chain. The 
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reasoning proceeds as before, but at the end there is a complication in case the unit at (k" , l + l) occurs in 
another chain C of the m-chain. If C contains no unit in column I, then one can resolve the difficulty by 
"genetic recombination" of C and C (a crossover of chains somewhat simpler than the one in figure 1): 
the unit at {k' , I) and any other units of C in columns < I combine with the unit at (k" , / + 1) and any 
other units of C in columns > I + 1 to form a new chain D, while the remaining units of C , which 
are in columns < I, combine with the unit at (k, I + 1) and any other units of C in columns > I + 1, to 
form another new chain D', which does not contain a unit in column I. This situation is schematically 
represented in figure 2, using the same conventions as in figure 1. If on the other hand C does contain 
a unit in column /, then one can repeat the same kind of argument to find another unit in column l + l, 
and so on; eventually there must be such a unit that is not in a chain with another unit in column I, 
and all the pieces of the chains involved can be pieced together to form an m-chain in M' (unlike for the 
previous lemma, any number of chains can be involved in the modifications here) . In order to proceed 
in an orderly fashion the formal proof shall start with untangling any chains that might be involved. 
Figure 3 schematically shows the two types of modification that may be required at the end of the proof 
below to reorganise the chains in the most difficult cases: a shift by several units in column I + 1 between 
chains, or a cyclic recombination of parts of several chains. 




Figure 2. Schematic view of a simple recombination of chains for a horizontal move in a binary matrix. 

Proof. Let the move effected by ej interchange bits Mk,i = and Mk,i+i = 1. Consider the subset S 
of chains in the m-chain that have a unit both in column / and in column l + l, and put r = #S > 0. 
Let / be the set of rows of those units in column I, and I' the set of rows of those units in column l + l. 
Since each chain of S contributes an i G I and an i' £ /' with i < i' , it follows after sorting that if 
I = {io, . . . , i r -x} and I' — {i' , . . . , i' r _x} with zq < • ■ • < V-i an d *o < ' ' ' < C-ij then ij < i'j for all 
j G [r]. Then we may assume (after recombining parts of the chains if necessary) that the units Mi jt i = 1 
and Mji j+i = 1 occur in the same chain C^' G S, for all t G [r]. 

If k ^ I', then either none of the chains of the m-chain has a unit at (k, I + 1) so that one can 
keep the m-chain, or that unit occurs in a chain without unit in column I so that it can be replaced 
in that chain by a unit at (k, I). Therefore assume that k = i' t for some t G [r], so that the chain C^> 
contains units at [k, I + 1) and at (it, I)', since Mk,i = one must have it < k. Now for any j < t, 
definition 1.3.1b gives 53i=;- Mij < $2»=i- ^i,l+i> an d each of the chains . . . ,C^* -1 ^ contributes 
equally to both sides, while contributes only to the left hand side. Then the right hand side either 
has a contribution not from any chain in S, or one from C^ -1 ' . The latter does not happen for j = 0, so 
there exists some i < k with M^+i = 1 and i ^ now put k~ = max{ i < k \ M^+i — l,i ^ I' } and 
s = min { j < t | ij > k~ } < t. The argument given shows that i s < k~ , and that ij < i'j_\ for s < j <t. 

Now if the unit at (k~ ,1 + 1) does not occur in any chain of the m-chain, then one can replace in 
the unit at (i' s ,l + 1) by one at (k~,l + 1), and for s < j < t replace in C^' the unit at (ij,l + 1) by 
one at (i'j_ 1 ,l + 1); the m-chain of size n in M so obtained is also an m-chain in M 1 . Finally if some 
chain of the m-chain does contain a unit at (k~ ,1 + 1), the this chain does not lie in S, and one can 
"cyclicly recombine" that chain C with the chains C^' for s < j < t: any units of C in columns < 
are combined with the unit at (k,l + 1) and any other units of in columns > I to form a new 
chain D; the units in columns < I of are combined with the unit at (kr , I + 1) and any other units 
of C in columns > /; and for s < j < t the units in columns < I of are combined with the units 
in columns > I + 1 of C" -1 ). We have obtained an m-chain in M of size n, in which the chain D 
contains a unit at (k, I + 1) that can be replaced by a unit at (k, I) to form an m-chain in M 1 of size n. □ 

We now consider the integral case, where the cases of m-chains and m-antichains must be treated 
separately, but moves in only one direction need to be considered. As we shall see, these two cases are 
like "thick" versions of those treated in lemmas 5.3 and 5.4, respectively. The thickness refers to the 
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Figure 3. Schematic view of a shift between chains by several units, 
and of a cyclic recombination of chains, for a horizontal move in a binary matrix. 

fact that many units can share a same position in the matrix; this seems to be a slight advantage in the 
case of chains (which are now more compact), but a disadvantage in the case of antichains (more effort 
is needed to tell them apart). To emphasise the indicated analogy, we shall deal with upward transfers 
for m-chains, and with leftward transfers for m-antichains. 

5.5. Lemma. Let an integral matrix M G M. contain an m-chain of size n for some to, n £ N, and let 
M 1 = el(M) for some ieN. Then M' contains an m-chain of size at least n. 

Proof. Let the transfer effected by e[ be from Mk+ij to Mjy. One may assume that one of the chains C 
of the m-chain contains a unit at (k + 1,1), for otherwise one can simply keep the same TO-chain. If any 
chain of the m-chain can accommodate a (new) unit at (k, I) while remaining a chain, then an m-chain 
of size n in M' can simply be constructed by removing a unit from C at (k, I + 1) and adding a unit 
at (k, I) to the chain in question. We shall show that one can always reach this situation by replacing if 
necessary the given m-chain in M by another one of size at least n. 

Let I' — max{ j | Ck+i,j > } > I be the column of the rightmost unit of C in row k + 1. If I' = I 
then C itself can accommodate a unit at (k,l). Otherwise, if none of the chains in the m-chain have 
any units at (k,j) for I < j < I' then we may replace C by a chain obtained by removing all its units 
at positions (fe + 1, j + 1) with I < j < I', while adding for those same values of j all Mj-j units of M 
at (k,j). The inequality X^=z iMk,j — Mk+i.j+i) > from definition 1.4.1a shows that the size of the 
new chain is at least that of C, so we have reduced to the case V = I. If some of the chains does have a 
unit at some (k,j) with I < j < V , then let jo be the minimal such index j, and C a chain of our m-chain 
with a unit at (k,jo)- If C has another unit in row k to the left of the one at (fc, jo), necessarily in a 
column < /, then it can accommodate an additional unit at (k, I). If C has no such unit, then we can 
replace the chains C, C by a pair of new chains D, D' formed as follows: the chain D contains the units 
of C in columns < I, all M^j units available in M at (k,j) for I < j < jo, and the units of C in rows < k, 
which are in columns > jo; the chain D' contains the units of C in rows > k, which are in columns < jo, 
and the units of C in columns > j . The inequality l^Lj (Mfcj — Mk+ij+i) > shows that replacing 
C, C by D, D' gives an m-chain of size at least n, and D can accommodate an additional unit at (k, I). □ 

We note that in the final case we could have included any units of C at (k + 1, jo) into D' , but the 
m-chain will be large enough without them; in fact we could have done without the units of C at that 
position as well. But there need not be any such units, and there will not be if n is maximal. 

5.6. Lemma. Let an integral matrix M 6 M. contain an m-antichain of size n for some m, n 6 N, and 
let M' = eJ~(M) for some I € N. Then M' also contains an m-antichain of size n. 

Proof. Let the application of e*~ result in a transfer from M^i+i to Mjy. Consider the subset S of 
antichains in the m-antichain the have a unit both in column I and in column 1 + 1, and put r = #S > 0. 
Let /, I' be the r-multisets on N such that the units in column I of the chains in S form the multiset 
§ (i, I) | i € I ]}-, and that the units in column Z + l of the chains in S form the multiset {[ (i, I + 1) \ i € I'J; 
write / = zo, ■ • ■ , ir-i} an d I' — f i' , . . . , iJ.-iS' with io < ■ ■ ■ < i r -i and i' < ■ ■ ■ < i' r _i- Since each 
chain of S contributes an index i to I and an index %' to /' with i < i', it follows that ij < i't for all j G [r]. 
We can then cut each chain in S into a part in columns < I and a part in columns > I, and recombine 
all the parts to form an r-antichain (AW^gu with the same sum as S, in which the antichain has 
units both at (ij,l) and at + 1). We replace S by this r-antichain, and set A = Ylje[ r ] ■ 

If Ak.i+i < Mk.i+i then either our m-antichain is already one in M' (if it does not use all units of M 
at (k,l + 1)), or it contains an antichain with a unit at (k, I + 1) and no units in column I, in which 
we can then replace its unit at {k,l + 1) by one at (k,l); in both cases we have found an m-antichain 
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of size n in M'. Now assume that Ak.i+i = Mk.i+i and let t = max{j | i'j = k} be the index of the 
last contribution to Akj+\- Definition 1.4.1b gives Y^lZij (-^i+i.M-i — ^f*,i) > 1 for all j < t, and the 
antichains . . . , give null contributions to the left hand side; therefore either M i ;+1 > A itl+ i holds 
for some i with ij < i < fc, or j > and > ij- Then we can set fc~ = max { i < k \ Aij + \ < M^j+i } 
and s = min{ j \ i'- > k~ } < t, and we have i s < k~ and ij < whenever s < j < t. 

If our m-antichain does not use all M k - i+1 available units at (k~,l + 1), then one can add such a 
unit to the antichain A^ to replace the one at (i' s ,l + 1), and then replace for all s < j < t the unit 
at (i'j, I + 1) in A^ by one at (i'j^iJ + 1) (some such replacements may give back the same antichain); 
thus our m-antichain is replaced by another m-antichain of size n in M that does not use all units 
of M at (k,l + 1), and is therefore also an m-antichain in M' . If our m-antichain does use all units 
at (fc~,Z + 1), then it contains at least one antichain A with such a unit and without any unit in 
column I. Any units of A in columns < I can be combined with the units of A^ in columns > I to 
from a new antichain B, while the units of A in columns > I can be combined with the units of A^ 
in columns < I; furthermore for s < j < t the units of A^ in columns < I can be combined with the 
units of A^- 1 ) in columns > I. Thus we have reorganised our m-antichain into another one of the same 
size in M, containing an antichain B with a unit at (k,l + 1) that can be replaced by one at (k,l). □ 

We are now set to prove theorem 5.2; the proof is just a formalisation of the argument given before 
our sequence of lemmas. 

Proof of theorem 5.2. If M, M 1 are a binary or integral k x I matrices, and M°, M'° are obtained from 
them by a rotation by a half turn (mapping the origin to (k — 1, 1 — 1)), then the relation M' = e\(M) 
with i e [k - 1] is equivalent to M'° = /fe_ 2 -i( M *)' wnilc M ' = e 7( M ) witn 3 e \ l ~ l \ is equivalent to 
M'° = f l ^ 2 _-(M' > ). Moreover, the rotation maps chains to chains and antichains to antichains. Then for 
each of the four preceding lemmas, one can derive the corresponding statement for the crystal operation 
in the opposite direction by applying the lemma to M°. Therefore each lemma may be strengthened 
to "M contains an m-(anti)chain of size n if and only if M' contains one" (lemma 5.5 speaks of an 
m-chain of size at least n, but this is only to make its proof easier to read; such an m-chain exists if 
and only if an m-chain of size n exists). Furthermore, applying rotation by a quarter turn to the two 
lemmas for binary matrices gives the corresponding statements with m-chains replaced by m-antichains, 
and with horizontal and vertical moves interchanged, while the two lemmas for integral matrices give 
the corresponding statements with horizontal and vertical moves interchanged by transposition. We 
may conclude both in the binary and the integral case that if M and M' are related by any crystal 
operation, then M contains an m-chain of size n if and only if M' does so, and the same for m-antichains. 
The theorem now follows from the fact that it holds for normal forms, as we had already observed. □ 
As the proofs above indicate, whenever M and M' are related by a crystal operation, one can 
obtain a to- ( ant i) chain in M' of maximal size by choosing an appropriate such TO-(anti)chain in M and 
either leaving its sum unchanged or moving one unit between the same squares as in the application of the 
crystal operation. Nevertheless, due to the fact that the moving unit may migrate from one (anti)chain to 
another and to the necessity of switching from one TO-(anti)chain to another with a different sum between 
successive crystal operations, the composition of the maximal TO-(anti)chains in M and the relations 
between them may be different from those in the normal form Diag [2I (A) or Diag(A) where A = ir(M). 
Notably the distribution of sizes of individual (anti)chains may be necessarily different: while in normal 
forms such TO-chains always decompose into chains of sizes A , . . . , A TO _i, an m-chain of maximal size for 
to > 1 need in general not even contain any individual chain of maximal size A . Here are examples of 
such a binary and integral matrix, both with A — (5,3,1), whose unique 2-chain of size 8 only decomposes 
into two chains of size 4, and fail to contain the unique chain of size 5: 

1 1 1 l\ /2 2 0\ 

1 £ M l2] and 1 e M (55) 

11011/ \0 2002/ 

(both examples were directly inspired by the one given in [Grcl, §3]). If one would trace the evolution 
of maximal m-chains from their normal form back to these matrices for to = 1 and to — 2, using the 
appropriate lemmas above, then it is clear that at some point the two cases must diverge; there is therefore 
no hope of improving those lemmas so that they keep track of complete decompositions into (anti)chains 
from which to- ( ant i) chains of appropriate sizes could be extracted for all to G N. The only thing one 
can say in general about the sizes of individual chains in a m-chain of maximal size is that they lie in 
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the range from A TO _i to Ao inclusive: they obviously cannot exceed the maximal size Ao of individual 
chains, and if any had size strictly less than A TO _i, then could remove that chain and be left with an 
(m — l)-chain of size exceeding X^e[m-i] wm ch is impossible. 

The property exhibited by the example allows another negative conclusion to be drawn: if Ai [2] and 
M are considered as partially ordered sets by 'C' (componentwise comparison) and V is considered so 
by inclusion of diagrams, then the implicit shape does not define a poset morphism. After all, the unit in 
the middle of the matrix is contained in any 1-chain of maximal size but not in all 2-chains of maximal 
size (in fact there is only one such m-chain in either case), so that suppression of this unit will decrease 
the initial part of tt(M) without decreasing the sum of its first two parts; consequently the diagram of 
7r(M) = (5, 3, 1) does not contain that of tt(M') = (4, 4) where M' is obtained from M by replacing the 
middle entry by 0. 

To prove that a given m-chain has maximal size it is not necessary to find all m-chains, but it suffices 
to produce any A m -antichain of maximal size. This is so because of the following general property, in 
which we define the "union" of two matrices asCliA = C + A — CnA = (max(Cij, Aij))ij e N'- 

5.7. Proposition. Let M be a binary or integral matrix, and A = ir(M). If for some k, I > and some 
sum Cofa k-chain in M and a sum A of an l-antichain in M one has \C<lA\ = kl and M = CU A, then the 
k-chain and the l-antichain are both of maximal size, and A& < I < Afc_i, in other words (k— 1, 1 — 1) G [A] 
and (k, I) ^ [A]. Conversely for any k, I satisfying Afe < I < Xk-i, the sum C of any k-chain of maximal 
size in M and the sum A of any l-antichain of maximal size in M satisfy \C fl A\ = kl and M = C U A. 

Proof. Since for any sum C of a fc-chain in M and a sum A' of an l-antichain in M one has C'flA' < kl, 
and therefore \C U A'| > \C'\ + \A'\ - kl, while on the other hand M D C U A' implies |C"UA'| < |M| = 
|C| + \A\ — kl, it is clear that in the first part |C| and \A\ are maximal among the sizes of fc-chains 
and of l-antichains in M, respectively. Now let C be the sum of any fc-chain in M of maximal size, so 
|C| = J2ie[m] an< ^ ^ et A be the sum of any l-antichain in M of maximal size, so \A\ = X^gm ^ji P ut 
d = kl — \ Cf)A\ and d' = \M — CUA\ (so in the case just considered one has d — d' — O). The expression 
giving \C\ can be interpreted as #[A]|[fe] X N (where [A]|x abbreviates [A] OX), the number of squares in 
the topmost k rows of [A]. Similarly the expression for |^4| can be interpreted as #[A]|nx[;]- Therefore 
the equation \C\ + \A\ = \C D A\ + \C U A\ can be written as 

#[A]| Wx n + #[A]|nx[I] =#([*] x [l])-d + #[X\-d'. (56) 

Now the only squares in N 2 with different contributions to the sets counted on the left and on the right are 
those in D = [k] x [i]\[A] or in D' = [A]|(isf_[/c]) x (N-m)) an d in both cases they contribute 1 to the right hand 
side and to the left hand side. Therefore d+d' = #D + #D', and we may conclude that d = d! = if and 
only if D = D' = 0, which is equivalent to (k — 1, 1 — 1) G [A] and (k, I) ^ [A]. The proposition follows. □ 

Note that although in the case of normal forms one even has d = #D and d' = #D' , this is not 
true in general: in the examples given the sum C of the 2-chain of maximal size is also the sum of a 
4-antichain, so in this case with (k, I) = (4, 2) and A — C one has d = but D = {(3, 1)}. 

While the proposition allows an easy certification of the maximality of a fc-chain and an l-antichain 
if an appropriate pair is given, it does not provide a procedure for obtaining such a pair. The best 
algorithmic method we know of to determine such witnesses for (partial information about) the implicit 
shape of a matrix M is to apply crystal operations to M until reaching a form in which such witnesses 
are evident (for instance a matrix encoding a straight tableau), and then trace back one's steps using 
the descriptions in the proofs of the lemmas of this section to find corresponding witnesses in M. This is 
somewhat laborious to do by hand, but still feasible for the matrices of our running examples. Here it is 
easiest to reduce to the form that encodes the straight tableau S of (31), namely is the matrix Q of (28) 
in the binary case, which is reached after 12 leftward moves, or the matrix P of (29) in the integral case, 
which is reached after 13 upward unit transfers. After this reduction we can find maximal fc-chains and 
l-antichains for every k and I by taking in the binary case the k topmost units (or all units if there are 
less than k) in every column of Q respectively the units in the I leftmost rows of Q, or in the integral case 
the k topmost rows of P respectively the I leftmost units (or all available units) in every row of P (in 
both cases these correspond to the k topmost rows and the I leftmost columns of the display of S). As 
an example, tracing that maximal 3-chain of size 13 and that maximal 4-antichain of size 16 back to M , 
we find in M the 3-chain coloured red, green and blue in the following displays 
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while the two 4-antichains one obtains in both cases can be described as follows by the sets of positions 
where their individual antichains have units: in the binary case these are {(0, 1), (1, 0), (2, 0), (5, 0)}, 
{(0, 3), (1, 1), (3, 1), (6, 1)}, {(1, 2), (2, 2), (4, 2), (6, 2)} and {(2, 5), (3, 3), (4, 3), (6, 3)}, and in the integral 
case these sets are {(1, 0), (2, 1), (3, 3), (4, 6)}, {(0, 0), (2, 1), (3, 5), (4, 6)}, {(1, 1), (2, 2), (3, 3), (4, 5)} and 
{(0, 2), (1, 3), (2, 4), (3, 6)}. The most complicated of these computations is for the 4-antichain in the 
integral case, where cases that require recombination of chains occur frequently; this explains the fact 
that the result for this case are rather tangled. There is also some liberty in interpreting the transitions 
prescribed in the proofs of the lemmas, since more than one antichain can contain a given unit; one could 
also perform some intermediate untangling of antichains, which would affect the final result as well. 

It is interesting to compare the results above with the original approach of [Grel]. There the roles 
of chains and antichains are played by increasing and decreasing subsequences extracted from a given 
finite sequence of numbers. In fact only sequences representing permutations are considered, but that 
restriction can be lifted if for chains one allows weakly increasing sequences while for antichains one 
requires strictly decreasing sequences; everything generalises immediately by treating equal entries as if 
they were increasing from left to right (a price to pay is that the theory becomes asymmetric, and the 
opposite convention would have been equally feasible, but we shall stick to this more usual convention). 
One can associate words to matrices in a way that is compatible with this interpretation of chains and 
antichains: for a binary matrix one takes concatenation of the strictly decreasing words whose weights 
are given by the columns of the matrix taken from left to right, while for an integral matrix one takes 
the concatenation of the weakly increasing words whose weights are given by the rows of the matrix 
from bottom to top. While this correspondence seems biased towards words that have this particular 
composition, one can in fact model words bijectively by matrices, by restricting the latter to those whose 
column sums (in the binary case) or row sums (in the integral case) are all 1 up to some point beyond 
which they are all 0: the mentioned concatenations then combine words consisting of a single letter. 

The invariance of the maximal size of m-chains or of to- antichains is then shown for a certain type 
of transposition of letters in the word, which generates an equivalence relation on words that is known by 
[Knu, Theorem 6] to characterise words with the same insertion tableau in the RSK correspondence. To 
be precise, the interchange of a pair of distinct adjacent letters is permitted in the presence of a letter I 
directly to the left or to the right of the pair, whose value is either strictly in between those of the letters 
being interchanged, or (in keeping with the interpretation of equal letters) if either I is to the left of the 
pair and its value equal to the larger of the letters being interchanged, or / is to the right of the pair 
and its value equal to the smaller one. If a, b, c designate letters with a < b < c, these transpositions 
can therefore be represented as transformations bac <-> bca, acb <-> cab, bab <-> bba or aba <-> baa of three 
consecutive letters; we shall call them elementary Knuth transformations. 

Modelling words by matrices, the fact that elementary Knuth transformations leave unchanged the 
RSK insertion tableau for the word is reflected by the fact that they can be realised by horizontal moves 
within the three columns (in the binary case) or vertical transfers between the three rows (in the integral 
case) corresponding to the letters concerned. We illustrate this in the binary case for transition bac — > bca: 

V Wf! Wf! Wfi 1 Wfi \, (58) 



V i/ V i/ \ i / \ i / V 1 

bac ba c ba c b ca bca 



and similarly for a transition aba — ► baa: 
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5 Characterisation of the implicit shape 

The word corresponding to the matrix changes only during the application of ff in the first case, and that 
of in the second. The integral case is similar; we show just these essential transitions for bac — > bca 
and for bab — > bba (keeping in mind that we must read from bottom to top): 

ac/1 1\ // a/1 \ , abfl IN f\ a/1 \ 

b[ 1 )^ bc{ 1 ij and b{ »( 2j- (60) 

It is clear that if these transformations are made inside a larger matrix M, then the one of the matrices 
P, Q associated to M in theorem 3.1.3 that is obtained by exhausting crystal operations in a direction 
parallel to those used here does not change; this matrix is Q in the binary case and P in the integral 
case. If M encodes a skew tableau T with at most one square in any row or column, then that matrix 
encodes both in the binary and the integral case the rectification S of T, which is also the RSK insertion 
tableau for the word read off from the display of T from bottom left to top right. 

It is less obvious that by repeating elementary Knuth transformations, one can transform the word 
corresponding to M into the word corresponding to the matrix that encodes S, which word can be read 
off from the display of 5* in the reverse Kanji (binary case) of reverse Semitic (integral case) order. The 
proof of [Knu, Theorem 6] is based on the fact that Schensted insertion can be emulated using elementary 
Knuth transformations: if T' = T <— x is the result of row-inserting an entry x into a semistandard Young 
tableau T, and t and t' are the words obtained by reading the displays of T and T' in reverse Semitic 
order, then t' can be obtained from tx by a sequence of elementary Knuth transformations (see also [Fult, 
§2.1]; the transpositions proceed in a single sweep from right to left). A similar result can easily be shown 
for column insertion x — ► T, reading the tableaux in reverse Kanji order. 

Since we have seen that elementary Knuth transformations can be realised by crystal operations, one 
may ask if one can conversely derive, from the exhaustion of crystal operations that transforms M into 
the matrix encoding S, a sequence of elementary Knuth transformations that realises the transformation 
of the corresponding words. Of course, since we remain inside one Knuth equivalence class, successive 
words obtained during the exhaustion of the crystal operations can be linked by elementary Knuth 
transformations, but the number of them required at each step turns out to be surprisingly large. It is 
shown in [Schii2, 2. 3. bis] how the change to the reverse Semitic reading of a tableau (with a hole) caused 
by the upward move of one entry occurring during an inward jeu de taquin slide can be realised by a 
sequence of elementary Knuth transformations (the goal of that statement was in fact to show that only 
the directions bca — > bac, bba — ► bab, acb — ► cafe, or aba — > baa are needed here; it should be noted that 
some transpositions might also allow an interpretation as a transformation in the opposite sense within 
another triplet). That case corresponds precisely to an upward unit transfer in an integral matrix. It is 
not quite clear how many transpositions are produced in the proof of the cited result, although it seems 
to be larger than necessary. We shall sharpen that result by providing a count of the transpositions used; 
experimental evidence strongly suggests that the number we give is in fact minimal. 

5.8. Proposition. Let M s M. allow an upward unit transfer from (io + 1, jo) to (io, jo) resulting in 
M' = e\ Q (M), and let the words w, w' be the concatenations of the weakly increasing words with weights 
given by the rows of M and of M', respectively, ordered from bottom to top. Then w can be transformed 
into w' by a succession of n 2 elementary Knuth transformations of type bca — > bac, bba — ► bab, acb — > cab, 
or aba baa, where n = J2j<j M i ,j + 52j>j M io+i,j- 

Proof. For brevity we take the designation bca — ► bac to include type bba — ► bab as well, and acb — > cab is 
similarly taken to include type aba — > baa. Put k = J2j<j ^io,j an d I = J2j>j -^io+iji so tnat ^ + ' = n - 
Write the rightmost n+1 letters of the subword corresponding to row io + 1 of M as c\, . . . , Ck, jo, di, • • • , di 
and the leftmost n letters of the subword corresponding to row i of M as oi, . . . , a^, b\, . . . , b t . The 
inequalities of definition 1.4.1a guarantee that the letters c x , . . . , c k and b\, . . . ,bi exist in those rows; 
the concatenation of these weekly increasing words forms a subword v = c\ ■ ■ ■ Ckjodi • ■ ■ d;ai • • • a^fei • • • bi 
of w in which all transpositions will take place, and we shall ignore all other letters. The same inequalities 
also imply that a, < Cj for 1 < i < k and fej < di for 1 < i < I, and if I > one has jo < fei; thus v 
corresponds to the reverse Semitic reading of a tableau at an intermediate stage of a jeu de taquin slide, 
when moving j upward in the situation 
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We must transformation v into v' = c\ • ■ ■ cud\ ■ ■ ■ dia\ ■ ■ ■ akjabi which we shall do explicitly by the 

following choreography. We proceed in n sweeps of n elementary Knuth transformations each, starting 
the successive sweeps with the triplets whose rightmost letters are successively a\, . . . , a^, b\, . . . , bi, and 
proceeding from right to left within each sweep. During sweep i with 1 < i < k, the letter is first shifted 
to the left n — i + l places by transformations of type bca — > bac, after which it ends up to the right of the 
letter a; that letter is then shifted to the left i — 1 places across the letters a\, . . . , dj_i by transformations 
of type acb — > cab, which results in a word c\ • ■ • Cjai • • • a^Ci+i • • • Ckjodi • • ■ diai+i ■ ■ ■ akb\ After 
these k sweeps one therefore has the word ci ■ • ■ Cfcai • ■ ■ akjodi ■ ■ ■ dib\ ■ ■ ■ bi. During sweep k + i with 
1 <i < I, the letter bi is first shifted to the left l—i places across the letters di+i, . . . , di by transformations 
of type bca — > bac in which it has the role of a, after which it finds itself to the right of the letter df, 
that letter is then shifted to the left k + i places by transformations of type acb — > cab, which results in 
a word c\ ■ ■ ■ c^d\ ■ ■ ■ rf»ai • • • akjobi ■ ■ ■ bid i+ i ■ ■ ■ dib i+ i • • • 6/. After these sweeps we have the word v'. □ 

For a leftward move in a binary matrix from (in, jo + 1) to (io,jo), there is a similar realisation of 
the change in the corresponding words, which is the change in the reverse Kanji reading of a tableau that 
occurs during a leftward step in a jeu de taquin slide, using n 2 elementary Knuth transformations of the 
same types, where n = J2i<i Mij + J2i>i ^ijo+i- The proof is the same, interchanging left and right 
in the word, and with the implied interchange of weak and strict inequalities. 

The proposition indicates that the number of elementary Knuth transformations corresponding to 
the exhaustion of upward crystal operations can be considerably larger than the number of such opera- 
tions; in our example we needed 13 upward unit transfers, which correspond to 150 elementary Knuth 
transformations. This can be improved a bit when there are successive transfers between the same pair 
of rows: if there are k such transfers then all involve the same value of n, and it appears that by an 
adaptation of the transposition procedure one can achieve the combined transformation of the word in 
n(n + k — 1) elementary Knuth transformations instead of kn 2 ; in the example this allows the total number 
of steps to be reduced to 122. This is still worse than the number needed for the traditional emulation of 
the RSK insertion, although it is of the same order of magnitude: a straightforward application of that 
procedure to the word corresponding to M gives a sequence of 120 transpositions, which can be reduced 
to 88 by observing that the first 32 transpositions bring back the word to its original form, due to the 
circumstance that the word corresponding to the three bottom rows of M is already the reverse Semitic 
reading of a straight tableau. It should be noted that our matrix M was chosen as the integral encoding 
of a skew semistandard tableau that is fairly close to its rectification, whence the numbers above are 
relatively small. These data suggest that for the purpose of finding maximal m- ( ant i) chains in a binary 
or integral matrix, there is some advantage in using the transformations described in our lemmas rather 
than the (somewhat simpler) ones corresponding to elementary Knuth transformations. 

Another approach to finding maximal m-(anti)chains would be to try to augment the inductive 
computation of n(M) via growth diagrams with additional data, keeping track of maximal to- ( ant i) chains 
as witnesses for the implicit shape of submatrices. It appears however that this can only be done for m = 1. 
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§6. Remaining proofs. 

We shall finally provide, for completeness, proofs of statements for which they were so far omitted. The 
first such statement was the fact that crystal operations can be defined in such a way for any semistandard 
tableau T that they correspond both to vertical crystal operations on its binary encoding M, and to 
horizontal crystal operations on it integral encoding N. 

Proof of proposition 2.1. We know that if cither eJ„(M) or e^(7V) is defined, then like M or TV it 
satisfies the tableau condition for A/ n (propositions 1.3.6 and 1.4.6), and is therefore the binary or integral 
encoding of some tableau T". Moreover T' can easily be seen to differ from T only in the partition \( m+1 \ 
for which some square (k, I) has been added to its diagram to form the diagram of the corresponding 
partition in X" (in terms of the display of T, an entry m + 1 in that square has been changed to m). We 
have to show that if the other type of encoding is taken of T and of T' , then the matrices obtained are 
also related by the appropriate raising operation. By interchanging the roles of T and of T' , this will also 
prove the statement about fa and fa. 

Suppose first that the binary encodings of T and T' are M and e| n (M), and let their integral encodings 
be N,N', respectively. Then with k,l as above, we know that the bits M m j — and M m+ u = 1 are 
interchangeable in M, and we must show that a unit transfer is possible from Nk, m +i to iVfc ;TO in N. 
Consider the left member ^i=fc 1 (-^ r i,m — A^+i.m+i) of the second inequality in definition 1.4.1b, which 
we must show to be nonnegative for any k' > k. It counts the number of squares of the skew diagram 
[A( m+1 )/A( m )] that lie in rows i G [k'] — [k], minus the number of squares of [\( m + 2 ) / \( m+1 )] in rows 
i + 1 G [k' + 1] — [k + 1]. Now these squares can also be described as those that lie in columns j e [I] — [I'] 
where V = Al™ +1 ^ < I = A^ m+1 \ since squares of [A^ m+1 ^/A^ m -'] in row k are in columns < I while 
squares of [A*- m+2 )/A^ m+1 ^] in row k are in columns > I, and similarly for k' and I'. It follows that 
J2te[k']-[k]( N hm - N l+1 . m+ i) = J2je[i]-[i']( M m,j ~ M m+i,j) which is nonnegative by hypothesis. For the 
first inequality in definition 1.4.1b the argument is similar. Here it is easiest to take the version of that 
inequality for N', which is equivalent to the one for TV: for any k' < k one puts V — A^T +1 ' (which 
satisfies V > I + 1), and one shows that T,ie[k]-[k']( N i+i,m+i ~ N i,m) = Ej 6 [l']-[l+i] ( M ™+i,j - M mJ ), 
which is nonnegative. This concludes the proof that N' = e^(7V). 

Now suppose that the integral encodings of T and T" are N and e^N), and let their binary 
encodings be M,M', respectively. The argument we have just given shows that, among the inequalities 
in definition 1.3.1a that must be shown to hold in order to justify M' — eJ rl (M), those for which 
I' e {A^ m+1) k' > k} U {1} U {A<™ +1) -1 I k' < k} are satisfied. This does not cover all inequalities, 
but the values in rows m, m + 1 of M that are incorporated into the summations on j between such 
values of I' have a simple structure, as they account for squares of [A( m+2 )/A' m+1 ^] in row k', and 
squares of [A^ m+1 ^/A( m '] in row kl — 1, for some fixed k' depending on the interval considered. In such 
an interval one has from left to right: a possibly empty sequence of columns ( 1 ), a possibly empty 
sequence either of columns Q) or (J), and a possibly empty sequence of columns (J); similarly the 
columns incorporated into the partial sums after the largest such value of V can only be (°) or (°). It 
follows easily from this description that the required inequalities also hold at the remaining values of V . □ 

The next statement left to prove is that the involutions defined by crystal operations on matrices 
satisfy the Coxeter relations of type A. We did not use this fact anywhere, but it does represent an 
interesting aspect of crystal operations, so we present the proof, even if this turns out to require quite a 
bit of work. 

Proof of lemma 2.3. This is a statement about the structure of crystals, so by proposition 2.1 it suffices 
to consider the involutions a" in the case of integral matrices. We have already seen that the are 
involutions (the first Coxeter relation), and the fact that cr~ and a" commute when \l—j\ > 2 (the second 
Coxeter relation) follows from the fact that they operate on disjoint pairs of rows. Hence the only non- 
obvious point is the third Coxeter relation, but here a direct proof would be difficult, because for instance 
the effect of applying e\ to M on the operation of e[ +1 depends on M in a complicated way. Instead we 
shall proceed by reducing the general case to one where a detailed analysis is feasible. In general, if M' 
is obtained from M by applying an invertible operation that commutes with erf and with crf +1 , then the 
equation afaf +1 af(M) = af +1 crf af +1 {M) holds if and only if afaf_ 1 af(M') = af +1 crf af +1 {M') does; 
we shall use such operations to reduce M to a simpler form. 

First of all we replace M by M Nx ^i +3 ^_^ (the inverse operation of adding back the columns outside 
the range {1,1 + 1,1 + 2} commutes with <rf and <7^ 1; which only look at the columns in that range), 
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and we may then as well assume that k — 0. Secondly we may assume that upward crystal operations 
applied to M, which commute with af and af_ x , are exhausted. Thus we have reduced M to the form 

(a b c\ 
d e G M with a>d>/anda + 6>d + e. (62) 
f) 

Although the proposition can be directly proved for this situation, as we shall indicate below, one can take 
a shortcut, as in the proof in [Loth], by further reducing the situation using the fact that the operations 
af commute with cyclic permutations of the nonzero rows of the matrix, as stated in lemma 6.1 below. 
Doing so to the matrix M of (62), moving row below the two following ones, one can again exhaust 
upwards crystal operations. In fact using e\ all a + b + c units that were cycled to row 2 can be transferred 
up to row 1, and then can be applied a + b + c — mm(d, c + f) times to obtain a matrix 

(a b + d c'\ 
e' where e' = min(d, c + /) and c' = c + f - e' + e (63) 

/ 

in which upwards transfers are once again exhausted. Thus the case of matrices of the form (62) has been 
reduced to the case of matrices of the form (63). The latter case is contained in the former, and we can 
repeat the steps the we applied to further reduce to the case of a one-line matrix M" = (a b + d c+e + f). 
Now for such matrices one has col(M") = Mq , and since we know that the effect of af on col(M) is to 
transpose the parts i and i + it is clear that af af af \M") = af afaf(M") = (c + e + f b + d a). □ 

We shall now state formally the commutation of operations af with cyclic permutations of rows that 
we used. It is valid for binary matrices as well as for integral ones, and it also has an analogue for af; 
since the statement is of independent interest we state it in its most general form. 

6.1. Lemma. Let n, m G N, d G [n] and d' G [m]. Denote by R either the set of matrices M G A4 [2] 
or of matrices M 6 M, such that M — Mr n i x r m i, and by ir r and ir c operations on R defined by the 
cyclic permutation of the rows { Mi \ i G [n] } by d places (so ~K T (M)i^ = Mi +C i mo d n _j for all M G R, 
(i,j) G [n] x [m]), respectively by the cyclic permutation of the columns { M| | j G [m] } by d! places (so 
TT c (M)ij = Mi j + d> modmj- Then for all k,l G N and M G R one has af(n Y (M)) — ir r (af (M)) and 
4(n c (M)) = n c (ai(M)). 

Proof. By symmetry it will suffice to consider the statement concerning af. Clearly only the columns 
1,1 + 1 that are affected by af are of any concern, so we may assume they are the only ones: I = 
and m — 2. Note that the effect of 7r r and af on row and column sums is determined without knowing 
more details of M, and these effects commute, so our only worry is that both sides of the equation 
might give different matrices with the same row and column sums. We shall start considering the 
special case of integral 2x2 matrices (so n — 2 as well) with d = 1, so that 7r r is the interchange of 
the two rows. We put row(M) = (s,t) = row(af(Mj), then M has the form M = (° s ff a b ). One 
has 5 — col(M)i — col(M) = s + t — 2(a + b) in (25), and the same value applies to tt i (M). By 
symmetry between M and w T (M) we may assume that s < t; wc distinguish the cases with s < a + b < t 
from those where either a + b < s or t < a + b. The condition s < a + b < t means on one hand 
that min(t — b,t — b — S) > a, so that in this case the transfer for computing af(M) takes place 
entirely in the bottom row, and on the other hand that s — a < min(6, b — S) so that the transfer for 
computing cr^(7r r (M)) takes place entirely in the top row of 7r r (M) = ( *Z ); as a consequence one has 
af(w r (M)) = ( b ~ s *7_ b +' 5 ) = ir r (af(M)). If instead one has a + b < s < t, then it follows that S > and 
in fact t-b > a > t-b-6, whence we may write af(M) = (ef) s (M) = (ef) s - {a+ ^ ((ef) t -^ a+b ^> (M)), in 
which expression (e5~)'~( a+b - ) operates in the bottom row, after which the factor (ef ) s ~( a+b ) operates in 
the top row. The computation of Oo*(7r r (M)) and can be performed similarly using s — a > b > s — a — <5; 
the two computations can be summarised as follows: 

M (a s-a\ (i) H " +l) / a s - a \ (efy-^f s - b b\ 

M= {b t-b) ^ [t-a a ) ^ [t-a a) =<J ^ M ^ 

(M\ ( b t-b\(efy- (a+b> ( b t-b\(efY-^ (t-a a\ ''^ 
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6 Remaining proofs 

Clearly (n r (M)) — tt t (<Jq (M)) in this case. In the remaining case t < a + b one has, with a' = t — b 
and b' = s - a, that a' + b' < s and M = (*~^ *',), so in fact M = o-jf(M') for a matrix M' = *~*') 
to which the previous case applies, which settles this final case. 

Now we return to the more general case of <7q* operating on a n X 2 matrix M (either binary or 
integral). Write M = (|j where S is dx2 matrix and T a (n — d) x 2 matrix, so that 7r r (M) = (^). Any 
horizontal crystal operation in M occurring inside S or T will also be valid in that part considered in 
isolation. But definition 1.3.1b or 1.4.1b implies that if conversely S — ► S' is a horizontal crystal operation 
that is valid in isolation, then — > (f,) is valid in M only if min(n^*(S'), n^*(S")) > n^{T), while such 
an operation T — » T' gives rise to one (^) — * (^,) only if n^* (5) < min(n5~ (T) , nj" (T')) . This means that 
the moves or transfers involved in the transition M — + <Jq{M) will be distributed between the parts S 
and T in the same way as for M = j^b) £ -M the transfers involved in the transition M — * aff(M) are 
distributed between the two rows, where a = n c f(S l ), s = a + nff (5*) , b = n${T), and t = 6+n^(T). So if 
ct ~(M) = il^ ), then one has crJ*(M) = (^), where either 5' = (e^) 5 " (S) and T' = (eo) Sl {T) 

with 5 ,(5i > 0, or S' = {fn*)- 5o (S) and T' = (Jq)~ &1 {T) with 5 ,*i < 0. We have seen that in all 
cases one then also has cr^( a *Z„) = (o+5o s-a-5o)' an< ^ reason i n g for 7T r (M) = (^) in the same way as 
for M, one finds cr^(^) = (^,) with the same values S',T', which proves o-^(tt t (M)) = n T (a^(M)). □ 

We note that in the basic case of 2 x 2 integral matrices, the identity <7Q*(7r r (M)) = n r (a^ (M)) comes 
about in a somewhat curious way when a + b < s or t < a + b, since the intermediate matrices traversed 
during the computation of erj* in both members no not correspond to each other under 7r r . The situation 
can be illustrated by drawing for fixed row sums (s, t) individual matrices as dots with coordinates given 
by their first column (a, b) € [s + 1] X [t + 1], and drawing both the ladders for e$ on one hand, and 
ladders for the conjugate operation 7r r o o 7r r on the other hand. This is done below for (s, t) = (5, 9), 
with the first type of ladders in red on the left and the conjugated ladders in blue on the right. We have 
also marked three pairs of matrices that are interchanged by crj* (as well as by 7r r o <rj* o 7r r ). One sees 
that outside the range given by s < a + b < t (delimited by the diagonal lines) the pairs of matrices 
exchanged form the only common points of a pair of ladders of the two types. 




As a consequence, the simplification in the proof of lemma 2.3 that we obtained by cyclicly permuting 
rows, while commuting with the operations crj* , does not preserve the structure of the crystal. Indeed in 
the case of one-line matrices ultimately reduced to, horizontal crystal operations between different pairs 
of columns commute among each other in the weak sense that e^(e^(M)) = e$ (ej~ (M)) whenever both 
members are defined; such a relation does not hold in general crystals. 

It is therefore instructive to present as an alternative a direct way to complete the proof of lemma 2.3, 
without using lemma 6.1: thus one will only have used reductions that preserve the structure of the crystal. 
This also provides a less indirect proof for the binary case: we used the isomorphism of the binary and 
integral crystals provided by proposition 2.1, but the operation of cyclicly permuting rows of an integral 
matrix is very hard to translate across this isomorphism (since it breaks out of the class of matrices 
corresponding to tableaux of a given shape). By contrast, the direct computations below, while given for 
integral matrices, translate easily into a computation for binary matrices. It suffices to observe that the 
matrix M of (62) corresponds, by the intermediary of a semistandard Young tableau with three horizontal 
strips, to a binary matrix B with three rows, of which the columns Bj can be described by associating 
to them three bit binary numbers Cj — 4i? j- + 2-Bij + B 2 j 6 {0, 1, . . . , 7}, which then form a weakly 
decreasing sequence, and whose frequencies m, = #{ j S N | Cj = i } can be deduced from the equations 
m 7 + m 6 + m 5 + 7714 = a, to 3 + m 2 = b, mi = c, m 7 + m 6 = d, m 5 + m 3 = e, m 7 = /, and 77147713 = 
(which holds since columns B* = (1,0,0) and Bj, — (0, 1, 1) cannot both be present). 
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6 Remaining proofs 



Before we do the computation, we make one more simplification: from M we can subtract / times the 
3x3 identity matrix, since / is the smallest diagonal entry, and the operation is easily seen to commute 



with the operations f 



and /f*(in terms of binary matrices this corresponds to removing the 
m 7 initial columns £?' = (1, 1, 1), which play a completely passive role). Thus we can assume / = 0, and 

write M={l b d 



G M. 



Since a> d, all units of the coefficient Mi i = d in M are fixed for the operations ej", whence 



a b c 
Ode 



b + d a — d c 
d e 



(65) 



Computing an expression for erf*(M), which will be of the form d , c e , ), is a bit more complicated. 
Assume first that min(6, e) = 0, then one will also have min(6', e') = 0; from the fact that moreover 
( b d , c e , ) must have rows sums (b + c,d + e) and column sums (c + e, b + d), one can deduce the expression 
c' = min(fe + c, b + d) = fe + min(c, d) for this case. We can reduce the general case to this one by removing 
min(6, e) units both from b and from e, which should leave c' unchanged; it follows that one must replace 
b in the expression for c' by 6 — min(6, e) for the general case. One then finds c' = b — min(6, e) +min(c, d), 
b' = c + min(6, e) — min(c, d), d' = e — min(6, e) + min(c, d) and e' = d + min(fe, e) — min(c, d), in other 
words 

a b c\ _ ( a c + min(6, e) — min(c, d) b — min(6, e) + min(c, d) 
d e) \0 e — min(6, e) + min(c, d) d + min(6, e) — min(c, d) 



(66) 



Now the verification that agcr^a^ (M) = o~i o~q (M) is a straightforward but laborious computation: 
one determines the appropriate expressions in a, b, c, d, e, which are piecewise linear ones built using the 
operators '+', '— ', and 'min', and one shows that the corresponding parameters of (Tq cr$* (Tq (M) and 
Oi <j± (M) are given by equivalent expressions. The relevant expressions are given in table 1. 



M 


/ a b c\ 
\0 d e) 


o-o(M) 






(b + d a — d 
\ d 


:) 




ara "(M) 


(b + d 

\ o 


c + min(a — d, e) — min(c, d) 
e — min(a — d, e) + min(c, d) 


a — d — min(a 
d + min(a 


— d, e) + min(c, d) \ 

— d, e) — min(c, d) J 




(c+e 

\ o 


b + min(a — e, d) — min(c, d) 
e — min(a — d, e) + min(c, d) 


a — d — min(a 
d + min(a 


— d, e) + min(c, d) \ 

— d,e) — min(c, d) J 


ar(M) 




(! 


c + min(6, e) — min(c, d) 
e — min(&, e) + min(c, d) 


b — min(6, e) 4 
d + min(6, e) - 


- min(c, d) \ 

- min(c, d) J 


<7 ~ar(M) 


c 


+ e 



a — e + min(6, e) — min(c, d) 
e — min(6, e) + min(c, d) 


b — min(6, e) + min(c, d) \ 
d + min(6, e) — min(c, d) J 


o-r cr o* <T r( M ) 


f c+e 


b + 
d- 


min(a — e, d) — min(c, d) 
min(a — e,d) + min(c, d) 


a — e — min(a 
e + min(a 


— e, d) + min(c, d) \ 

— e, d) — min(c, d) J 



Table 1. Verification of a^a^a^(M) = a^cr^a^(M) for M = 



a b 
d 



Since all row and column sums are known beforehand, only one coefficient not in the leftmost column 
really needs to be computed at each step, but we have retained for each coefficient a form close to the one 
obtained from application of the defining formula, possibly after some simplification. The most important 
simplifications take place when (66) is applied to compute (a^ (M)); for instance the coefficient at 



59 



6 Remaining proofs 



index (0, 1) is computed as c' + min(6', e') — min(c', d!) where ( 



a' b' c 
d> e> 



/) = cr (T cr r( M )> giving 



b — min(6, e) + min(c, d) 



+ min(a — e + min(6, e) — min(c, d) 



d + min(6, e) — min(c, d)) 
min(6, e) + min(c, rf)) 



— min(6 — min(6, e) + min(c, d) , e 
= b — min(6, e) + min(c, d) 



+ (min(6, e) — min(c, d)) + min(a 
= b + min(a — e, d) — min(c, d). 



e, d) + (min(6, e) — min(c, d)) 



— min(6, e) 



Since the result is identical to the corresponding coefficient of <Jq a~ <Jq (M) , all other coefficients must be 
equivalent as well, which is also easily checked explicitly. This computation demonstrates the potential 
complications in composing piecewise linear expressions and deciding their equivalence; here the exercise 
remained doable because there is effectively just one component defined in a piecewise manner, with 
the other components depending linearly on it. The general situation in the statement of lemma 2.3 
for integral matrices is much more complicated than that of (62), so exhausting operations e[ during 
the initial simplification, while preserving the structure of the crystal, does greatly simplify the problem 
viewed as an equivalence of piecewise linear expressions. 

There is one more proof that we have omitted, namely that of theorem 4.3.2, which states that for 
binary matrices the computation of implicit shapes using French normalisation leads to the row-insertion 
shape datum for (V, t— , [2]). We did suggest that the assiduous reader could check this, but we did 
not prove her existence. So here is a proof; those who already found one may skip this. 

Proof of theorem 4.3.2. We return to the French normalisation process described in the display (43), 
and analyse the exhaustion of { // \ i £ [k] } applied to A, on the left. To the left of the vertical line, the 
submatrix A[ fe+1 ] x [;j, consisting of NF(A, k) with below it v l — A' in row k, is transformed into NF(f, fc+1), 
moving down one place every unit in a column j with z/| — A' = 0. The exhaustion of the // proceeds by 
weakly decreasing i, and after operations with i < i' < k have been exhausted, row i + 1 will contain 
a unit in a column j < I if and only if vj — Aj = 1 and vj > k — i. Such a unit could block a potential 
downward move in column I if in addition the bit directly above it is A^j — 0; such a pair (J) will exist 
if v* = k — i and Aj = k — i — 1, in other words when (k — i — 1, j) £ [^/A]. Indeed the existence of at 
least one such pair will rule out any downward move between rows + 1 in column /, since there can 
be no pairs (J) to the left of column I that could be blocked instead. In the example of (43) we see this 
phenomenon for i = 1: since [^/A] has squares in row k — i — 1 = 3 (namely in columns and 1), the 
unit at (6, 1) is blocked, although the unit at (6, 2) has been moved down by before f[ is applied. 

We now consider column / to see how its initial value of revfe(^ — A) followed by a bit Akj = M^j is 
transformed into revfc + i(K — v), thereby determining n. If after exhausting /{ the bit at (k — i,l) is '1', 
then we shall say there is a carry into row i; initially there is a carry into row if and only if Mk y i = 1. 
Now if there is a carry into row i, then [k/v\ will certainly have a square in row i, and if moreover [/x/A] 
already had a square in row i, there will be a carry into row i + l. If there is no carry into row i and [///A] 
has no square in row i, then neither will [k/u], nor will there be a carry into row i + l (here column / 
contains ([j)). If there is no carry into row i but [fi/X] does have square in row i (there is a pair (*) in 
column I), then that square is also in [k/v] if vi = Xi (the unit in column I moves down), or else it gives 
rise to a carry into row i + l (if [v/\] has any squares in row i, so that the pair in column / is blocked). 

Now comparing with the description of the row-insertion shape datum for (V, [2]) given before 

theorem 4.3.2, any square G [a*/A] fl [V/A] is an element of S that is absent from [A], and it assures a 
carry into row i + l; if there was also a carry into row i, then the square of [k/v] in row i is some (i,f) G T 
with j' > j. If there is a carry into row i but no square of [/it/ A], then (i, Vi) G T is a square of [k/v], and 
conversely a square G T can only lie in [k] if there is a carry into row i, because if there is not, any 
square of [k/v] in row i is also in [/x/A] and therefore not in T. Thus k is as given by the shape datum. □ 

Alternatively one could determine k by analysing what happens during the French normalisation of 
j4[fc]x[*+i]- This is interesting and proceeds rather differently than in the given proof, but this time we 
really leave this for the enthusiasts to discover for themselves. 
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Conclusion 



Conclusion. 

After introducing crystal operations on binary and integral matrices, we have performed constructions 
with them related to various known facets of the RSK correspondence. We have found links with all of 
the facets considered, and in many cases obvious questions suggested by the results found have led us to 
investigate in detail the precise relations with the existing theory. This required some meticulous work, 
but the result has shed new light on many aspects of that theory, and has resulted in a framework in 
which those aspects are more tightly and explicitly integrated with each other than they were before. In 
fact we have been surprised by the extraordinary versatility of these crystal operations, whose original 
humble vocation was to streamline the proof of a version of the Littlewood-Richardson rule found in 
[vLee5], and by the dimensions this paper has taken as a consequence. To close, we mention one more 
obvious question that we shall not answer here, namely whether for the notions of crystals in other types 
than A n there exists a similar notion of double crystals. 
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